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HIGHER-ORDER FOCUSING IN THE 17/2 8-SPECTROMETER'! 


G. E. LEE-WHITING AND E. A. TAYLOR 


ABSTRACT 


The expansion coefficients of that magnetic field which gives best focusing in 
the high-aperture double-focusing spectrometer have been determined to the 
sixth order. The analytical calculations have been checked by numerical inte- 
gration of the equations of the orbits. ‘The third- and fourth-order coefficients 
published by Stoker et al. (1954) are not in agreement with our results. 


1. INTRODUCTION 


It seems to be admitted universally that the double-focusing principle 
introduced into 8-spectrometry by Svartholm and Siegbahn (1946) is most 
advantageous. The flat field of the original 180° type is replaced by an inhomo- 
geneous field having the correct slope at the central orbit to produce first- 
order radial focusing after an angle of r+/2 radians. For this particular value 
of the slope, a = —} in the notation of Siegbahn (1955), first-order axial 
focusing also occurs, bringing the great benefit of a much reduced image 
height. The question of the best value of the coefficient of the quadratic term 
in the expansion of the magnetic field, 6 in the notation of Siegbahn, is not so 
easily answered. The choices 6 = }, 3, and } have been suggested. We shall 
show that both 4 and # are superior to } when a rectangular aperture is used, 
except possibly at very high transmission. The decision between $ and ? is not 
so clear-cut, though there are several minor advantages to be gained by 
adopting the latter. For either choice, but for no other, it is possible to bring 
about a great improvement in transmission by employing higher-order 
focusing. Higher-order focusing requires a more accurate control of the field 
shape. Though the accuracy required may not be possible with an iron-cored 
instrument, there is no obvious reason why focusing to say the fourth order 
should not be possible with air-cored coils. In addition, iron-free spectrometers 
have other well-known advantages. Moussa and Bellicard (1954) and Siegbahn 
and Edvarson (1956) have succeeded in building double-focusing spectrom- 
eters with different types of air-cored coils. 

A spectrometer with higher-order focusing based on the choice 6 = } will 
have a much reduced radial aberration for thos: orbits lying near the plane 
of symmetry. The theory of such spectrometers has been given by many 


1Manuscript received September 11, 1956. 
Contribution from Chalk River Laboratory, Atomic Energy of Canada Limited, Chalk 
River, Ontario. Issued as A.E.C.L. No. 873. 
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authors; we give references in Section 3. A spectrometer with higher-order 
focusing based on the choice 8= } favors orbits lying near the right cylinder 
standing on the central orbit. Our main purpose in writing this paper is to give 
the optimum field parameters as far as the sixth order. Though it may not be 
possible to fit the field beyond the fourth order, the fifth- and sixth-order 
terms should at least be useful in calculating the residual aberrations. Our 
parameters, which do not agree with the fourth-order results of Stoker et al. 
(1954), are confirmed by numerical integration of the equations of the orbits. 
The optimum field form is very close to r~, the shape approximated by the 
solenoids of Siegbahn and Edvarson (1956). 

In undertaking the complex problem of a higher-order orbit calculation it 
is well to begin with a systematic treatment of the properties of the various 
expansion coefficients of the magnetic field. Such a discussion is given in 
Section 2. The first- and second-order properties of the spectrometer :zre 
reviewed in some detail in Section 3 in order to permit us to give a reasonable 
discussion of the advantages of the different choices for 8. The results of the 
higher-order calculations are presented in Sections 4 and 5. 

Following Svartholm (1946) we eliminate the time from the equations of 
motion and obtain differential equations in the spatial variables for the electron 
orbits. These differential equations are solved by a method of successive 
approximations. 


2. DESCRIPTION OF THE MAGNETIC FIELD 


We restrict our discussion to magnetic fields having an axis of cylindrical 
symmetry perpendicular to a plane of symmetry. Assume a cylindrical co- 
ordinate system (r, ¢, 2) with its z-axis in the axis of symmetry and its origin 
in the plane of symmetry. For the magnetic fields under discussion only the 
$-component of the vector-potential is non-zero; call this component A. The 
radial and axial components of the magnetic field are then 


0A 

(1) G.(r,2) = -F,, 
1 OrA 

2) G.(r,2) = 5 
In a shaped-field 6-spectrometer the source and detector will be located on a 
particular circle in the plane of symmetry centered at the origin. We call this 
circle the optic circle. Its radius, a fixed parameter of the spectrometer, we 
represent by the symbol 7o. Stationary motion on the optic circle can exist 


for electrons of a given wholly tangential momentum, say po, if the field at 
the optic circle, Ho, is adjusted to satisfy the relation 


(3) H.(ro, 0) = Ho = —Polero)-. 
We shall use differential equations for the electron orbits in the form given 


by Svartholm (1946). The properties of the magnetic field enter these equations 
only through the product rA. A method of finding the coefficients in a certain 
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expansion of this product is the chief business of this section. It is convenient 
to transform from the coordinates (r, 2, @) to the coordinates (n, 7, ¥) defined 


by 

(4) 2 =(r—1)/r, t=2/r, w= o/V72. 

We write the vector-potential in the neighborhood of the optic circle in the form 
(5) rA= roA ot+ Doro’ = mat” n ’ 


in which ¢o99 = 0. Since z = 0 is a plane of symmetry, the coefficients with m 
odd are necessarily zero. We demand that ¢o; = 1, in order that the field on 
the optic circle derived from (5) be consistent with stationary motion on the 
optic circle of electrons with tangential momentum pp, i.e. in order that (3) be 
satisfied. The differential equations for the orbits depend on the initial con- 
ditions. Let us consider an orbit for which 7 = h, r = t, n' = H, and 7’ = T 
at y = 0. The angle between the initial momentum vector and the tangent 
to the optic circle at y = 0, , is given by the relation 


(6) cos y = [1+3(1+4)7(#?+7%}7. 


Let the absolute value of the momentum of the electron in question be 
Pp = po(1+e). We shall need the function F defined by 


(7) F = (1+e)(1+h)cos y+ >> Cma(r"n"—t"h" ), 
and its derivatives as defined by 
oo wie 
(8) P = on ’ Q = Or ¥ 
The differential equations for the orbits are 
n” P 
(9) +2 1+ 2r— 2 aly = aa ey 


—2(1+n)* 5 g. 


(10) v4 2 +2 — |e 


For e = 0 these equations are identical with Svartholm’s equations (18). 

The coefficients Cm, introduced in (5) are not all independent. Of the two 
Maxwell equations applicable to a static magnetic field in vacuo, one is auto- 
matically satisfied by a field derived from a vector-potential. Under the sym- 
metry here imposed the other reduces to 

99, _ 99: 

By substituting (1) and (2) derived from (5) in (11) and equating coefficients 
of 7" from the two sides of the equation, we obtain 


(12) (m+2)(m+1)[¢m+2 n+1 + Cm+2 n] = —(n+2)[(n+3)cm n+3 + Cm n+2] 


valid for n > 0, and 
(13) (m+2) (m+1)Cm+2 = —2m2+Cmi- 
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It follows from (12) and (13) that, given Cn for fixed m and all n, we can 
compute ¢Cm+2 n for all 2. For, after using relation (13) to find cm+2 9 in terms of 
known quantities, one can employ (12) as a recurrence relation to increase 
the value of ” in unit steps. Thus we have shown that the matrix of coefficients 
Cmn is completely determined when we fix the row Con. Since Maxwell’s 
equations have been satisfied, the coefficients Co, are independent. 

The variation with r of the field intensity in the equatorial plane, i.e. the 
function §,(r, 0), is customarily used to describe the magnetic field of a 
B-spectrometer. To show that a knowledge of this function is equivalent to a 
knowledge of the set of coefficients con, we expand the function about the 


point r = 7: 
(14) $.(r, 0) = God a, 1". 


In order that (14) should satisfy (3), we must make the definition a) = 1. 
(The coefficients a, defined by (14) are identical with those used by Stoker et al. 
(1954). The a and 8 of Siegbahn (1955) are a; and a2 respectively; the a and 8 
of Shull and Dennison (1947) are —a, and az respectively. The coefficient a3 
used by Svartholm is our ¢o3.) By equating expressions (2) and (14) for §,(r, 0) 
it is not difficult to show that 

(15) Non = An—-1+An-2 


for > 2; the relation is true for 2. = 1 too, if we make the definition a_; = 0. 
From (15) the set co, may be calculated if we are given the set an. Conversely, 
since @) = 1, we may use (15) as a recurrence relation to generate the a, from 
the con. Hence the properties of the magnetic field in the neighborhood of the 
optic circle are completely and uniquely defined by either set of independent 
coefficients, the con or the dn. 

There are several relations between the cm» and the a, which are useful in 
checking results obtained by the equations already given. In cylindrical co- 


~ 
ordinates the equation div § = 0 becomes 


OO: = 
(16) ~ = VG,)+ 
Differentiation of (16) with respect to z and the use of (11) leads to 
(17) =(A-+~ ) FS: _ +(1+n ) 28 


We shall use this relation in the special case r = 0. By substituting the §, 
obtained from (2) in the left-hand side of (17) and that of (14) in the right-hand 


side we obtain 

(18) —2Con = nant(n+1)an41. 

Simple expressions analogous to (15) and (18) do not seem to exist for m > 2. 
Relation (19), the inverse of (15), is easily obtained either directly from the 
definitions or from (15): 


(19) a, = : (—)" "7 (r+1)co 741. 
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In the design of spectrometers built with iron-cored coils the set of co- 
efficients a, is probably the most useful description of the field. For air-cored 
coils, on the other hand, convenient expressions are available for the vector- 
potential itself. We have found that the independent set of coefficients which 
is most easily computed from coil dimensions is not the Con, but is the one 
defined by Cmo, Cm1; the rest of the matrix and the a, may be obtained by using 
(12), (13), and (15) as recurrence relations in n. 


3. SECOND-ORDER SOLUTION 
To solve equations (9) and (10) we use the method suggested by Svartholm 
(1946). It is assumed the H, 7, h, t, and ¢ are all small compared to unity. 
If we retain terms in (9) and (10) as far as the first order only, then the approxi- 
mate equations (20) are the result: 


(20a) n +w;? = 2e, 

(206) rt’ +wer = 0. 

We have made the definitions w;? = 4¢92 and we? = 4¢29. The solutions of (20) 
corresponding to orbits originating in a point source on the optic circle at 
y = Oare: 

(21a) 7 = wrt sin(wiy) +2ew1-*[1 —cos(wiy)] ’ 


(210) tT = wo"'T sin(woy) . 


Double, i.e. two-directional, focusing occurs in the plane ¥ = Yo, if 
(22) wo = wo = m. 


The coefficients ¢o2 and ¢2o, upon which w; and w2 depend, are not independent, 
being related through 
(23) Coot+coo = 3, 


which relation is readily obtained from (13) in the special case m = 0. The 
unique solution of equations (22) and (23) is cor = C20 = 3, Wo = w. From (15) 
it is easily shown that 

(24) 2Co2 = 1+a. 


Thus double-focusing occurs, as Svartholm and Siegbahn (1946) discovered, 
for a; = —}4; the focal plane is situated at ¢o = Yor\/2 = 1/2. It may be seen 
from equation (21a) that the radial separation in the focal plane of the images 
of electron groups with fractional momentum difference « is 4ero. 

The results which we have obtained so far are well known. Before pro- 
ceeding to the second order, we point out that the first-order aberrations which 
occur when a; is not quite equal to —4 are easily found from equations (21) ; 
we give them, in expressions (25), for later use: 


(25a) n(x) = —xr(ai+3)H, 


(255) t(x) = w(ai+34)T. 


In the remainder of this work it is to be assumed that a; is exactly equal to —}. 
Also, since we are not interested in higher-order corrections to the dispersion, 
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we shall keep « = 0 henceforth. On the other hand, since we must deal with 
sources of finite size, we shall no longer restrict our discussion to orbits origi- 
nating on the optic circle. The appropriate first-order solution is given by the 


equations 
(26a) n = Hsiny +h cosy, 
(26d) t7=Tsiny +t cosy. 


As defined in the preceding section, h and ¢ are the initial values of » and 1, 
while H and T are the initial values of n’ and r’. If h = t = 0, then H and T 
are also the extreme values of 7 and 7 respectively ; these extreme values occur 
for Y = 2/2. Though h and ¢ are not actually zero, they are sufficiently small 
that H and T may be considered to be good approximations to the first-order 
maximum departures from the optic circle. In our later work we shall consider 
only those orbits which satisfy the conditions 
-H,<H< HM, -I,i<T<Th, 


and refer to the permitted range of H and T as the aperture. The required 
selection of orbits could be made, to the first approximation, by a rectangular 
baffle with an opening of dimensions 2Horo X27 oro placed in the plane y = 3z. 
The transmission -7 defined as the percentage of all electrons emitted by 
the source accepted by the spectrometer, is simply related to the size of the 
aperture. We define ¢, as the angle between the projection of the initial 
velocity vector on the plane of symmetry and the tangent to the optic circle 
at the point y = 0. Similarly, ¢, is the angle between this same tangent and 
the projection of the initial velocity vector on the plane which contains that 
tangent and a line parallel to the z-axis. It is not difficult to show that 


(27) /2 tan ¢, = Hand v2 tan ¢, = T. 


A sufficiently good approximation to the transmission is 
T = 2($r)o.2(b2)0-(44) X 100%, 


(¢r)o and (¢,)o being the maximum values of | ¢, | and | ¢,| respectively. 


Using (27) we find 
(28) TF = HoT (24) X 100%. 


We are now prepared to obtain the second-order solution. Following 
Svartholm, we write equations (9) and (10) in the form 
(29a) a’ +n = R(n, 7, 9’, 7’, H, T, h, 0), 


(290) +r = S(n, T, n’; '. H, Re h, t). 


The right-hand sides of equations (29) can be obtained as power series in all 
the parameters indicated ; there are no terms of degree lower than the second. 
We call the equations obtained by discarding third- and higher-order terms 
the second-order orbital equations. These equations are then approximated by 
the pair of equations formed by substituting the first-order forms of n and 7 (26) 
in the second-order terms. The simple, linear, inhomogeneous equations thus 
obtained are easily solved, subject to the initial conditions stated in the 
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preceding section. For the moment we give only the results for Y = 2, which 
we term the aberration; the solutions for general values of y are included in 
the lists of results in the next Section. 

(30a) n(x) = —h+3(1—8a2)H?+($a2—1)T? 


+32(1—2a2)h? + ($a2.—1)2, 
(300) T(r) = —¢t+2($a2.—1)HT+§azth. 


The relations (30) are precisely equivalent to the formulae obtained by Shull 
and Dennison (1947) by solution of the time-dependent equations of motion. 
Svartholm gave the terms in H and T only. 

Because the dispersion is in the radial direction, choice of the best value of 
a2 must be dictated primarily by the need to minimize the radial aberration. 
The second- and higher-order terms in h may be ignored, because the presence 
of the first-order terms in (30a) demands that h be exceedingly small. The 
limits of the first two terms remaining in (x) depend on the size of aperture, 
that of the last on the height of the source. To maximize the efficiency of the 
spectrometer we must minimize the aberration caused by each of the three 
terms. By proper choice of a2, namely az = §, %, or ?, we can make any one 
of these three terms zero, but not all three simultaneously. The least significant 
of the three terms is that proportional to #2. For a spectrometer of reasonable 
size the dimensions of sources available will usually be such that ¢ <1. Ifa 
source of greater height is to be used, it can be curved in such a way that the 
second-order aberration in é? is cancelled by the first-order aberration in 4; the 
central line of the source should fall along the curve 


(31) —h+(4a2.—1)2 = 0. 


The cancellation will hold, of course, only on the central line of the source, 
leaving an aberration proportional to the width of the source. 

The relative importance of the two remaining terms of the radial aberration 
depends upon the shape of the aperture. Ideally one ought to find the shape of 
aperture that, for fixed resolution, maximizes the transmission. The aperture 
shape so determined, and the corresponding value of the transmission, would 
be functions of a2 and of the desired resolution; the best value of a2 would be 
that which gives the largest transmission. Because of the neglect of the higher- 
order terms, the procedure just outlined is probably not worth the effort 
required. Instead we shall arbitrarily restrict the problem to that of finding 
the best aperture of rectangular shape. For simplicity we shall fix the basal 
line-width, ¢, rather than the resolution. Now we may write 
(32) ¢= pH?+qTe 


where p = }| 8a2—1 | and g = | $a2—1]. It is not difficult to show that the 
maximum transmission, under the constraint of constant ¢, is attained when 
(33) He = o/2p and T,? = o/2q. 


The axial and radial dimensions are given equal weight in making up the 
permitted line-width o; note that this does not occur for equal radial and axial 
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dimensions unless p = g, which equality implies a. = 3. With the optimum 
rectangular aperture we find 
(34) FT = o(4x)—(pq)-* X 100%. 


The quantity pg is roughly equal to 64/9 a,? for large values of | a2| ; it 
vanishes at a2 = 4 and at d2 = }; between the roots it reaches a local maximum 
at a2 = 3}. Thus the best choice of az is either $ or 3; values of a2 very much 
smaller than $ or very much larger than } give much poorer transmission ; in 
the range } < az < } the worst possible choice is az = }. This last statement, 
at first sight, seems to contradict. the findings of Rosenblum (1947), who 
showed that, with a given circular aperture, the smallest basal line-width is 
obtained for ag = }. The arbitrary restriction to circular apertures is respon- 
sible for the paradox. We can see no good reason for making such a restriction, 
except possibly when an instrument with large higher-order terms is used at 
high transmission. It is interesting to note that Rosenblum’s calculation per- 
formed for a square aperture gives a constant line-width for az in the range 4 
to 3; outside this range the line-width increases rapidly. 

The analysis of the preceding paragraph has shown that the best trans- 
mission for given resolution, with a variable rectangular aperture, occurs when 
a2 is either $ or 3. Consider the former case first. A large transmission is ob- 
tained because, for p = 0, Ho is permitted to become very large; the limit on 
Hy will be set by the third- and higher-order terms, which we have omitted 
from (32). These terms will limit Hy to values sufficiently small that ex- 
pression (28) for the transmissicn remains applicable. For this choice of az the 
value of g is 3; the value of J», fixed by ¢, is much smaller than that of Ho. 
Terms in (32) containing high powers of Ty or products of T,? with powers of ' 
Ho may be ignored. The dominant terms in o are the powers of Ho. Since 
Maxwell’s equations allow us one free parameter for each order used in the 
description of the magnetic field, we can make the coefficients of the powers of 
Hp in the expansion of ¢ all zero by proper choice of our free parameters. If, 
in an actual spectrometer, a few of the a, beyond az were fitted to their opti- 
mum values, a much larger value of Hy could be used. Such an instrument we 
shall call a wide-aperture spectrometer. The required values of the coefficients 
Gn, up to m = 4 or 5, have been published by several authors: Verster (1950), 
Pavinski (1954), Saulit (1954), and Huster et al. (1955). 

When we look at the other choice of az, viz. 3, we see that the roles of Ho 
and 7» are reversed. Here the residual second-order aberration is 4H,?; it is 
now Hp that must be kept small. The dominant terms in the expansion of T 
are now those proportional to T,?Ho, Tot, To*Ho, To®,.... Again there is one 
large term for each of the a,. We can choose the a, so that the expansion co- 
efficients can be made zero; a large value of JT) may then be employed. Such 
an instrument we shall term a high-aperture spectrometer. Calculation of the 
required optimum values of the a,, up to = 6, was the main purpose of the 
work which we are now reporting; these coefficients are supplied in the fol- 
lowing section. Their calculation is much more difficult than that of the 
corresponding coefficients for the wide-aperture type, because orbits not 
confined to the plane of symmetry must now be considered. 


LEE-WHITING AND TAYLOR: SPECTROMETER FOCUSING 9 


Though it is certain that either the wide-aperture or the high-aperture 
instrument, i.e. either ag = $ or a2 = 3, can be superior to one built with any 
other choice of az, it does not seem possible to ascribe an overwhelming superi- 
ority to one over the other. In the matter of transmission at given resolution 
the two spectrometers should be almost equally satisfactory, if the same 
number of coefficients a, is properly fitted for each. The odd-order dominant 
aberrations of the wide-aperture and of the high-aperture types of spectrom- 
eter are proportional to H,?"*' and to HoT," respectively; the latter has 
a slight advantage because it contains the smaller factor Hy. Nevertheless, 
superiority in transmission cannot be claimed for either type without reference 
to the method of providing the required magnetic field. For a given method 
of generating the magnetic field the highest transmission will be obtained with 
the spectrometer for which the largest number of a, can be fitted to their 
optimum values. 

There are, however, three minor advantages in favor of the high-aperture 
type: 

(1) By good luck the choice az = } also gives second-order axial focusing, i.e. 
the coefficient of HT in (306) vanishes; for the appropriate value of a; the 
third-order terms are small. For the choice a; = $, on the other hand, the 
second-order axial aberration is —$HT. The better axial focusing for the 
former case, though probably not important at apertures corresponding to 
0.01% resolution, might be useful at apertures corresponding to 1% resolution. 
For example, for the high-aperture case with H = 0.25, T = 0.4, we find that 
7(x) = 0.0088; for the wide-aperture case with H = 0.4, T = 0.25, we find 
that r(x) = 0.13. 


(2) For az = } the coefficient of £2 in n(x) is —}; for az = } this coefficient 
is —%. Hence linear sources 30% higher may be used for the former case with 
the same resolution. 


(3) The required magnetic field for the high-aperture type has a shape very 
near to that of the function 1/+/r. This means that the spectrometer will have 
good focusing properties for ro’s chosen from the whole fitted range of r, and 
therefore that it will focus electron groups of considerably different momenta 
simultaneously. This property is useful in that it enables one, with photo- 
graphic detection, to examine an appreciable part of the spectrum in one 
operation. 


4. HIGHER-ORDER CALCULATIONS 


Having decided to use a; = —} and az = 3, we next proceed to work out 
the optimum values of a, for m > 2. The first step is to solve equations (29) to 
the desired order of approximation by a generalization of the procedure by 
which we found the second-order solution. Because the effect of the finite size 
of the source can be estimated with sufficient accuracy in the second-order 
approximation, we now put 4 = ¢t = 0. The approximate solution of order 
(n+1) is obtained by substituting that of order m in the right-hand sides of 
the equation (29). The simple, linear, inhomogeneous differential equations 
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thus obtained may be solved in a straightforward manner; the amount of 
labor required mounts rapidly with the order. In this fashion we generate 
approximate expressions for 7 and r which are polynomials in H and T, the 
coefficients being functions of ¥. Though it would be difficult to prove that the 
approximate method is mathematically legitimate, it is true, as Svartholm 
points out, that the form of the solution obtained is that which is of interest 
practically ; the excellent agreement of the analytical formulae with the results 
of the numerical integrations, to be described later, shows that the procedure 
for obtaining the approximate solutions is indeed a very good one. 
We write 


(35a) 7 = a Fnn(y)H"T", 
(356) t= D0 Gnn(y)H"T". 


Because of the symmetry of the magnetic field, only even powers of T occur 
in 7 and only odd powers in +. All terms in the approximate solutions (35) as 
far as the fourth order have been obtained ; because of their possible usefulness 
in baffle calculations, we give these terms in full. 


Fy = sin y Fo. = —isin’y Fo = 2. sin'Y+3 (cos y—1) 
. 1 21 : 
Fy. = (i C22— 21) in ile C22— ry cos y—sin yp) 
F3o = ba a 1+(3eoe- 1g cos y—sin y) 
+2. sin ¥(cos y—1) 
Fu = ae +22 cots a 1) 
0 = \i5 C41 = css = cw 340 cos y— 
bar : 5on +59 tut 2 n'y 
+2 we Bo ; Feo 13h) in'y 
a cot zag y sin 2y 
Fas = (8¢62—-curtcrs)— 2) (cos ¥=1) 
2=\5 C2s— Coat C22) — F565 cos y— 
+(8(cas—cou) +4 Co2— — 81) in v 
2 
i ncnetem C22— — ein ¥ 


1 7 13 1 
+{den- 1) siny cos y+( 28 1) ¥ sin 
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+(1 on 3 ysin 2y 
Fy = (w Coo a) (cos y— »+(8 Coo— = 2 cou 355; aig y 
(Bvt on~ an (on~ shi cn 


(ie a ” +(3 ait sy “apy 
Go = sin vy Gu = $ sin *y 


Go3 = (cu! )sinty+(3 cot Ev cos y—sin wy) 
Ga = (3 cnt 1)in'y 
aed) —sin ¥)+3 sin ¥(cos ¥—1) 
“F 9 O82 G4 (y cos y—sin y g sin ¥(cos y 


128 64 32 11 
Gi; = (2 Cott. Cu “set (cos y—1) 


15 

+ nt Bea 10 en) 
+35 ots bad ca 35s" 
mi. cot’ ona sin 2y 


Ga = (2 carten—cu) +715) (000 y—1) 


16 77 bits. 
By eaten Co2— 0 cua) sin'y 


133 7 
+ = cout t39 Coa to €23— -81.\sinty 


+(1 otc )ain" y cos y— (: ontds)y sin Ww 


+(3 cute 5 aye sin 2y 


In addition, we give the values at y = = of those fifth- and sixth-order terms 
which affect the focusing significantly. 


121 1237 -3 





§ 
Fy4(1) _ Aafeo 48 Cu 384 ° 2 cutie 


145 a) 
49,152 


+o © cts 2 ast + Fates 
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Fos(n) = 2276 ¢,, 289 ,,, 4 3888192, 64. 
eS Te hae” ce ee 
4 12 32 8192 512 
+35 Cos — 35 Ct 35 Ciot 175 C40 +775 C22C41 
256 2 1024 4096 e 
175 “2? 175° Coxlaot— 175. Caveat e55 
Pe 2 
~4 (3 cote) 
a 357 105 15 
Gos(m) = dea 64°22 + jos Cot Cx —— 4 “tT6 394 a) 


The coefficients ¢mn upon which the functions Finn(y) and Gmn() listed depend 
can, as is shown in Section 2, all be expressed in terms of the coefficients an. 
Though putting a, equal to its optimum value before proceeding to work out 
terms of order (x+1) would have lightened the algebraic labor enormously, 
we did not do so because we wished to calculate orbits in fields other than the 
optimum, principally that of Stoker et al. (1954) and that which varies like 
1/v/r. 

The possession of formulae for Fin(y) and Fon(W) enables us to determine 
the optimum values of a, for the high-aperture spectrometer. From the 
equations 


Fin(r) = Oforn = 2 and 4, 
Fon(r) = Oforn = 4and6, 


we have calculated the required values of a3, a4, as, and as. These values are 
listed in the second column of Table I. In the third column of this table we 





TABLE I 
VALUES OF MAGNETIC FIELD EXPANSION COEFFICIENTS 
n Optimum 1l/vV/r Adn 
0 1 1 —_——— 
1 1 
1 “% —* —0.0013 
3 3 
8 5 0.0004 
43 5 
3 -—a = —0.2986 . ~ie™ —0.3125 0.0063 
553 35 
4 3304 = 0.2400... 128 = 0.2734... —0.0058 
69,709 63 
5 _ 345,600 = —0.2017 eee > 256 = —0.2461 ne —0.050 
60,299 35x? ” 231 
6 345,600 + 147,456 = 0.1768... 1024 = 0.2256... 0.075 
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give the expansion coefficients for a field which varies like 1/+/r. The similarity 
of the two sets of numbers is striking, and is probably responsible for the good 
focusing properties of instruments built with a field varying like 1/+/r. The 
significance of the numbers in the fourth column will be explained in the next 
section. 

By substituting the values of the cm», derived from the optimum values of 
the a, in the expansion for Fyn(x) and Ga,(r) we obtain formulae for the 
aberrations, viz. 


_ 28, 8 pars, 8 
(36a) a(x) = —“3-+7-H T'+57 Ht, 
—_ _% 7s_14 pp 
(365) r(x) = at 135 ff T+0.0196357". 


The relations (36a) and (365) are complete as far as the fourth order, and 
contain the important terms, for the high-aperture spectrometer, as far as 
the sixth and fifth orders respectively. The coefficients of H? in the former and 
HT, H°T, and HT? in the latter vanish accidentally. 

Because the values of a3 and a4 which we obtained did not agree with those 
published by Stoker et a/. (1954), namely —5/24 and 37/192, we thought it 
advisable to check our solutions by some other means. The unexpected occur- 
rence of x? in ads, which could not be removed by persistent recalculation, also 
excited suspicion that something might be wrong. Accordingly, we decided to 
compare the aberrations calculated analytically with aberrations obtained by 
a direct numerical integration of equations (9) and (10). We are much indebted 
to Mr. G. A. Chisnall* for carrying out these numerical integrations on 
FERUT, the Ferranti Computer at the Computation Centre of the Uni- 
versity of Toronto. The method used in the integration was that devised by 
Adams and Bashforth. The results obtained by numerical integration with 
the optimum field appear in the columns under the symbol MN in Table II. 


TABLE II 


COMPARISON OF NUMERICAL AND ANALYTICAL CALCULATIONS OF 
THE ABERRATION IN THE OPTIMUM FIELD 


—n(x) X108 —r(r) X10 
r H — 
N A N A 
0 +0.1 6.637 6.637 0 0 
-0.1 6.637 6.637 0 0 
0.1 +0.1 6.622 6.619 0.140 0.141 
0 0.000 0 0.131 0.131 
—0.1 6.618 6.619 0.119 0.121 
0.2 +0.025 0.412 0.413 1.040 1.041 
0 0.000 0 1.041 1.041 
—0.025 0.411 0.413 1.040 1.041 
0.3 +0.025 0.408 0.407 3.477 3.486 
0 —0.003 0 3.481 3.486 
—0.025 0.401 0.407 3.481 3.486 





*Now at the Department of Astronomy, University of Manchester, Manchester, England. 
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Aberrations calculated from the formulae (36) appear under the symbol A. 
The agreement is excellent, the small discrepancy which appears for larger 
values of T being of a size easily explained by the higher-order terms neglected 
in (36). Similar comparisons have been made with the same sort of agreement 
for the field of Stoker et al., and for that which varies like 1/+/r. We believe 
that we are justified, therefore, in claiming that our analytical solutions and 
our set of optimum values for the coefficients a, are correct and that those 


of Stoker et al. are incorrect. 


5. THE IMAGE OF A POINT SOURCE 


The expansions of the aberrations (36) in powers of the aperture dimensions 
are more likely to give a good description of the image when the aperture used 
is small, as it is at high resolution. We shall examine the image at an aperture 
corresponding to a resolution of 0.01%. First some suitable value of Hy) must 
be assumed. Put the residual second-order aberration, }H,?, equal to the 
separation in n(7) of two electron groups differing by 0.01% in momentum. 
In the absence of other aberrations and of finite source width, this procedure 
corresponds to a complete resolution of the two lines. From our knowledge of 
the dispersion of the instrument we can say that the separation in n(7) of the 
two groups is 4X10-*; hence Hp = 0.025. In Fig. 1 (a), (6), and (c) we show 
the images of a point source on the optic circle made by our optimum field, by 
the field of Stoker et al., and by the field of the form 1/+/r. Lines labelled 
H = +0.025 or H = 0 are the loci of the traces on the focal plane of orbits 
with the indicated values of H. Similarly the lines labelled with a value of T 
are the loci of the traces of orbits having that value of T. The quality of the 
fit to the optimum field which we have been able to obtain with a system of 
pairs of air-cored coils suggests that values of T up to J) = 0.3 are not un- 
reasonable. Inspection shows that the image made by the Stoker field is much 
inferior to that made by our optimum field. The image made by the 1/+/r 
field, though appreciably wider than that of the optimum field, is still very 
good. 

The transmission attainable at a resolution of 0.01% depends on the maxi- 
mum allowable value of 7». Fig. 1(a) shows that no appreciable broadening 
has occurred, for the optimum field, for values of T as large as 0.3. Though 
FERUT calculations indicate that the broadening of the image is negligible 
for T as large as 0.5, in any real instrument the effect of the higher-order, un- 
cancelled aberrations would probably be very serious at such large values of T. 
We shall content ourselves with pointing out that for JT) = 0.3 the trans- 
mission is about 0.12%. 

At a resolution of 1% the criterion for finding Hy which we used in the last 
paragraph gives Hy = 0.25. If values of To as large as 0.4 can be used, then a 
transmission of 1.5% should be achieved. 

The radial defocusing occurring when a, is not quite equal to its optimum 
value can be calculated from the formulae which we have given. For the first 
and second orders use (25a) and (30a) respectively; for higher orders use the 
formulae for Fim(x) and Fom(x). To give some idea of the closeness of fit 
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H=-—0.025 
H = 0.025 


- 





-4 -2 ° -4 =2 °o 
7m x 10° 7 x 10” 7 x 10° 


Fic. 1. Image of a point source on the optic circle; (a) refers to the optimum field, (6) to 
that of Stoker et al., and (c) to that of the form 1/+/r. 


required we give the numbers in the last column of Table I. The quantity Aa, 
is defined as the increment in a, producing an increment in n(7) of 10-4 at 
H = 0.025, T = 0.3; all earlier a, are assumed to have their optimum values, 
and the effect of the increment in a, is calculated in the mth order only. Under 
these conditions the relation between (m) and a, is approximately linear. 
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HALF-LIFE DETERMINATION OF SOME RADIONUCLIDES! 


W. F. Merritt, P. J. CAMPION, AND R. C. HAWKINGS? 


ABSTRACT 


The half-lives of the following nuclides have been determined by absolute 
counting techniques using a 4 beta proportional counter. The period over 
which the observations were made is indicated in brackets. The errors quoted 


are the statistical counting errors. 
Na® —2.58+0.03 years (2 years), 
Ru —1.02+0.01 years (5 years), 
Cs!4 —2,19+0.02 years (5 years), 
Ce! —285 +2 days (4 years), 
Pm"’—2.64+0.02 years (4 years), 
TI? —3.56+-0.05 years (4 years). 


A review of previous determinations is included. 


1. INTRODUCTION 


Since the introduction of the proportional 8 counter with 4r-geometry to 
this laboratory in 1950, accurate determinations of the disintegration rate of 
solutions of various radioisotopes have been made. Several of these nuclides 
have half-lives of the order of 1-5 years and a redetermination of the dis- 
integration rate of the solutions provides a direct method of measuring their 


half-lives. 
2. METHOD 


The counting techniques used have been fully described elsewhere (Hawkings 
et al. 1951). In most cases, the original sources were available for counting. 
Otherwise, new sources were made from the solutions which had been stored 
in pyrex bottles with greased glass stoppers. These were deposited on the same 
type of source mount used in the previous determination, usually a laminated 
zapon-formvar film of ~100 yug./cm.? superficial density, with a conducting 
layer of gold vacuum distilled on to both sides. The same counters were used 
in both series of observations. 

The following checks on the counting system were made at frequent intervals 
during the determinations: 

(i) Reference sources of Sr®° and of Ra D + E were counted at regular 
intervals to ensure that the system remained stable. 


(2) The resolving time of the system was measured several times during 
the course of the observations. This was done by observing the apparent 
counting rate of a high activity (35,000 c./s.) Sr®® source. The true counting 
rate of this source was previously determined in a counting system whose dead 
time had been measured directly by observing the decay of 24.99 minute I'8. 


41Manuscript received August 17, 1956. 
Contribution from the Nuclear a. I ae. Atomic Energy of Canada Limited, 
1 


Chalk River, Ontario. Issued as A.E 
*Present address: Development Chemistry. Branch, Atomic Energy of Canada Limited, 


Chalk River, Ontario. 
16 
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(3) Where necessary, the sources were counted in a low background labora- 
tory remote from large sources of radioactivity. This was done for sources with 
a lovy counting rate where the background correction was appreciable. 


To demonstrate the validity of the counting techniques, the disintegration 
rate of a sample of Cl* which had been used in a measurement of its half-life 
in 1953 was redetermined (Bartholomew eé/ al. 1955). This agreed with previous 
results to within 0.6%, well within the quoted standard error. 

In addition, sources of Ra D + E (Pb?!*) have been observed over a period 
of 53 years. These give a value for the half-life of 22.4+0.4 years, in agree- 
ment with the values quoted by Curie ef al. (1931) of 22 years, but in disagree- 
ment with Tobailem’s (1955) value of 19.40 years. 


3. CALCULATIONS 


From the equation 
(1) N = Noe? 


it follows that 
(2) Ty = 0.69315 T/In(No/N) 


where N, and N are the disintegration rates at the beginning and end of the 
period of observation. 7} is the half-life and \ the disintegration constant. In 
considering the errors involved in the determinations, it can be shown that by 
differentiating equation (2), the error in the half-life is related to the error in 
the counting rates by 
(3) AT, _ _-—T  ANo/N 

T, 0.693157 No/N ~ 


Fig. 1 shows the fractional error in Tj plotted against the period of observation 
measured in half-lives for a 1% error in No/N. This curve shows that the 
period of observation should be at least equal to one half-life and preferably 
more if reasonable accuracy is to be obtained. 


1.4 


12 


° 


ERROR (%) 
@ 


e  £.,%. 5. 6. .6..8:.8. + 
TIME (HALF-LIVES ) 


Fic. 1. A plot of equation (3) showing the fractional error in half-life as a function of the 
period of observation assuming an error of 1% in No/N. 
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For those nuclides where the original sources were available, the errors 
quoted are the standard deviations from the mean. In the case of those samples 
for which new sources had to be prepared from the original solutions, the errors 
quoted are calculated from the above equation using an estimate for the maxi- 
mum error in No/N. This estimate includes statistical errors and the spread in 
disintegration rate values as observed from several sources. 

4. MEASUREMENTS 
Na” 

Laslett (1949) observed a sample of Na” over a 3-year period with a 
Lauritsen electroscope and reported a half-life of 2.60 years. By irradiating 
fluorspar with polonium a-particles, Saha (1939-41) prepared a sample of Na” 
which he observed with a point counter. From his measurements he estimated 
a half-life of 2.8 years. 

We received a sample of Na” from the National Bureau of Standards, 
Washington, in 1954. Four sources prepared from this solution were counted 
immediately and recounted after a period of 2 years. A half-life of 2.59+0.02 
years was obtained. 

At about the same time a shipment of Na” was received from the Atomic 
Energy Research Establishment, Harwell. Because this sample had a high 
solids content, it was purified by elution from a 1 meter long column of Col- 
loidal Dowex 50 ion exchange resin, following the method of Beukenkamp and 
Rieman (1950), and its disintegration rate determined. A redetermination of 
this rate yielded a half-life value of 2.55+0.05 years, in good agreement with 
the N.B.S. sample. 

Because the original sources were available for the N.B.S. sample and be- 
cause the A.E.R.E. sample was of high solids content we quote a weighted 
average of 2.58+0.03 years. 


Ru'% 
Glendenin and Steinberg (1951) reported a half-life of 1.0 year for Ru'®, 
In 1951 we counted three sources prepared from a solution of Ru! which 
had been carefully purified by a distillation technique. A redetermination of 
the disintegration rate of these sources after 5 years showed a decay corre- 
sponding to a half-life of 1.02+0.01 years. 


C5134 

Kalbfell and Cooley (1940) prepared a sample of Cs'** by bombarding cesium 
with cyclotron-produced neutrons and deuterons. They obtained a value for 
the half-life of 20-1 months. Giencdenin (1951) irradiated cesium with thermal 
neutrons in a nuclear reactor, foiiowed its decay for 13 months using G-M 
counters, and reported a half-life value of 2.3+-0.3 years. 

We irradiated ‘‘Specpure’’* cesium chloride in the Chalk River NRX 
Reactor and removed the ch!vrine with an anion exchange resin. Five sources 
were recounted after a 5 year interval and a half-life of 2.19+0.02 years was 
obtained. 


*Trademark of Johnson and Matthey Company Limited, London. 
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Ce!4 

In 1944 F. Joliot reported a value of 290 days. Schuman and Camilli (1951) 
followed a sample of aged fission product Ce™ for two half-lives in a flow type 8 
proportional counter and reported a half-life of 285+3 days. Burgus et al. 
(1951) followed a sample for 660 days using G—M counters and found a half- 
life of 275 days. 

Sources were prepared from a solution of aged fission product Ce and 
counted in 1952. Four years later the solution was restandardized by the 
preparation and counting of fresh sources. A decay corresponding to a half-life 
of 285+2 days was obtained. 


Pm" 

Inghram et al. (1950) estimated the half-life of Pm" from the yield curve of 
the fission products as 2.26 years. Schuman and Camilli (1951) followed a 
source for 1 year and reported a half-life of 2.6+0.2 years. Melaika et al. 
(1955) measured the half-life of Pm from the growth of Sm with respect to 
Sm" using a mass spectrometer. Observations over a period of 5.6 years gave 
a value of 2.52+0.08 years. 

In 1952 we received a solution of Pm'’ which had been purified by ion ex- 
change techniques. Five sources prepared from this solution were counted 
immediately, and recounted after a period of 4 years. A half-life of 2.640.02 
years was obtained. 


T/2% 


The wide spread in reported values for the half-life of Tl? is shown in 
Table I. 


TABLE I 
Period of Source of 
Date Workers Apparatus observation material Half-life 
1941 Fajans and Bubble type G-M 300 days Cyclotron irradiation 3.5+0.5 
Voight counters, Lauritsen of thallium years 
electroscope 
1953 Locket and TQQ 8B electroscope 8 months 2.71+0.05 
Thomas years 
1953 Reynoldsand G-—M and propor- 15months Pile irradiation of 2.7 years 
Brooksbank __ tional counters thallium 
1953 Harbottle Flow type 8 pro- 2.9 years Oak Ridge thallium 4.0+0.1 
portional counter pile produced years 
1954 Horrocks and _ Lauritsen electro- 3 years Oak Ridge thallium 4.1+0.1 
Voight scope pile produced years 
1955 Cheng etal. Ionizationchamber 10 years Cyclotron irradiation 2.50+0.03 
of thallium years 
1955 Tobailem and Differential ionization 45 days U.S.A.E.C. TE 4.26+0.06 
Robert chamber years 





We have measured the decay of two solutions of Ti?, one for 44 years, the 
other for 3 years. The first solution was prepared from thallium metal which 
had been irradiated in the NRX Reactor in 1949. Subsequent spectroscopic 
examination of the metal has shown no significant impurity. This sample was 
dissolved in strong nitric acid and counted in 1951. The other solution was 
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produced by neutron irradiation of ‘‘Specpure”’ thallium carbonate in the NRX 
Reactor in 1953 followed by dissolution in strong nitric acid. The gamma-ray 
spectrum for this sample was examined with a 5 in. X4 in. Nal(T1) crystal for 
possible impurities. In the range 32 kev. to 1.7 Mev. only the characteristic 
X-ray of mercury was observed. 

Recently these solutions have been restandardized by the preparation and 
counting of fresh sources. A decay of the older and younger solutions corre- 
sponding to half-lives of 3.54+0.05 years and 3.59+0.05 years respectively 
was obtained. We quote a mean value of 3.56+0.05 years. 
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SPALLATION YIELDS FROM HIGH ENERGY PROTON 
BOMBARDMENT OF HEAVY ELEMENTS! 


J. D. Jackson 


ABSTRACT 


The Monte Carlo calculations of McManus and Sharp (unpublished) for the 
prompt nuclear processes occurring upon bombardment of heavy elements by 
400 Mev. protons are combined with a description of the subsequent neutron 
evaporation to determine spallation cross sections for comparison with experi- 
ment. The model ceed is a schematic one which suppresses the detailed 
characteristics of individual nuclei, but gives the over-all behavior to be ex- 
pected. Many-particle and collective effects such as alpha particle emission and 
fission are ignored. The computed cross sections are presented in a variety of 
different graphical forms which illustrate quantitatively the qualitative picture 
of high energy reactions first given by Serber (1947). The calculations are in 
general agreement with existing data when fission is not an important effect, but 
the agreement does not imply a very stringent test of the various features of 
the model. 


1. INTRODUCTION 


The qualitative features of nuclear reactions at high energies (over 100 
Mev.) were first described by Serber (1947) nearly ten years ago. Since that 
time a large amount of experimental work has been done on high energy 
spallation and fission yields for various types and energies of bombarding 
particle, and for targets of different mass numbers. Essentially all the work 
has been radiochemical in nature. It will suffice to refer to a few of the more 
recent papers from which earlier references can be obtained (Bennett 1954; 
Nervik and Seaborg 1955; Lindner and Osborne 1956). The experimental 
data confirm in broad outline the behavior predicted by Serber. At any given 
bombarding energy, a large number of different reactions proceed with 
comparable yield, while individual yields vary relatively slowly with energy. 
For all but the heaviest elements, the most likely reactions are those which 
give rise to product nuclei lying less than 10 or so mass units below the target 
mass, although nuclei 20 or more mass units below the target are formed in 
significant amounts. In the heavy element region, fission competes with 
spallation, giving a broad mass yield distribution centered somewhat below 
half of the target mass. The fission contribution increases rapidly as a function 
of atomic number—for tantalum bombarded by 340-Mev. protons fission 
occurs less than 1% of the time (Nervik and Seaborg 1955), while for uranium 
it occurs in more than 70% of the interactions (Jungerman and Steiner 1956). 

Serber described the impact of a high energy nucleon or group of nucleons 
with a target nucleus as a two-stage process. The first process is a series of 
quasi-free, two-body collisions between the incident particle and the nucleons 
within the nucleus in which, on the average, relatively sma!l amounts (~30 
Mev.) of energy are transferred. Commonly the incident particle either passes 
through the nucleus without any interaction, or makes a few collisions, 
knocking some of the nucleons out of the nucleus and perhaps emerging 


1Manuscript received August 30, 1956. 
Contribution from the Radiation Laboratory, McGill University, Montreal, Quebec. 
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itself. In this way a few prompt particles, carrying away varying amounts 
of kinetic energy, are emitted, leaving the residual nucleus with a wide range 
of different excitations. The second stage of the reaction is the dissipation of 
this excitation energy by the emission of further particles (evaporation 
particles) in the manner well known at low energies. These two steps give 
rise to the general behavior of the cross sections outlined above. 

A number of theoretical calculations have been made on different aspects 
of the problem of high energy reactions. Horning and Baumhoff (1949) 
discussed stars in photographic emulsions induced by 190 Mev. deuterons, 
using an analytic approach. Goldberger (1948) was the first to use the so- 
called Monte Carlo method to describe the prompt processes occurring upon 
collision of a high energy nucleon with a nucleus. More elaborate Monte 
Carlo calculations have been made by Bernardini, Booth, and Lindenbaum 
(1952) for high energy nucleons incident on medium weight elements, and 
in particular by McManus and Sharp (McManus, Sharp, and Gellmann 1954; 
McManus and Sharp, unpublished). The evaporation process has been 
treated by LeCouteur (1950, 1952), and Fujimoto and Yamaguchi (1950), 
among others. Recently the present writer has discussed (p, xm) reactions in 
heavy elements, giving a simplified description of neutron evaporation 
(Jackson 1956). 

The present paper presents a theoretical description of spallation yields 
from heavy elements bombarded by 400 Mev. protons. The Chalk River 
Monte Carlo calculations (McManus and Sharp, unpublished) are used for 
the description of the prompt processes, while the evaporation stage is de- 
scribed in a simplified way which suppresses the detailed properties of indi- 
vidual nuclei. Fission as a mode of de-excitation is ignored, although some 
qualitative remarks are made about it (see Section 5). 

The aim of the work has been to combine the theoretical descriptions of 
the prompt and evaporation processes to yield an over-all, schematic picture 
of high energy spallation reactions in heavy elements for explicit comparison 
with experiment. The statistical accuracy of the Monte Carlo computations 
is not high (see Section 2). Consequently all of the calculated cross sections 
are uncertain to some extent, aside from the inadequacies of the model em- 
ployed, and the less likely processes are given only as to order of magnitude. 
Nevertheless, it is hoped that the present crude.calculations will indicate the 
behavior expected from a model which ignores many-particle and collective 
effects, such as alpha particle emission and fission. 

A summary of the results of the Monte Carlo calculations is presented in 
Section 2. The model for the evaporation process is outlined in Section 3, 
while the calculated cross sections are given in Section 4. Comparison with 
experiment appears in Section 5, and a summary and conclusions in Section 
6. Mathematical details of the evaporation model are described in an Appendix. 


2. PROMPT PROCESSES 


The Monte Carlo calculations for the prompt cascade used in the present 
work are those of McManus and Sharp (unpublished). These authors em- 
ployed a weighted two-dimensional model using the observed nucleon- 
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nucleon scattering cross sections to describe the internal collisions and a 
Fermi gas model for the nucléons in the nucleus. They distinguished between 
neutrons and protons, and included Coulomb and centrifugal barriers for the 
emerging particles by means of cutoffs. The validity of the cutoff procedure 
was checked by a detailed computation of barrier penetration factors for a 
sample of the computer runs. Meson production was neglected. The com- 
putations were performed on the Ferranti computer at the University of 
Toronto Computation Centre for both neutrons and protons incident, and 
for several energies and target elements. The specific results used here are 
these for 400-Mev. protons incident on uranium, but because of the relative 
insensitivity of the prompt processes to small changes in nuclear charge and 
mass, they should be valid for any heavy element. 

The computations yield the number, energies, and directions of nucleons 
emerging from the nucleus for each particle incident. From this information 
the charge, mass, and excitation of the residual nucleus can be inferred. For 
the purposes of evaluating reaction yields only the number of neutrons and 
protons emerging and the state of excitation of the residual nucleus are 
relevant. The results of computer runs for 200 protons incident at 400 Mev. 
on uranium are shown in Figs. 1 and 2. 

Fig. 1 presents a three-dimensional histogram of the frequency distribution 
of prompt emission of protons and neutrons. The most probable event is the 
emission of one neutron and one proton, corresponding generally to a single 
n—p collision within the nucleus. The distributions are not sharply peaked, 
however. The mean number of prompt protons (neutrons) emitted per inter- 
action at 400 Mev. is 1.3 (1.6), while the r.m.s. deviation from the mean is 
0.9 (1.3). The interaction cross section is 0.84 times the geometrical cross 


section of the nucleus. 





NEUTRONS ——~ 


Fic. 1. Calculated frequency distribution of numbers of prompt protons and neutrons 
— from a heavy nucleus bombarded by 400-Mev. protons (McManus and Sharp, un- 
published). 
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Fic. 2. Calculated distributions of residual excitation for several types of prompt events. 
The pair of numbers (i, 7) correspond to ejection of 7 prompt neutrons and 7 prompt protons 
(McManus and Sharp, unpublished). 


In Fig. 2 are shown the residual excitation distributions for various types 
of prompt events. Only the more probable are presented; even then the 
statistical accuracy is such that it is necessary to average over 40-Mev. 
intervals in excitation energy. The events are characterized by two numbers 
(i, j) corresponding to the prompt emission of i neutrons and j protons. 
A noteworthy, if expected, feature of these distributions is the peaking towards 
low energies for those events in which a single proton, two protons, or a 
neutron and a proton emerge. By and large these events correspond to one 
dominant nucleon-nucleon collision within the nucleus, with the collision 
partner perhaps escaping and perhaps not. Relatively small amounts of 
energy are imparted to the residual nucleus. For interactions in which several 
neutrons and/or protons are ejected, the excitation distributions are broader 
and extend to higher energies, corresponding to the development of a size- 
able internal cascade. 


3. NEUTRON EVAPORATION 


Dissipation of the excitation energy remaining after the prompt processes 
occurs by the successive evaporation of one or more nucleons or groups of 
nucleons in the well-known manner first elaborated by Weisskopf (1937), 
(Blatt and Weisskopf 1952). In general the detailed calculation of the relative 
probabilities of evaporation is very difficult, and only approximate solutions 
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on the average behavior and the fluctuations have been obtained (LeCouteur 
1950, 1952; Fujimoto and Yamaguchi 1950). In heavy nuclei, however, a 
simplifying assumption can be made. Because of the high Coulomb barriers 
the emission probabilities for charged particles are very small compared 
to those for neutrons, and one can to a good approximation neglect all but 
neutron evaporation. 

As will be seen in the next section the averaging of the evaporation probabi- 
lities over the broad »xcitation distributions of Fig. 2 means that only the 
mean behavior of the evaporation process is important. The details of a 
simplified evaporation model (Jackson 1956), which assumes a constant 
nuclear temperature, and its averaged properties are presented in an Appendix, 
along with the average behavior of neutron evaporation assuming a nuclear 
temperature varying as the square root of the maximum available energy. 

For comparison with experiment it is necessary to know neutron binding 
energies and empirical nuclear temperatures. Very few neutron binding 
energies are known experimentally, especially on the neutron deficient side 
of the stable valley, and one must estimate them from semiempirical mass 
formulas. In the region of bismuth as a target the mass formula implies that 
the average neutron binding energy (apart from the odd-even alternation 
of about +0.6 Mev.) increases by about 0.13 Mev. per neutron emitted, 
and decreases by about 0.2 Mev. per proton emitted for fixed neutron number. 
For simplicity, in the calculation of the cross sections it was assumed that 
the average binding energy of successive neutrons was constant and equal 
to 7.5 Mev., regardless of the mass and charge of the residual nucleus. In 
view of the statistical uncertainties in the Monte Carlo calculations, as well 
as lack of knowledge of the nuclear temperature, this simplification will not 
add appreciably to the uncertainty of the results. 

Experimental data on nuclear temperatures at high energies are few. 
Measurements on the energy spectrum of evaporation neutrons from 190- 
Mev. proton bombardment of various targets (Gross 1956) imply an average 
temperature of 1.90.1 Mev. and an average excitation energy of 85+15 
Mev. for gold and uranium targets. In a similar experiment with 40-Mev. 
protons Heckman (1956) inferred temperatures of about 1.4 Mev. for heavy 
elements. On the other hand the shapes and positions of (p, xm) excitation 
functions for lead and bismuth from 10 to 90 Mev. imply a nuclear tem- 
perature of around 1.8 Mev. (Jackson 1956). A reasonable theoretical expec- 
tation is that the excitation energy E should be roughly proportional to the 
square of the nuclear temperature 7. Lang and LeCouteur (1954) have made 
an analysis of the available data at low and high energies to determine an 
empirical relation between temperature and excitation for various mass 
numbers. In rough approximation their empirical equation can be replaced 
by the expression JT? = ¢ E, where the level density parameter e has the value 
¢ =~ (8.5/A) Mev. The data of Gross and Heckman for A ~ 200 imply 
e =~ 0.04 Mev. in good agreement with the empirical expression, while the 
(p, xn) excitation functions of lead and bismuth at 30 to 40 Mev. would 
imply e ~ 0.10. 
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Fic. 3. Mean excitation energy for the evaporation of a given number of neutrons as a 
function of the number of evaporated neutrons. The heavy straight line corresponds to a 
constant neutron binding energy of 7.5 Mev. and a constant nuclear temperature of 1.9 Mev. 
The other curves correspond to constant neutron binding energy B, but a temperature varying 
as the square root of the excitation energy, T? = eE. The various assumed binding energies 
and coal density parameters are indicated on the individual curves in Mev. 


The essential results of the neutron evaporation calculation outlined in 
the Appendix are presented in Fig. 3, where the mean excitation energy 
for the evaporation of a given number of neutrons is plotted as a function 
of the number of evaporated neutrons. The various curves illustrate the 
assumption of constant nuclear temperature, of the square root of energy 
variation of temperature, and of different values of neutron binding energy 
and level density parameter. It will be seen that for excitations up to 200 
Mev. the various curves shown (corresponding to a variation of 1 Mev. in 
binding energy and a factor of two in level density parameter) do not differ 
by more than +10 Mev. in excitation energy for fixed neutron number or 
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+1.5 in neutron number for fixed excitation. In view of the broadness of the 
distributions of residual excitation shown in Fig. 2 these various possibilities 
will not lead to qualitatively significant differences in cross sections, although 
the detailed behavior will depend on the particular dependence assumed. 
In the calculations of Section 4 two examples are considered: Assumption 
A—constant binding energy (7.5 Mev.) and constant temperature (1.9 Mev.); 
Assumption B—constant binding energy (7.5 Mev.), and level density pa- 
rameter, ¢ = 0.075. The differences in calculated cross sections for these 
two cases are illustrated in Figs. 5(@) and 5(d). 


4, SPALLATION CROSS SECTIONS 


The calculation of the individual reaction cross sections follows closely the 
pattern of the author’s paper on (p, xm) reactions (Jackson 1956). Let 
Q(i, y, E)dE be the probability per inelastic interaction for the prompt 
emission of 7 neutrons and y protons, leaving behind a residual nucleus with 
excitation energy between E and E+dE. These probabilities for 400 Mev. 
protons incident are shown in Fig. 2 for the more likely events. The total 
probability for the ejection of 7 prompt neutrons and y prompt protons, 
illustrated in Fig. 1, is given by 
(1) q(i, y) = JOG, y, E)dE. 


The probability P(j, E) that a nucleus of initial excitation E will evaporate 
exactly 7 neutrons is described in the Appendix. With the assumption of no 
charged particle evaporation, the probability per inelastic collision that 7 
prompt neutrons, j evaporation neutrons, and y protons will be emitted is 
evidently 


(2) p(i,j,y) = [QG, y, E)P(j, E)\dE 
and the cross section for the emission of x neutrons and y protons is 
(3) o(xn, yp) = 0 2) P(i,4,9) 

{+ jaz 


where a; is the total inelastic interaction cross section (approximately 0.84 
times geometrical for 400-Mev. protons on heavy elements, from the Monte 
Carlo calculations). With the neglect of the ejection of aggregates of nucleons, 
a(xn, yp) is the cross section for formation of the nuclide of mass number 
(x-+-y—1) less than, and atomic number (y—1) less than, the target nuclide. 
Owing to the comparatively sharp peaking of the neutron evaporation 
probabilities around the mean energies E, (see Appendix) relative to the 
excitation distributions of Fig. 2, it is legitimate to approximate (i, j, y), 
Eq. (2), by making a Taylor series expansion of Q(i, y, E) around E,. Then 
Eq. (2) reads 
(4) pli,j,y) ~fPG, EdE [Q(i, y, By) +4(4E,? Oi, y, B+ «.-1- 
The integral f P(j, E)dE can be interpreted approximately as the mean total 


energy (binding plus kinetic) carried off by a neutron evaporated from a 
nucleus of excitation E,. This is exact for the constant temperature model 
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described in the Appendix. If this energy is denoted by (B+27) ,, the probabil- 
ity p(t, j, y) can be written approximately as: 


the second derivative of Q(i, y, E) with respect to F being ignored. 

Using Eqs. (3) and (5) and the information contained in Figs. 2 and 3 
the cross sections were evaiuated for the two assumptions about neutron 
evaporation described at the end of the previous section. The results are 
presented in Figs. 4 and 5. The ordinate scale in both figures is millibarns, 
based on a total inelastic cross section at 400 Mev. of 1.8 barns. The statistical 
fluctuations in the Monte Carlo calculations give rise to corresponding 
fluctuations in the spallation cross sections. The scattered computed points 
and smooth curves drawn by eye through the points are both shown in Fig. 4 
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Fic. 4. Calculated spallation cross sections o(xn, yp) for 400-Mev. protons on a heavy ele- 
ment (A~200) as functions of x for various values of y, making assumption A (see text). 
Note the difference in ordinate for parts (a) and (b). The points represent the numerical values 
obtained with the Monte Carlo calculations. The curves are the smoothed cross sections 
Rey ed eye through the points. The ordinate scale is based on a total inelastic cross section 
of 1.8 barns. 
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to indicate the statistical accuracy involved. In Fig. 5 only the smoothed 
curves are shown. The cross sections in Fig. 4 are for assumption A (constant 
temperature), while those in Fig. 5 are for assumption B. The manifestations 
of the differences in the evaporation process implied by assumptions A and B 
are also presented in Fig. 5, where the cross sections for assumption A (dotted 
curves) are shown for comparison. It will be seen that the over-all appearance 


eo OTT eee 


| 
4 


j ‘. : (a) 





(mb) 





O (xn, yp) 


(mb) 





O (xn,yp) 


Fic. 5. Calculated spallation cross sections o(xn, yp) for 400-Mev. protons on a heavy 
element (A ~200) as functions of x for various values of y, making assumption B. The dotted 
curves show the corresponding results for assumption A. The ordinate scale is based on a 


total inelastic cross section of 1.8 barns. 


of the two sets of cross sections is very similar, the only essential difference 
being that the cross sections for B are shifted slightly towards larger neutron 
number in the region of appreciable cross section. This behavior is what would 
be expected from Figs. 2 and 3. Fig. 2 shows that most of the excitations lie 
below 160 Mev., and Fig. 3 shows that at around 100 Mev. assumption B 
implies emission of roughly one more neutron for a given excitation than 


assumption A. 
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It should be noted from Figs. 4 and 5 that for reactions involving the 
é€mission of a small number of nucleons the only cross sections which are 
appreciable are those in which one or two protons are ejected. The reason 
for this, as <liscussed at the end of Section 2, is the relatively small amounts 
of excitation remaining after the prompt (p, p’), (p, pm), and (p, 2p) pro- 
cesses, which correspond mainly to single internal two-body collisions. For 
y = 1 and 2, the dotted curves in Fig. 4(a@) and 4(b) show the decomposition 
of the cross section into contributions from these single collisions and from 
the larger internal cascades. Once the ‘‘single collision’ part is subtracted, 
the residues of the y = 1 and y = 2 cross sections behave very similarly to 
the other cross sections, which peak broadly around 12 nucleons as the most 
probable number to be emitted. 

The calculated cross sections of Figs. 4 and 5 are shown in somewhat dif- 
ferent form in Fig. 6. The sum of individual cross sections for fixed mass 
number gives the so-called “chain yield” cross section. This is plotted in 
Fig. 6(a) as a function of mass number below the target. For about 10 or 












































Fic. 6. (a) Calculated cross section for the production of nuclei of a given mass number 
by 400-Mev. protons on a heavy element (A -~200) as a function of mass number, A, below 
the target mass number, Ao. The total inelastic cross section is assumed to be 1.8 barns. 

(6) Contour diagram on the N-Z plane of calculated spallation cross sections for 400-Mev. 
protons on a heavy element (A ~200). The shaded square represents the target nucleus plus 
a proton. The dashed line is the line of stability, while the solid contours are based on Fig. 5. 
The dotted 10-millibarn contour shows the possible effect of proton evaporation. 
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12 mass numbers below the target the yield is roughly constant, but it then 
falls off rapidly to about one fifth of its maximum by 20 mass units below 
the target. 

The cross sections are displayed as a contour diagram on the N-Z plane 
in Fig. 6(b). The ‘‘compound nucleus” (target plus a proton) is shown shaded. 
The line of stability is given by the dashed curve. The solid contours are 
based on the cross sections of Fig. 5. To illustrate the possible effects of 
neglected proton evaporation an estimated 10-millibarn contour, based on 
the calculations of Fujimoto and Yamaguchi (1950), is dotted in on Fig. 
6(5). For less than 10 mass units below the target there is no significant effect, 
but at 20 mass units below there is an appreciable swinging of the contour 
downward towards the line of stability. Unfortunately, it is unlikely that 
an experimental check can ever be made on the individual yields so far from 
the stable line. 


' 5. COMPARISON WITH EXPERIMENT 


The calculated cross sections shown in Figs. 5 and 6 can be compared with 
available experimental data of Bennett (1954) on the spallation of bismuth 
by 400-Mev. protons, and of Nervik and Seaborg (1955) on the spallation of 
tantalum by 340-Mev. protons. The effect of fission can be observed by com- 
parison with the work of Lindner and Osborne (1956) on non-fission, inelastic 
interactions of thorium and uranium with 340-Mev. protons. 
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Fic. 7. Comparison of observed and calculated spallation yields for 400-Mev. protons 
on bismuth, shown as a function of atomic number for various fixed mass numbers. The 
dashed curves are the calculated results of Fig. 5. The solid points are the data of Bennett 
(1954). Only the shapes of the experimental yield curves are significant, not their absolute 


values. 
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Bennett’s cross sections, while quoted in millibarns, are not really absolute 
yields, but are only relative among the different isobars of a given mass 
number. Yields of different mass number chains are not interrelated, no 
correction having been made for counter efficiency. A comparison between 
Bennett’s data and the calculated cross sections is shown in Fig. 7, where 
the yields of various sets of isobars are plotted as a function of atomic num- 
ber. For all mass numbers, Bennett’s quoted cross sections have been multi- 
plied by a factor of four before being plotted as the solid points. This is not 
an unreasonable factor on the average to compensate for counter efficiency, 
although it might be grossly incorrect in specific instances. For A = 205, 
the open circles are Bennett’s values multiplied by 20. The dashed lines 
represent the calculated cross sections from Fig. 5. It will be seen that the 
shapes of the calculated yield curves are generally in good agreement with 
the experimental data, and that the over-all correction factor of four gives 
fair absolute agreement. The chief exceptions are for A = 205, where the 
shape agreement is excellent, but there is a factor of 20 between Bennett's 
values and the calculated ones; and for A = 200, where the experimental 
yield peaks at Z = 82 as compared to Z = 83 for the calculated yield. The 
first exception cannot be regarded as significant until something more is 
known of the decay scheme and radiations involved. The second disagreement 
may well reflect the effect of the neglected proton evaporation in the calcula- 
tions, although the apparent shift in peak with atomic number is somewhat 
larger than expected. 

Nervik and Seaborg’s work on the bombardment of tantalum with 340 
Mev. protons is devoted to both the spallation and fission regions, although 
the fission cross section is only about 4 millibarns. Aside from the somewhat 
lower bombarding energy, their spallation data are not very suitable for 
detailed comparison with the present calculations because, on the one hand, 
there are a number of stable nuclides in the region of presumably high yield, 
and on the other, an appreciable fraction of their data are cumulative chain 
yields only. Their Fig. 6 (Nervik and Seaborg 1955) shows the estimated 
cross section for formation of a given mass number as a function of mass 
number. For mass numbers less than about 10 mass units below the target 
the yield is relatively flat at around 150 millibarns; by 20 mass units below 
it has fallen to about 40 millibarns and it is assumed to fall rapidly beyond 
that. This behavior is roughly consistent with the results shown in Fig. 
6(a), although the experimental estimate is 50% larger than the calculated 
one. 

The effect of fission as a mode of decay can be observed clearly by com- 
paring the present calculations, which ignore fission, with non-fission spalla- 
tion yields from nuclei with appreciable high energy fission cross sections. 
Such data are recently available from the work of Lindner and Osborne 
(1956). These authors observed spallation yields from 340-Mev. protons on 
thorium and uranium. For thorium, fission occurs in roughly 45% of the 
interactions, while for uranium it occurs more than 70% of the time. The 
observed spallation cross sections for thorium are shown in Fig. 8, along 
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Fic. 8. Comparison of observed (Lindner and Osborne 1956) and calculated spallation 
yields for 340-Mev. peotene on thorium. The open circles are for the Pa (Z = 91) isotopes, 
the solid circles for Th (Z = 90), the crosses for Ac (Z = 89), and the squares for Ra (Z = 
88). The solid curves are the calculated cross sections which neglect fission as a mode of 


decay. 


with the calculated cross sections of Fig. 5. The somewhat smaller binding 
energy for neutrons in thorium as compared to bismuth and lead compen- 
sates crudely for the reduction in bombarding energy from 400 Mev. to 
340 Mev., so that the cross sections of Fig. 5 would be expected to give a 
reasonable description of these data if fission were unimportant. The signi- 
ficant feature of Fig. 8 is the fall-off of the yields of the Pa, Th, and Ac isotopes 
with decreasing mass number as compared to the calculated curves. The 
yield of the heaviest isotope of each element is in rough agreement with 
calculation, but the yield of the lighter isotopes falls off rapidly. For the 
radium isotopes, all the data are in qualitative accord with the calculated 
curve. 

As has been emphasized by Lindner and Osborne, the observed behavior 
finds an easy explanation in terms of competition between fission and neutron 
evaporation. The probability of fission, apart from fine effects, depends 
sensitively on the parameter Z?/A. This means that for a given mass, the 
nuclei of higher atomic number will be more prone to fission, and for a given 
element, the nuelei of smaller mass will be more fissile. On the other hand, 
for a given excitation, the probability of neutron evaporation will remain 
essentially unchanged for a group of neighboring nuclei. After the prompt 
processes a number of different residual nuclei are left with a wide range 
of different excitations. As neutrons are evaporated from these nuclei, the 
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values of Z?/A increase and the probability of fission as a mode of decay 
rises accordingly. This competition from fission will tend to deplete the 
spallation yields of the lighter isotopes of a given element, and will be the 
more important the larger the atomic number of the element. Both of these 
effects are found experimentally in thorium (see Fig. 8). While the observed 
yields of the heaviest isotopes of all elements are in rough accord with the 
calculated cross sections, the yields of the lighter isotopes tend to be less, 
the effect being progressively more pronounced as the atomic number goes 
from Z = 89 (Ac isotopes) up to Z = 91 (Pa isotopes). For the Ra isotopes 
(Z = 88), the atomic number is apparently sufficiently small that fission is 
not important, at least for the mass numbers shown in Fig. 8. In uranium, 
the observed fall-off of the spallation yields is much more pronounced, as 
expected from its larger fission cross section. 


6. CONCLUSIONS 


The present paper contains an attempt to combine the Monte Carlo calcula- 
tions of McManus and Sharp (unpublished) on the prompt nuclear processes 
occurring when 400-Mev. protons are incident on a heavy nucleus with a 
description of subsequent neutron evaporation to yield spallation cross 
sections for comparison with experiment. The model employed is one that 
suppresses detailed characteristics of individual nuclei, and gives an over- 
all schematic picture. The agreement with existing experimental data can be 
considered adequate where fission is not an important effect. It should be 
made clear, however, that this agreement does not afford a very sensitive 
test of the various features of the model. All that can be said is that the 
Monte Carlo calculations, based on Serber’s original ideas, seem correct in 
broad outline, at least, and that the dissipation of residual excitation pro- 
ceeds mainly by neutron evaporation. There may be some evidence for 
proton evaporation where more than 10 nucleons are emitted, but apparently 
no need to postulate the evaporation of aggregates of nucleons (e.g. alpha 
particles) as a significant mechanism of de-excitation. 
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APPENDIX 
The details of a simplified model of neutron evaporation, assuming a con- 
stant nuclear temperature, as well as the average behavior of a model in 
which the temperature varies as the square root of the available energy, 
are outlined in this appendix. 


Constant Temperature Model 


If it is assumed that the energy spectrum of evaporated neutrons is given 
by U exp(— U/T), where U is the kinetic energy of the neutron and T is the 





JACKSON: SPALLATION YIELDS 35 


constant nuclear temperature, then it can be shown that the probability 
that a nucleus of initial excitation energy E will emit exactly 7 neutrons 
is: 


(Al) P(j, E) = (Ay, 27-3) —I(Ayui, 27-1) 
where I (z, 2) is Pearson’s incomplete gamma function 
(A2) (2,0) = 4, [ste ae 

nN: Jo 


j 
and A, = (e- ¥ 3.) / ? 
t=1 


is the energy (in units of JT) above threshold for the emission of j neutrons, 
B, being the binding energy of the ith neutron. The expression (A1) is in 
error by terms of order exp(—B/T) where B is a typical neutron binding 
energy. For 7 = 1, 2 it agrees with the results of Weisskopf (Blatt and 
Weisskopf 1952). The behavior of P(j, E) is illustrated in Fig. 2 of the paper 
by Jackson (1956) for a typical choice of parameters. 

In the arguments of Section 4 it is useful to have certain integral properties 
of the probability function P(j, EZ). First of all the integral of P(j, E) over 
all energies is: 

(A3) SPGQ, E\dE = By +2T, 


which is just the average energy (binding plus kinetic) carried off by the 
(j+1)th neutron. For simplicity of discussion in what follows it will be 
assumed that the successive neutron binding energies are all constant and 
equal to B. Nothing essential to understanding is lost by this assumption, 
and the various formulas are considerably simplified. The mean excitation 
energyE, for the emission of 7 neutrons is: 
= JEP(j,EdE . (B’—2T”) 
A BE, = —— = }(8B+-27)+-—__| 
(A4) L™ ipa mar "PHD, 
As expected, this is linearly proportional to the mean energy carried off per 
neutron, apart from the small finite intercept at 7 = 0. The energy spread 
in the distribution P(j, E) can be measured by the r.m.s. deviation AE, from 
the mean energy E;: 

(AE,)* = f (E—B,)’P(j, E)dE/f P(j, E)dE. 
Performing the necessary integrations using Eqs. (Al) and (A2) gives the 
result: 


2_ 1 (BY+8B°T+12B°T°—127") , ...,.70 


Again, apart from the small constant term, this expression 1s just what one 
would naively write down from the fact that the mean squre deviation 
for the individual spectrum, U exp(— U/T), is equal to 27”. 

The table below illustrates numerically the mean energie: and r.m.s 
deviations for B = 7.5 Mev., T = 1.9 Mev., the values used in “ig. 3. It 
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j 2 4 6 ~ 10 12 14 16 18 20 
E; (Mev.) 25 47 70 92 150 ol 160 182 205 228 


AE; (Mev.) 4.7 6.0 ooh 8.0 8.8 9.6 10.4 ° 02.0" EG es 





will be noted that the r.m.s. deviations are very small compared to the mean 
energies; the probability distributions peak relatively sharply around their 
mean energies. 


Mean Behavior of T? = eE Model 

Elementary arguments treating the nucleons in a nucleus akin to a degen- 
erate Fermi gas lead to the expectation that the nuclear temperature, 7, 
should vary approximately as the square root of the excitation energy, E. 
Other models lead to similar power laws with slightly different exponents. 
The mean behavior of a model based on T? = e£, where ¢ is a ‘‘level density 
parameter’’ with dimensions of energy, can be obtained approximately by 
analogy with the constant temperature model just described. The level 
density parameter ¢ will be assumed to be essentially constant for the various 
nuclei involved in the neutron evaporation process. Actually, € varies slowly 
with mass number on the average, and has rapid fluctuations depending on 
the odd—even character of the nucleus. 

If a nucleus has a large initial excitation E, the maximum energy available 
for evaporation of the first neutron is E,; = E—By,, where B, is the first 
neutron binding energy. The temperature for the evaporation spectrum is 
therefore 7, = (e£,)'/2, and the mean total energy carried off by the first 
neutron is U; = B,+27). On the average then, the maximum energy avail- 
able for second neutron emission is E, = E;—U,—Bz2, the average tem- 
perature is 7, = (e£,)!”, and the mean total energy carried off by the second 
neutron is U,; = B,+27>. The emission of successive neutrons proceeds in 
this way until the maximum energy available at some stage vanishes. If the 
mean initial excitation energy for the emission of j neutrons is desired, one 
must solve the simple set of j simultaneous algebraic equations: 


r A-1 
(A6) E, = E-— > B,—2>. T,(E,) 

u=1 u=1 
where \ = I, 2,..., 7. For the last neutron emitted there is some ambiguity 


as to whether the energy available should be equal to zero or some finite 
amount. Taking the constant temperature model as a guide, one should put 


(A7) E, = 2T (E;) 
as the last of the 7 equations (A6) in order to give the intuitive result, 


3 
E,= » (B,+2T7), 


when constant temperature is assumed. Choices other than (A7) are possible, 
but the particular value assumed for the energy available for the last neutron 
is of little consequence as far as the mean excitation energy is concerned, 
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except for the first few neutrons. In practice, it is simplest to start at the 
“bottom” and find successively the mean energy for the emission of one 
neutron, two neutrons, three neutrons, and so on. For various choices of 
parameters, the results of a computation based on Eqns. (A6) and (A7) are 
presented in Fig. 3, along with the mean energy for the constant temperature 
model, given by Eq. (A4). 
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THE EFFECTS OF GRAIN BOUNDARIES ON THE PLASTIC 
DEFORMATION OF ZINC CRYSTALS! 


G. B. CRAIG AND B. CHALMERS? 


ABSTRACT 


The tensile plastic deformation of single-crystal and tricrystal specimens of 
zinc was investigated by analysis of the external change in shape of the specimens, 
and of the changes in X-ray diffraction patterns. It was demonstrated that the 
single crystals deformed by slip on the basal plane, but pyramidal as well as basal 
slip occurred in specimens containing grain boundaries. 


INTRODUCTION 


The purpose of the present investigation was to analyze the mode of de- 
formation of zinc single crystals and to determine how the deformation process 
was affected by the presence of grain boundaries. 

The production, deformation, and calculations of the mode of deformation 
of zinc single crystals presented little difficulty. The selection of a specimen 
which would demonstrate the effects of grain boundaries was more difficult. 
Ideally the study should be conducted on a crystal which is completely en- 
closed by neighboring crystals, but the problem of producing and measuring 
the change in shape of such a crystal has not been solved. Therefore, as a first 
approximation to an enclosed grain, a crystal with two neighbors was chosen. 

Zinc was used because only one slip plane (0001) is reported at room tempera- 
ture. Thus the advent of any other deformation mechanism should be readily 
apparent. 

In order further to simplify the problem it is desirable to eliminate any 
deformation due to twinning, and also to minimize the internal strains set up 
by the anisotropy of thermal expansion of the individual zinc crystals. These 
requirements are satisfied if the angle between the specimen axis and the basal 
plane is nearly 45 degrees because: 

(a) The shear stress resolved on the basal plane is high, 


(6) Twinning occurs most readily when the stress axis is perpendicular 
to the C-axis, 


(c) The thermal strains which arise during cooling down after solidification 

of the crystals from the melt are minimized. 

To this end, a “‘tricrystal”’ test piece consisting of three crystals running 
longitudinally and comprising a rectangular bar 6 in. long, ? in. wide, } in. 
thick was evolved. The two outside crystals had identical orientations with 
the C-axis parallel to the top surface and inclined at an angle of 45 degrees to 
the specimen axis. The center grain was also oriented with the C-axis parallel 
to the top surface and at 45 degrees to the specimen axis, but this latter rota- 
tion was in the opposite sense. 





‘Manuscript received September 21, 1956. 
Contribution from the Department of Metallurgy, University of Toronto, Toronto, Ontario. 
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EXPERIMENTAL TECHNIQUES 


Single-crystal and tricrystal specimens were produced by means of the 
technique described by Chalmers (1940). 
Zinc was supplied by the New Jersey Zinc Co. Ltd., Palmerston, N.J., U.S.A. 


It had the following analysis: 
Zinc — 99.999% 


Pb — 0.00013 
Cd — 0.00005 
Fe — 0.0006 
Cu — 0.00009 
As — 0.000009 


The surfaces of the specimens were polished electrolytically using a technique 
described by Vernon and Stroud (1938). 

The standard Laue back-reflection technique was used to determine crystal 
orientations before and after deformation. 

It was suggested by Orowan (1952) that more detailed examination of the 
Laue spots and also a measure of lattice disturbances over a wide area could 
be obtained by a modification of the Laue technique due to Berg. This con- 
sisted of using a slit (2 mm.X7.2 mm.) to define the line focus of the X-ray 
beam. The remarkable detail shown by these pictures is shown in Figs. 8 and 9. 

In order to measure the change in external shape during deformation a 
reference network of lines was photoengraved on the top and bottom surfaces 
of the specimen. The ‘‘Pitman Deep Etch Process’, a standard photoengraving 
technique, was found to be most satisfactory for this purpose. Fig. 7 shows a 
portion of a specimen with engraved grid. 

The mechanism of deformation is deduced from the external changes in 
shape of the specimens. A system of grid lines parallel and perpendicular to 
the specimen axis (tension axis) permits an analysis of this change in shape 
to be made. 

The measurements made, and the co-ordinate system used for identifying 
angles during the present investigation, are shown in Fig. 1. All measurements 







7 face4_ ic 
<* topsurface 7 


| ‘ 
<> right side 


L 


Fic. 1. Measurements used to calculate distortion in metal crystals. 


except thickness were made with a travelling stage and microscope and are 
accurate to +0.01 mm. The thickness was measured with a micrometer with 
spherical tip attachments to provide point contact. This measurement was 
made to the nearest 0.0001 in. but is subject to the same inaccuracy as the 
other measurements because of the etched nature of the surfaces. (Etching for 
10 seconds was necessary during the photoengraving process.) The lengths 
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“b” and “‘c” along the grid lines shown in Fig. 2 are calculated after taking 
into account the angles y, 8, and X, e.g., distances were measured parallel and 
perpendicular to the specimen axis and the actual length of ‘‘c’’ is governed by 
the angle y. 

The angles y, 8, and X were calculated from measurements of the top and 
bottom grid lines. The co-ordinate system used for the top surface of the 
specimens was maintained for the bottom surface by using a specially machined 
jig with the same origin for co-ordinates on the top and bottom surfaces. 

The measured co-ordinates of the top and bottom grid lines produce directly 
the extension (e = final length/initial length), the dimension ‘‘c’”’, and the 
angle y. These co-ordinates plus the thickness measurement enable the di- 
mension ‘‘b” and angles 6 and X to be calculated as shown in Fig. 2. 





_Fic. 2. Angular relationship between y, 8, X, and the dimensions ‘‘t”’ and “b". Co-ordinate 
distances are in millimeters, e.g. A—(0, 0) origin, D—(—0.1, —0.2) from origin, —microm- 
eter thickness. 


The Taylor analysis determines the directions in the crystal which are un- 
changed in length during the deformation. Now if the deformation is by means 
of slip on a certain crystallographic plane this fact will be evident immediately, 
for all the directions in this plane will be unchanged in length. In fact that 
which is being determined is the cone of intersection between the strain 
ellipsoid and the undeformed sphere, hence the term ‘‘unstretched cone” is 
used frequently in this investigation. The cone degenerates into two planes if 
the simple slip process occurs. The plane on which slip actually took place 
can only be determined by further deforming the crystal and noting that one 
of these planes is crystallographically unchanged while the other varies. Ii is 
concluded that the unchanged plane is the slip plane. Further consideration 
leads to the conclusion that the slip direction is 90 degrees from the inter- 
section of the two planes and is in the slip plane. 

The most convenient way of picturing the analysis is by means of the stereo- 
graphic projection and this method is used throughout the investigation. The 
front side of the specimen shown in Fig. 1 is the plane of projection and the 
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polar co-ordinates of the directions in the specimen which are unchanged in 
length during the deformation are calculated using the intersection of faces 1 
and 4 as the origin for 6, and face 1 as the plane where = 0. Thus @ is the 
angle that the projection of the undistorted line makes with face 1, and @ is 
the angle that the line makes with the intersection of faces 1 and 4. 

In order to carry out the type of analysis planned it was essential that no 
bending or twinning of the crystals occurred. These conditions were fulfilled 
in all cases, except for one specimen in which twinning occurred. 


EXPERIMENTAL RESULTS 
(a) Taylor Analysis 


Measurements of the changes in external shape of Number 7 single crystal 
are given in Table IA. The data in Table IA have been used to calculate the 


TABLE I 


NUMBER 7 SINGLE CRYSTAL 


A. DIMENSIONS OF MEAN PARALLELEPIPED 
BEFORE AND AFTER DEFORMING 


Before After 
b 5.80 mm. 5.23 mm. 
c 2.53 mm. 2.53 mm. 
l 2.53 mm. 2.79 mm. 
e = 1.1028 
r 87.8° 85.6° 
B 86 93.2 
y 90 91.9 


6b = length of cross grid on faces 1 and 3. 

c = length of cross grid on faces 2 and 4. 

1 = length of longitudinal grid on faces 2 and 4. 

e = final length/initial length. 

dX = angle between faces 1 and 4. 

B = angle between specimen axis and imaginary 
grid joining equivalent points on faces 2 
and 4. 

y = angle between specimen axis and cross grid 
on face 4. 


B. POLAR CO-ORDINATES OF UNSTRETCHED DIRECTIONS 





Unstrained material Strained material 

n %; A 6,’ D2 02 62’ 
—5 —67.2 57.1 112.2 —71.3 54.6 116.0 
—2 —42.1 38.2 116.4 —46.2 35 123.9 
-—1 —23.9 32.3 117 —26.6 28.9 126.4 

0 0 29.7 114.3 0 26.6 124.7 

1 23.2 31.5 107.8 24.9 29.2 117.4 

2 40.1 36 101.1 41.9 34.5 109.1 

5 63.6 49.7 90.7 63.8 49.3 95.5 





polar co-ordinates (, 0) for the directions which are unchanged in length 
during deformation. These co-ordinates are recorded in Table IB and plotted 
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on stereographic projections in Figs. 3 and 4. The stereographic projections 
show that in this case the unstretched cones have degenerated into two planes. 

A similar analysis was made on the center crystal of Number 8 tricrystal. 
The external changes in shape of this crystal are given in Table IIA. The polar 
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Fic. 3. Number 7 single crystal. Stereographic projection of positions in unstrained ma- 
terial of cone of directions which are unstretched after material has been elongated 10.3%. 
Theoretical slip direction 90 degrees from intersection of the two planes is marked. Pole of 
0001 plane and slip direction a; determined from Laue back-reflection patterns are also indi- 
cated. 
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Fic. 4. Number 7 single crystal. Stereographic projection of positions in strained material 
of cone of directions which are unstretched after material has been elongated 10.3%. Theo- 
retical slip direction is marked. Pole of 0001 plane and slip direction az determined from back- 
reflection Laue patterns are also plotted. 
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TABLE II 


NUMBER 8 TRICRYSTAL—CENTER CRYSTAL 


A. DIMENSIONS OF MEAN PARALLELEPIPED 
BEFORE AND AFTER DEFORMING 


Before After 
6 6.32 mm. 6.23 mm. . 
c 2.53 mm. 2.29 mm. 
l 2.52 mm. 2.81 mm. 
e= 1.115 
d 87.5° 89.5° 
B 100.9 100.5 
¥ 89.6 84 


B. POLAR CO-ORDINATES OF UNSTRETCHED DIRECTIONS 


Unstrained material Strained material 

n ?, 6; 6,’ P, 02 6,’ 
—2 —39.6 57.6 109.5 —36.5 55.1 113.6 
-1 —22.2 71.3 107.3 —20.3 69.3 109.5 

0 0 78.8 116.3 0 75.6 117.7 

1 21.5 72.3 131.5 20.1 66.6 133.8 

2 37.7 67.8 138.6 36.2 60.5 141.4 

4 56.2 65 143.8 55.4 56 147.3 


co-ordinates of the unstretched directions are recorded in Table IIB, and the 
unstretched cones are plotted in Figs. 5 and 6. It is at once evident that the 
unstretched directions no longer lie on two planes. 





Fic. 5. Number 8 tricrystal. Stereographic projection of positions in the unstrained 
material of cone of directions which are unstretched after material has beri: -longated 11%. 
Position of the basal and pyramidal planes as determined by back-reflection Lave patterns 
is also plotted. 
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Face | 
2:0 





Fic. 6. Number 8 tricrystal. Stereographic projection of positions in the strained _ma- 
terial of cone »/ directions which are unstretched after material has been elongated 11%. Posi- 
tion of basal a »d pyramidal planes as determined from Laue patterns is shown. 


(b) X-Ray Diffraction 

X-ray diffraction patterns representative of those obtained on deformed 
single crystals and the center crystal of deformed tricrystals are shown in 
Figs. 8 and 9. The Laue back-reflection patterns from the center crystal of 
Number 8 tricrystal show a most pronounced form of streaking after de- 
formation. The Berg patterns from the center crystal are also greatly trans- 
formed as shown in Fig. 9. 


DISCUSSION 


Single Crystals 

The Taylor analysis of the change in external shape of Number 7 single 
crystal, as plotted on the stereographic projections in Fig. 3 for the undeformed 
sample, and Fig. 4 for the deformed sample, reveals that in this case the un- 
stretched cone has degenerated into two planes. When the pole of the basal 
plane, as determined from back-reflection Laue patterns, is plotted on these 
projections it is at once evident that one of the undeformed planes is the 
(0001) plane. The theoretical slip direction, which is 90° from the intersection 
of the two undeformed planes, and the position of the close packed direction 
(a2), determined from Laue patterns, are also plotted. These directions are seen 
to be within 6 degrees of each other. The analysis of the external change in 
shape, therefore, shows that all the observed deformation can be accounted for 
by means of simple slip or glide on the basal plane in a close packed direction. 
This confirms the simple picture of the deformation process arrived at by 
earlier investigators. 

The Laue X-ray patterns also support this conclusion, i.e., that the de- 
formation proceeds by means.of simple glide on the basal plane, in two ways: 
firstly, the basal plane rotates towards the specimen axis in the close packed 
direction dictated by the maximum resolved shear stress law; secondly, the 
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i Fic. 7. Grid engraved on zinc tricrystal specimen. G. I. Taylor analysis (Taylor and Elam 
' 1922-1923). 
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Laue spots themselves are sharp after deformation, indicating little or no 
disturbance of the regular crystal symmetry. Fig. 8 shows that deformation 
has not altered the sharpness of the diffraction spots. The undistorted nature 
of the Berg pattern shown in Fig. 8 is also evidence that simple slip has oc- 


curred. 


Tricrystals 

The observations on the mode of deformation of the center crystal of 
Number 8 tricrystal will now be considered. If we plot the position of the basal 
plane, as determined by Laue patterns, on the stereographic projections, 
Figs. 5 and 6, which show the unstretched directions in this crystal, it is seen 
that there is some relation between the major axes of the cones and that of 
the basal plane. It is not at once evident what this relation is or why it exists. 
Since the basal plane intersects the cone at a direction of close packing, i.e., a1, 
we are led to consider the possibility of a second slip plane which intersects 
the basal plane in this line. A consideration of lowest energy dislocations 
(Chalmers and Martius 1952) shows that the next operative slip plane should 
be a type I pyramidal plane (0111). All these pyramidal planes intersect the 
basal plane along the closest packed directions (a;, a2, a3). That pyramidal 
plane which has the maximum resolved shear stress acting on it should operate. 

The (0111) plane has the maximum resolved shear stress and, therefore, 
should operate if pyramidal slip occurs. This pyramidal plane is plotted in 
Figs. 5 and 6 and we see that the cone, the basal plane, and this pyramidal 
plane all intersect near the a; direction. This strongly suggests that the ob- 
served deformation is due to slip on the basal plane plus slip on the pyramidal 
plane, since the analysis of the external change in shape shows that this 
direction, i.e., the line of intersection, is unchanged in length during the de- 


formation. 
X-Ray Evidence for Pyramidal Slip Plus Basal Slip 

The zone streaking which appeared on the Laue back-reflection patterns 
from the deformed crystals supports the view that pyramidal plus basal slip 
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Fic. 10. (a) Stereographic projection showing initial position of Laue spots (. dots). 
Calculated position of Laue spots due to basal slip (O circles). Calculated position of Laue 


spots due to pyramidal slip (A triangles). : J 4 
(b) Stereographic projection of the streaking observed in Laue back-reflection 


patterns from the center crystal of the deformed tricrystal. 
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occurred. The streaks were plotted on a standard stereographic projection, a 
quarter section of which is shown in Fig. 10(6). Now if we consider the rota- 
tions of the lattice planes which take place, first for basal slip and then for 
pyramidal slip (Fig. 10(a)), we see that the lines joining this orientation range 
correspond to our streaking X-ray pattern. In a single crystal the whole 
lattice rotates uniformly from initial to final position, whereas in this case 
portions of the lattice are rotated to a greater or lesser extent in the slip 
direction. This is, of course, a special case of deformation bands and .perhaps 
will lead to a better understanding of the nature of these bands. ° 





Fic: 11. Stereographic projection showing how the “unstretched cone” degenerates into 
two planes when the rotations of the lattice due to pyramidal slip are reversed. First rotation 
(P 37°) makes (0111) plane the basic circle. Second rotation (P 7.5°) reverses the rotation 
due to pyramidal slip. C-axis and basal plane (dashes) are plotted from X-ray data. 


On the assumption that pyramidal slip is operative, the component of the 
lattice rotation arising as a result of this slip system can be calculated to be 
7.5° (Fig. 10). If we reverse the deformation (theoretically) due to pyramidal 
slip we arrive at a position of the unstretched cone which is due to basal slip. 
This operation has been carried out, and the results are shown in Fig. 11. We 
see that the cone has reverted to two planes and our X-ray data show that one 
of these planes is the basal plane. The theoretical glide direction is also shown 
and we see that it differs by only 7 degrees from one of the directions of closest 
atomic packing (a2) in the slip plane. In view of the errors inherent in the 
experimental techniques this is regarded as satisfactory agreement. 

Only slip lines due to basal glide were visible on the surfaces of the deformed 
crystals. One explanation for this is that the etched surface masked the 
pyramidal glide bands. 

There is, however, direct X-ray evidence in both the Laue and Berg patterns 
that disturbances of the lattice in directions other than those accountable for 
by basal plane glide took place. The streaking in the Laue patterns, and the 
pronounced break-up of the Berg patterns, are so markedly different from those 
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obtained on Number 7 single crystal which deformed by pure slip that there is 
not the slightest doubt that some other mechanism is operative. 

The work of Rohm and Kochendorfer (1950) shows that homogeneous shear 
deformations up to 20-25% in long aluminum crystals produce no macro- 
scopic slip bands and no asterism in Laue back-reflection X-ray patterns. 

It is, therefore, considered that the absence of visible slip bands due to 
pyramidal slip does not contradict the argument for deformation by this 
mechanism. 

The analysis reported was carried out on a very simple crystal configuration 
and it is certain that in normal polycrystals many further complications arise. 
Two cases of interest in this respect were encountered in the course of the 
present investigation. 

First was a tricrystal similar in orientation to Number 8 tricrystal. It had, 
however, one more complicating factor as far as the deformation process is 
concerned. This was the wavy or irregular nature of the two grain boundaries. 
This so altered the stress distribution that twin crystals formed in the center 
crystal. These twin crystals made a Taylor type analysis of the external 
change in shape impossible and also spoiled the X-ray patterns as far as meas- 
uring the rotation of the specimen axis was concerned. 

The presence of large angle striation* boundaries (2-10 degrees) in a single 
crystal also altered the mode of deformation. A complete analysis of this case 
was not made but it appears likely that the operation of a second slip system 
could account for the observed changes in the deformation process. 


CONCLUSIONS 


1. It is concluded that the presence of grain boundaries in zinc crystals 
modifies the simple basal slip mechanism which is observed when striation-free 
single crystals are deformed. The modification in the case investigated was 
the operation of an additional pyramidal slip plane. This result could not be 
predicted from tests carried out on single crystals. 


2. The observed deformation of a striation-free single crystal (Number 7) 
can all be accounted for by means of slip on the basal plane in the close packed 
direction. 
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NUMERICAL ANALYSIS OF A UNITARY PARTICLE MODEL' 


F. A. KAEMPFFER AND EpitH M. KENNEDY 


ABSTRACT 


A classical unitary field model describing a family of electrically uncharged 
heavy particles is analyzed numerically. A reasonable mass spectrum: is ob- 
tained. Limitations of the model are stated. 


1. INTRODUCTION 


Ever since Mie (1912) published his famous program, unitary theories of 
the so-called elementary particles have remained in what might be called 
a state of animated suspension. New attempts in the spirit of Mie keep 
appearing in the literature, but although they raise high hopes by establishing 
the feasibility of a unitary theory, none of them has yet succeeded in pre- 
dicting the outcome of a single experiment. 

The aim of the investigation reported in this paper was not to enlarge on 
the general feasibility of unitary particle theory, but rather to study the 
numerical feasibility of a known model, which appears attractive from a 
general point of view, and, if possible, to develop a feeling for the order of 
magnitude of the parameters characterizing this model. 

In particular, the family of heavy particles, consisting of nucleons and 
hyperons, was taken as the reality to which the model under consideration 
should correspond. The reason for this particular choice is that all unitary 
models to be discussed below describe a ‘‘particle’’ in terms of self-maintained 
solutions of certain nonlinear differential equations, which allow for the 
existence of a ‘‘ground state’ together with ‘“‘excited states”, a situation 
resembling most closely the relation between the various members of the 
heavy particle family. 

It will be found that one can obtain reasonable values for the mass spec- 
trum. The description of the decays of the various particles and the predic- 
tion of their lifetimes, however, is beyond the reach of present unitary theories. 


2. EXISTING UNITARY THEORIES 
Existing unitary theories may be separated naturally into three groups: 
(1) Purely classical field theories, in which the cohesion of particles is 
provided by either 
(la) coupling of the field to itself, or 
(1b) coupling of two fields, each acting as the glue which holds the other 
together. 
Models resulting from this approach have been studied by Rosen (1939), 
Finkelstein (1949), and others. 


‘Manuscript received August 13, 1956. 
roe from the Department of Physics, The University of British Columbia, Van- 
couver 8, B.C. 
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(2) Quantum field theories, with coupling of the field to itself. This ap- 
proach has been followed by Heisenberg (1954) and others. 

(3) Field theories in which a classical ‘‘source’’ is coupled to a quantized 
“‘glue”. The feasibility of this intermediate type of model was first pointed 
out by Heber (1955). 


The interpretation of models resulting from the approach (2) meets with 
great difficulties, as has been shown by Kita (1956), and so far no model of 
type (3) is known which will correspond to a particle with positive mass. 
Therefore the work reported in this paper was confined to a model of type 
(1). 

3. SUMMARY OF FORMALISM 

The field equations governing a classical unitary theory may be derived, 

as usual, from a variational principle requiring that the action integral 


(3.1) I =fZfd% 


be an extremum, where “= 7 (y,, dy;/dx,,...) is the Lagrangian density. 
The conservation laws for momentum, energy, and electric charge follow 
from invariance properties of the Lagrangian (see e.g. Wentzel 1943) and 
take the form 

(3.2) OT ,,/dx, = 0, 0s,/dx, = 0, 


where energy momentum density tensor 7T,, and charge current density 


vector s, are defined by 
__ af oa 
O(Op:/Ox,) Ox, 
fa ee ‘) 
sit ditties Yi 5(ay ax) ) 
Now any solution of the field equations following from (3.1), which is free 
from singularities, is interpreted as a particle of mass 


(3.4) m= —* f Tud'x 
and charge 


(3.5) e= -* J sd°x. 


It is evident from the definition of the charge current density vector that a 
theory with real fields ¥; only describes an electrically uncharged field and 
therefore neutral articles, and that charged particles may be obtained 
only if at least one of the field variables is complex and time dependent. 


4. THE SPECIFIC MODEL 


For the numerical investigation a specific model was chosen, characterized 
by the Lagrangian density 


= = +Z%6,,, 


(3.3) 
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* 
(4.1) - - 3 


_ 96 96 


* 2 
OX, ax, EY % ’ g > 0. 


eyy* 
The corresponding field equations are 
(4.2) Oy—Ky = —gyv¢’, 


(4.3) O¢ = —gw*¢. 


By analogy with orthodox quantum field theory one may think of these as 
describing a “‘bare’’ nucleon field y of ‘‘bare’’ mass m = xh/c coupled to a 
massless meson field ¢. The main reason for this choice is its simplicity and 
the fact that particles of positive mass result from it. It will also be shown 
later that assignment of a ‘“‘bare’’ mass yp to the meson field ¢ does not alter 
appreciably the features of this model, provided p/m <1. 

To make the nonlinear simultaneous system of equations (4.2) and (4.3) 
accessible to numerical treatment only static spherically symmetric solutions 
of the type 


(4.4) v=Vv(r), &=9(7), real, 
are considered. By the substitutions 
(4.5) x = xr, “i= Vz ry, v= Vare 


the field equations reduce to 


(4.6) oli ie) . 
(4.7) -o(%) 


which should yield particle-like solutions if they are subjected to the boundary 
conditions (see Rosen 1939) 


u=0 
v=0 


u/x—-0 


ane eee ©, 


(4.8) ate =0, 


By inspection of the differential equations (4.6) and (4.7) it is seen that 
u" < Oifv/x > 1, u” > Oifv/x <1, 
v’ < 0forallx > 0. 


(4.9) 


It is therefore to be expected that the differential equations and boundary 
conditions will in general allow solutions with u having no nodes, one node, 
two nodes, etc. If each such solution is interpreted as a particle, one has a 
model which allows, in principle, the calculation of a mass spectrum of various 
members of an electrically neutral particle family described by these equa- 
tions, using definition (3.4) with 


(4.10) Tu = —{(’—g¢?)¥+(V¥)?+(V¢)*} 


in this case. 
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5. APPROXIMATION PROCEDURE 


Since an exact solution of the equations (4.6) and (4.7) seems out of reach, 
the following approximation procedure was adopted. Let, in crude approxi- 
mation, 

_ Jv, = ax forx < xo, 

(5.1) mer . = aXxo for x > Xo, 
where a and x» are constant parameters to be determined later. Then equation 
(4.6) becomes 

» _ Ju(1—a’) for x < Xo, 
(5.2) mo cc (axo/x)*} for x > xo. 
For a < 1 one has from (5.2) u”’ > 0 for all x, and for a = 1 one has u” = 0 
for x < x) and u” > 0 for x > xo. In either case u/x does not vanish as x >, 
in contradiction to the boundary condition. Therefore one requires 


(5.3) a>l. 
Equations (5.2) now yield, together with the boundary conditions, 


(5.4) aie 12 sin(\/a*—1 x) for x < xo, 


C exp(—~x) for x > 1. 


This solution is approximated again by writing 


bh 1" = B sin(./a?—1 x) for x < Xo, 


(5.5) Us = C exp(—x) for x > Xo, 


with B and C constants, and requiring that u and u’ be continuous at x = Xp. 
This latter requirement yields 


asl 
(5.6) aa C exp(—%o) : a tan” ( Va 1) 
sin(/a?—1 Xo) a—1 

where one must have 

4/2 <Va?—1 x < x for the nodeless solution, 
(5.7) 34/2 < Va—1 xo < 2m for the one-node solution, 

54/2 <Va?—1 x9 < 3m for the two-node solution, etc. 
From (5.6) one can eliminate B and a, using (5.7) for the respective solutions. 
The model is therefore characterized by the four parameters C, xo, g, x. For 


the numerical calculation it was found convenient to express C in terms of 
the normalization constant m, defined for dimensional reasons by 


a 
(58) Sve ==, 
so that is a pure number and is an energy density. The parameters fixing 
the model are thus 7, xo, g, x. Substituting then the solutions into (3.4), using 
(4.10), and performing an integration by parts, one obtains 


(5.9) H = mc = festy'ate at meme r:) 
g Y; 


where 
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2 
Ye me eee se 1 


Va?—-1 sin’(\/a’—1 Xo) 


x {\/a?—1 xo—cos(r/a?—1 x9) sin(y/a?—1 x)}, 


Y; = 2a*x9—4(axo)” exp (2x0) J exo(=4) du, 
2z0 


1 
sin?(./a?—1 Xo) 
X {\/a?=1 xo—cos(+/a?—1 x0) sin(x/a?—1 x0)}, 


with a to be determined numerically from (5.6) and (5.7) for the nodeless, 
one-node, two-node, etc. solution respectively. 


Y,= 1+ 


6. NUMERICAL RESULTS 


The dimensionless quantity Hg/nx was calculated numerically for a range 
of values of the parameter x) for each of the three cases corresponding to 
no node, one node, and two nodes in y respectively. Solutions with more than 
two nodes in y were not considered, although such solutions do, of course, 
exist. The results of this calculation are tabulated in Table I and plotted in 
Fig. 1. 


Hi 
se} Ahn 


2 


10 is 20 
Xo 


Fic. 1. The variation of Hg/nx with xo. At the beginning of each curve is shown the 
number of nodes in the corresponding y-field. 





KAEMPFFER AND KENNEDY: UNITARY PARTICLE MODEL 53 


TABLE I 


VALUES OF Hg/nx AND CORRESPONDING VALUES OF X9 FOR THE NODELESS, ONE-NODE, AND 
TWO-NODE SOLUTIONS 


Nodeless solution One-node solution Two-node solution 
Xo Hg/nx Xo Hg/nx Xo Hg/nx 
0.70 3.377 1.50 4.249 2.50 4.409 
1.20 2.336 2.75 2.546 4.00 2.708 
2.00 1.950 4.50 2.033 6.00 2.153 
2.40 1.9186 5.00 1.9942 8.00 2.001 
2.45 1.9148 5.50 1.9777 8.50 1.9911 
2.50 1.9185 5.75 1.9759 9.00 1.9884 
2.60 1.921 6.00 1.9778 9.50 1.9921 
3.00 1.950 6.50 1.9906 11.00 2.030 
4.50 2.211 7.00 2.013 18.00 2.479 
7.50 2.987 10.00 2.270 
10.00 3.713 


The graph of Hg/nx versus x) shows that there is one value of xo in each 
case for which Hg/nx, and hence the energy, H, is a minimum. The existence 
of these minima is the most remarkable feature of this model. The authors 
argue that if this model has any relation to reality, one should consider the 
particular solutions which yield the minimum energy in the three cases as 
the ones representing three possible neutral particles in their normal states 
and use the corresponding values of r = xo/x as a measure of the size of 
these particles. The mass ratios are independent of n, x, and g within the limits 
of the approximations used and are found to be 
(6.1) My 1M: m2 = 1: 1.032 : 1.038. 


The ratios of the radii are also independent of m, x, and g in the same sense 
and are found to be 
(6.2) 79271 27%. = 1: 2.34 : 3.67. 


The values of xo are those found by approximating Hg/n« by a parabola in 
the vicinity of the respective minima. 
By comparison the mass spectrum of the known neutral members of the 


heavy particle family is described by the ratios 

(6.3) My : My, : my = 1:1.18 : 1.26. 

Expressed in terms of n, x, g the rest mass my and the radius 7 of the smallest 
particle described by the model are 

(6.4) my = 2.13X%10-*! nx/g gram, 79 = 2.45«—' cm. 

By equating these expressions with the values found experimentally for the 


nucleon, namely my = 1.67X10-* gram and 7% = 0.8X10-'* cm., one obtains 
the further relations 


(6.5) x = 3.06X10" cm“ 
and 
(6.6) n/g = 2.56X10-'" c.g.s. 


It is interesting to note that the value (6.5) for « corresponds to assigning the 
‘“‘bare’’ nucleon field y a rest mass 








54 
(6.7) 


The theory does not yield m and g separately by comparison with experi- 
mental values. In view of the simplicity of the model and the crudeness of the 
approximations used it is felt that the theoretical values (6.1), (6.2), and (6.5) 
can be considered as reasonable. 

One may try to improve the theoretical mass spectrum (6.1) by assigning 
a “‘bare’’ mass not only to the y-field but to the ¢-field as well. This has been 
investigated and it was found that the general method developed in Section 5 
will still work provided the ‘‘bare’’ mass yu of the ¢-field is small compared 
with the “‘bare’’ mass m of the y-field. In particular, by choosing the ratio 
n/m to be 1/6, assuming thus that the ‘‘bare”’ mass ratio reflects approximately 
the observable mass ratio of the corresponding ‘‘dressed”’ particles, it was 
found that the resulting modification of the first half of the mass spectrum 
(6.1) is 
(6.8) my my’ = 1 : 1.034. 
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It appears thus that such modifications will not yield any significant new 
information and this matter was, therefore, not pursued any further. The 
assumption that ‘“‘bare’’ mass ratios reflect the observable mass ratios of 
“‘dressed’’ particles is highly dubious in any case. 

It is obviously beyond the capacity of the specific model treated in this 
paper to describe instability and predict the lifetime of certain particles. 
Quantization of electric charge also cannot be grasped within the framework 
of such a purely classical model. In fact, if ¥(r) is a solution of the field equa- 
tions, then Yexp(iwt) is also a solution for any value of w, and it is easily | 
verified that the minima of the curves Hg/nx versus x» occur for the same 
values of xp as in the case w = 0. By definition (3.5) solutions with w #0 
describe an electric charge 


(6.9) e = ewn/cag, a = P—w/ce 







carried by the particle, where ¢ is a scale factor. Clearly, e remains com- 
pletely undetermined in the classical theory. One may entertain the hope, 
however, particularly in view of Heber’s (1955) interesting results, that the 
minima of the energy, occurring for certain values of the spatial extension 
of the particles, are a feature of the theory which will survive incorporation 
of quantum postulates in this type of model. 
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MULTILAYER PROBLEMS IN THE SPHERICAL 
HARMONICS METHOD! 


B. DAvISON 


ABSTRACT 


In applying the spherical harmonics method to multilayer problems it is neces- 
sary to invert certain matrices. It is shown, for the cases of plane and spherical 
symmetry, that these inverse matrices can easily be written down explicitly, so 
that there is no need for numerical matrix inversion. 


1. INTRODUCTION 


In the present note we are concerned with neutron transport problems in 
one-velocity theory for multilayer systems with plane or spherical symmetry. 
One of the more powerful methods for dealing with neutron transport theory 
is the spherical harmonics method. For the case of a single medium its theory 
is well known (Mark 1944, 1945; Davison 1956), but the technique of its 
application to multilayer problems has not yet received sufficient attention. 
An important and, unless properly handled, the most laborious part of the 
calculations in applying this method to such systems is to express the ‘‘arbi- 
trary’ constants (constants of integration) for each consecutive layer in 
terms of those for the layer from which one starts the calculation. If the 
constants of integration for one layer are known, those for the next layer can, 
of course, always be obtained by solving a system of (N+1) simultaneous 
algebraic equations, where NV, assumed odd, is the order of approximation. It 
turns out, however, that solving this system of (V+1) equations numerically 
is unnecessary, since their solution is always given by a simple explicit ex- 
pression. In other words there exists a simple explicit expression for the inverse 
of the ‘fundamental matrix’’. This fact is already known for problems with 
plane symmetry (Davison 1956). The purpose of the present note is to show 
this for the case of spherically symmetric assemblies. 

For the reader’s convenience we first summarize, in Section 2, the main 
relevant formulae of the spherical harmonics method and recapitulate, in 
Section 3, the derivation of the inverse of the fundamental matrix for the plane 
case. The derivation of the inverse of the fundamental matrix for the spherical 
case is given in Section 4. Finally in Section 5 we illustrate the application of 
the formulae by considering the critical size problem. 

For the sake of simplicity we restrict ourselves in what follows by the fol- 
lowing assumptions: 

(A) Scattering is isotropic in the laboratory system of coordinates. 
(B) Gaps and layers of strongly absorbing material are treated in the same 
way as any other layer. 


(C) There are no sources. 


1Manuscript received September 10, 1956. : 
Contribution from the Computation Centre, McLennan Laboratory, University of Toronto, 


Toronto, Ontario. 
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(D) The number of secondaries per collision is not equal to unity in any of the 
layers. 

The modifications needed if any of these assumptions is violated are easily 
found. In so far as (C) and (D) are concerned, the modifications do not in any 
way weaken the results derived below. In particular if (D) is violated, the 
number of secondaries per collision, c, in any layer being unity, one has simply 
to carry out the analysis for c ¥ i and pass to the limit with c tending to unity 
in the result. 


2. THE BASIC FORMULAE OF THE SPHERICAL HARMONICS METHOD 


In the Py-approximation spherical harmonics method, in both the plane 
and the spherical case, the angular distribution can be expressed as 


(2.1) ¥(E, uw) = goede Ont1) vale) Pals) 


+4(N+1 


) 
(2.2) vat) = DY {M(r/l,c)}iA; 


jul 
where, in the plane case, r is to be replaced by the Cartesian coordinate x, u is 
the cosine of the angle between the x-axis and the direction of neutron travel, 
and 


(2.3) {M(r/l, c)}a = {My (x/l, c)}a = Gr(vs, c) exp(rx/l), 


while, in the spherical case, r is to be replaced by the distance r from the center 
of symmetry, yu is the cosine of the angle between the radius vector and the 
direction of neutron travel, and 


(2.4) {M(r/l,c)}Z = {Men(r/l, c)}2 = Galvy, 6) ~/ —2lv4/ wr Kary (—vyr/l). 


Throughout the above formulae the order of approximation N is assumed 
odd, P,x(u) are Legendre polynomials, c is the number of secondaries per 
collision, / is the total mean free path, the A; are constants of integration, and 
the functions G,(», c) are defined by the recurrence relationships 


2n+1— 
v 


(2.5) — (n+1) Gara) + Cet Gy, c)-+nGes(r, ) = 0 


(don being the Kronecker symbol) together with 
(2.6) Go(v,c) = 1. 


The »; are the roots of the equation 

(2.7) Gwi(v;,c) = 0. 

These roots come in pairs: if vy = »; is a root of (2.7), so is y = —»,;, and for 
c # 1 they are all different. They will be assumed to be numbered in accor dance 
with the conventions 

(2.8) vij~=— rj, Via > 977 


together with 
(2.9) sign v; = sign j for v7 > 0. 





- i ee 
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Finally K,(z) is the modified Bessel function of the second kind; the Macdonald 
definition of these, 

(2.10) K.(z) = $n {I_,(z) —I,(z)} cosec tx, 

being used. We shall have to use later the following properties of these func- 
tions (Whittaker and Watson 1950) : 


(2.11) Kiyi1(z) = 2tK ,(2)/2+Ki-1(2), 

(2.12) K_,(s) = K,(z) = \/ 4/22 exp(—2), 

and 

(2.13) K,(z) —exp(mit)K,(—z) = wt exp(wit)I;(z) . 

It also follows from the definitions of G,(v, c) that 

(2.14) Gr(—yv,c) = (—1)"Galy, c). 

The A; are to be chosen so that, as we pass from each layer into the next, 
(2.15) ¥n(f) is continuous (nm = 0, 1,... NV) 


where y,(r) is interpreted as Wn, (x) in the plane case and as Wn spn(7) in the 
spherical case. 

When considering several layers at once we shall write cy, 4, vjx, and A; for 
the values of c, /, v;, and A; in the kth layer and shall denote by r;(x, or 7,) the 
value of r(x or 7) at the interface between the kth and the (&+-1)th layer. 
Combining (2.2) and (2.15) we readily see that where the A;,4; are known, 
the values of A, are given as the solution of the system of equations 


+4(N+1) +4(N+1) : 
(2.16) 2d, {M(ri/le, cr} n= 2, {M (fe/lesr, Ce+1) }nA tnt 
j= 


and the solution in question can be obtained at once, provided for our / and ¢ 


we have constructed {M-'(r/l, c)} such that 
N 


(2.17) LD {Mi (r/l, c)\{M(e/l, c)}a = 84s. 


n=0 
This question will occupy us in the next two sections. 
3. PLANE CASE 
We start by considering the quantity 
(3.1) Sa = (m+1) [Gm4i(? ue+t, Cu+.1)Gm(Yjx, Cr) —Giii( Pixs Cre) Gm(¥ e415 Ce+1)] . 


Using the recurrence relationships (2.5) to eliminate the Gm4i we easily 
obtain, provided m > 1, 
Sm = (1/rj——1/v 441) (2M+1)Gm( vin, Ce)Gm(V setts Ce+1) + Sm, 


and repeated application of this result gives 
Sn = (1/vp—1/v 41) a (2n+1) Galva, Cr) Ga(¥ e+1) Cri) + So. 


Putting m = N in this, noting that by (2.7) Sy = 0, and substituting for So 
from (3.1), with Go and G; given by (2.6) and (2.5), gives 
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N 


(3.2) D> (2m+1—crdon) Ga(v se, Ce) Gu (Penta, Cet) = Ve (Ce—Cet1)/(Vnsi—P ye): 


Recalling that under the assumption (D) all vj, for a given k are different, 
passing to the limit with c,.4: tending to c,, and noting that vg4, will then tend 
to vm, we arrive at (on suppressing the subscript k) 


(3.3) 2d (2n+1—Ccbom) Ga(vs, ¢) Galvs, c) = —vb4,/(dv,/dc). 


Comparing this identity with (2.3) and (2.7) we see at once that 
(3.4) {Myr (x/1, c)}) = — (dvj/dc)(2n+1—cbon) Ga(vs, ¢) exp(—vyx/l)/v, 


so that the solution of (2.16) for the plane case is given by 
+4(N+1) N 


Ay = x Aas 2, {Myi- (xn/les Cx) Vf Moi (Xe/lev1, Cx+1) }a 


toot 


or, on using (3.2) to simplify the inner sum, by 


(3.5) Ay = (Ce—Ce41) (dv »,/dcy) exp(—rpXx/lz) 
+4(N+1) 
x 2, A epi CXP (¥ n41%e/L) / (Vpr—vess)- 


We may note that (3.4) can be rewritten, on using (2.5) to eliminate the 
explicit appearance of 1/»;, as 


(3.6) {Myr (x/1, c) = { (M+1)Gusi (vy, €)-+NGa-1 (v5, €) jexp(—vyx/l) (dv,/dc). 


4. THE SPHERICAL CASE 
4.1 Proof of an Identity 
Analogy with (3.6) suggests that in the spherical case { Mgpn~"(r/I, c)}, may 
be proportional, apart from some simple factors, to 
((#+1)Gn4i(v;, €)Kn4a/2(vjr/l) +nGra(vj, ¢)Kn4(vjr/l)](dv;/dc) . 


In order to establish this result we shall consider an auxiliary quantity 
a N = = = - 
(4.1) To = (29/4) J vv 2 (Gp Kn4yGn—1Kn4+GarKnyG, ~1Ky-4)n 


where we have put for brevity 
Gr = Gal vj, Ck); G, = Gal iegr, Ce+1), 


(4.2) Kast: = Knsay(—yvjx), Kas Po Knsy(yv nti), 
y = re/ly. 


We introduce also the further auxiliary quantities 7,,(m > 1) defined as 
follows: Tm is obtained by replacing, in Tm1, Ky-misn and Ky—misp by 
Ky-m-—4 and Ky_m_4 respectively, while leaving all other factors in each term 
exactly as they were in 7-1. 

Let us assess now the actual numerical difference between T,, and Tn-1. 
Notice that Ky4, and Ky4; each appear only once in To, namely in the term 
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with n = N; Ky-4 and Ky-4 each appear twice in 7), namely in the terms 
with » = N—1 and n = N; Ky-sp and Ky_3, each appear three times in 
T2, namely in the terms with n = N—2, n = N—1, and nm = N; and so on. 
Notice also that, in To, Ky44, apart from other factors, is multiplied by Ky-+4, 
while Ky+4 is multiplied by Ky_;; in 7), apart from other factors, Ky-4 is 
multiplied on each occasion by Ky-3,2 while Ky_ is multiplied on each occasion 
by Ky-s, and so on. Thus we have 

(4.3.1) %—-T) = (Ky-32—Kw44)Kn-yNGyiGyt+ “adjoint terms”, 


(4.3.2) T2—Ti = (Ky-se—Ky-4)Kn-32[NGy+(N—1)Gy-2]Gy- 
+“‘adjoint terms”’, 


(4.3.3) T3:—T2 = (Ry-12—Ky-sp)Kn-sn{ NGyiGy 

+ [(N—1)Gy_1+(N —2)Gy_s]Gy_2} + “adjoint terms”, 
and so on, where ‘‘adjoint terms”’ is understood on each occasion to mean the 
expression similar to that written out explicitly, but with the barred and non- 
barred quantities interchanged. 

Recalling the abbreviations (4.2) and using the recurrence relationship (2.11) 
and (2.5) together with (2.7) one can readily see that the term written out 
explicitly in (4.3.1) is equal to 

(2N— 1)Ky-4Ky4(2N+ 1 —crdon)GuGu/Yv ikgiVik 
while the other term (the ‘‘adjoint term’’) gives 
—(2N—1)Ky-yRy-4(2N +1 —ce4150n) GuGw/Yr inv ines 
so that 
T,—T) = (2N- 1) Ky—yKy—4(Ce+1— Ce) Gu Gu Son /Yv jx? ins 

Transforming in a similar manner the formulae (4.3.2), (4.3.3), etc. we obtain 

T.—T, = (2N —3)Ky-3Kw-sa(ce+1— Ce) Gy—1Gy—150n—1/ YP jR¥ ik, 

T3;—T2 = (2N—5)Kw-spKw-s(Ce41— Ce) (GuGw Son + Gy—2Gn—250n—2)/VV 5K? ik-+1s 

T.-T3 = (2N —7)Kw-712K-10(Cet1—) (Gy-sGy—150n—1 
+Gy—sGy—s5on—s)/VjEY ik41, 


and so on. 

Examining these formulae one can readily see that (Tm41— 7m) vanishes for 
all m < N and hence 
(4.4) To =@T, =...= Ty. 


On the other hand it is evident from the construction of Tm that, of various 
Kn44 and Kn44, Ty can involve only Ky, K_4, Ky, and K-_, so that, in view of 
(2.12), 


Ty = (2y/4)\/—v pv asi KyRy > 2(GpGr-1+GnGa_1) 


N 
=I 


N 
= exply(ra—var)] 2 Ga{ (n +1) GaritnG,-1}. 
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Combining this with (4.4), recalling (4.2), and using (2.5) and (3.2), we can 
rewrite the last result as 


(4.5) To = exp[y(vje— vik+1)] (Ceti — Ce) / (ine — Mik), 
or substituting for To, etc., from (4.1) and (4.2), and slightly rearranging the 
resulting expression, we have 
(4.6) > Gav jes Ce) Knsg (ter ge/Tle) { (+1) Gata (Vest, Cot) Kats 2 (rev ats /le+1) 
+nGr-i (Vint, Crp) Ky (rev est /Le+1)} (2re/4)V (— V je t+1/Ueles1) 
= exp {re (¥je/be—v ts /let1) } (Cori — Ce) / Vari — Ye) 
+ terms vanishing for },41 = J. 


4.2 The Final Inversion Formula 


Passing in (4.6) to the limit with /,41 and c¢,4: tending to , and c, respectively, 
suppressing the subscript k, and multiplying through by dv;/dc we get 


(4.7) (2r/ml) (dv,/dc) Vand, Gu(vy, C)Kn+y(—vyr/l) 
X { (a +1) Gari (vs, C)Knys2(vir/l) +nGu_i(vs, C)Kny(vir/l)} = 815, 
while comparing this with (2.4) and (2.17) we see that 
(4.8) {Magn (1/1, €)} = (r/1) "(dv /de) \/ 2k / wr { (n+ 1) Garr Py €)Knsare(r7/2) 
+NGy-1(¥ 4, C)Kn-y(vyr/l)} 


and the solution of (2.16) for the spherical case is then given by 
+}(N+1) 


(4.9) Ayn = (27, / mlz) WV lets x/1x (dv p/dex) a A sett V — Pay 


i=t 


N 
x a n{Ga(vyr, Cr) Kn+4 (vere /l)G, ~1(V eps, Cus) Kany (—v eegite/le+1) 


+G,-1 (Vsx; Cr) Kn-y (Vpate/lk)Gn (Vas, Cet) Knsy ( i Vieeite/Te+1) } . 


4.8 Miscellaneous Remarks 

If /.41 is different from /,, no appreciable simplifications arise in (4.9). But 
if Apu, = , ie. if the kth and the (k+1)th media differ only in the number of 
secondaries per collision, their total mean free paths being the same, (4.6) 
shows that (4.9) reduces to 


(4.10) Ay = (Cx—Ce+1) (dv /dcy) exp(—vyre/lx) 


+4(N+1) 
xX 4, A et CXP (Ving ite/le) / (¥e—¥ 41) 


(spherical case with /,4: = ) which is identical with the corresponding 
formula, (3.5), for the plane case. 

Two remarks may be added in connection with this last result. Suppose the 
mean free path is the same throughout a spherically symmetrical system. 
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Then comparison of (4.10) with (3.5) shows that the Aj, (all 7 and e) for this 
system will be equal to those for the appropriately chosen plane case. Also 
comparing (2.3), (2.4), and (2.12) we see that 


r{Mepn(r/l, ¢)}o = {Mp(r/I, c)}3. 


Then, by (2.2), we have 
rWo,spn(7) = ¥o,pi(7) 


provided the mean free path is the same throughout the system. 

This result for the exact solution has been known for some time (Davison 
1956). We have now shown that it also holds in any odd-order approximation 
in the spherical harmonics method. The other remark arising in connection 
with (4.10) is as follows. Let Mp,(x/I, c), Mpi(x/l, c), Mapn(r//, c), and 
M,pn~(r/l, c) be the matrices whose elements are given by {M,(x/I, c)}Z, 
{ Myi7(x/l, c)}, etc. Then the result of comparison of (4.10) with (3.5) can 
be expressed as 
(4.11) Mopn—"(re/Lis Ch) Mi spn(re/les Ces) = Moi '(re/les e)Mopi(re/ley Ces) 


for l;, = Vest 


and this suggests that M,,,(r//, c) should be of the form 
(4.12) Mopn(r/l, c) = W(r/1)M,\(r/I, c) 


where W(r//) is a certain matrix whose elements depend on r// only. The 
formula (4.12) can be proved directly but the proof is rather lengthy and we 
shall omit it. 
5. CRITICAL SIZE PROBLEMS 

To illustrate the application of the above formulae we derive below the 
determinantal equation for the critical radius of a sphere in a multilayer 
reflector. We start, for the sake of simplicity, with a single-layer reflector 
contiguous to the core and extending to infinity, and shall distinguish the 
quantities appertaining to the core and the reflector by the subscripts 1 and 2 
respectively. Since the reflector extends to infinity, we should obviously take 
co < 1 and then, with the assumption (D) and the convention (2.9), the con- 
dition that there should be no supply of neutrons from infinity can be ex- 
pressed as 
(5.1.1) Aj. = 0 forj > 0. 

At the same time, using (2.2) and (2.4), one can readily see that the con- 
dition of regularity at r = 0 implies 
(5.2) AptA-p =0, j=1,2,...3(N+1). 
Expressing here Aj; and A_;; in terms of Yn.spn(7) (w = 0,1...) by means 
of (4.8) and recalling (2.8), (2.14), and (2.13) we easily obtain 


(5.3) > Vavepn (7) { (m+1)Gasi(van, €1) Lng-3/2(7¥:/l1) 
+nGy_i (vin, C1) [ny (7¥10/h1)} = 0 
(r in the core, i = 1,2... 3(N+1)). 
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Applying this at the interface between the core and the reflector, i.e. at r = 1, 
say, where 1; is the critical radius of the sphere, recalling (2.15), and using (2.2) 
and (2.4) for the reflector, together with (5.1.1), we easily obtain 


N 
x A 2V —lvye/n 2» G,(v 32, C2) Knsy(—vyor1/le) 


X { (m+1)Gayi(vaa, 61) Ing32(riva/h) +nGva, C1) Ln-y(rivir/hi)} =0 


or, which is the same, 
N 

x A_» Jlvp/n dX (—)"n{Ga(vyo, C2)Gu—1(Via, C1) Knog (11 2/l2) Iny (rir /hr) 
J n= 


—G,-1(v32, €2)Ga (Via, €1)Ky_y (1710 j2/l2) Ing (rv /h) } = 0. 


Eliminating between those equations the constants A_j;2: we obtain the deter- 
minantal equation 


Xu X12 eee X yw+1) 
(5.4) X21 X22 eee X24w+) = 0 


Xyw+1 Xiywen2 see X ys pyvsy 


where 
N 
(5.5.1) Xi = d (—1)"{Gy(v 52, C2)Ga—1 (vir, C1) Kngy (119 2/ le) In (71¥11/h1) 


—G,-1(v 32, C2)Ga(va, 61) Ky_4 (1710 32/12) Tn44 (1a /lr) 2 
Suppose now we have a two-layer reflector, with the inner reflector con- 
tiguous to the core and the outer reflector, and the outer reflector extending to 
infinity. We shall distinguish the quantities appertaining to the core, the 
inner reflector, and the outer reflector by the subscripts 1, 2, and 3 respectively, 
so that (5.1.1) is now replaced by 
(5.1.2) A;; = 0 for 7 > 0. 
We proceed exactly as above, except that, before substituting into (5.3), we 
should express A 2 in terms of A;3 by means of (4.10). In this manner we arrive 
at the determinantal equation of exactly the same form, (5.4), as before, except 
that the X,; are now given by 
ENED A 
(5.5.2) Xy= Do 3 (vy | 
Smt c=ce 


1 
N 
xX x n’ {Gay (v 53, €3)Gur—1 (¥972,€2) Kne44 (rev 43/ls)Ky—y(—12vy'2/l2) 
+Gy—1(¥33, C3)Gy’ (v 52, C2) Ky-4 (rev y3/ls)Knr44(—1rovy2/le) } 


x > (—1)"n{G,(vy2, €2)G,-1(¥ 11, C1) Kas (rivy2/le) Insy (ria /hi) 
— n—1(Py"2, €2)G,(¥ a1, ai)K, 3 (rv y2/le) Ing (rva/h)} 
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and similarly for any larger number of layers. If any two adjacent layers 
differ only by the value of c, the values of / being equal, the formula (5.5.3) 
can be easily simplified using (4.6). 

We have been assuming that the outermost layer of the reflector extends to 
infinity, but this is no restriction since any finite system can be always re- 
garded as embedded in an infinite purely absorbing medium. The assumption 
(D) can always be eliminated, as pointed out in the Introduction, while the 
assumption (C) will be automatically satisfied in critical size problems. Thus, 
provided the assumptions (A) and (B) are satisfied,? the ultimate determinantal 
equation for the critical radius is given by a determinant of the order }(N+1) 
(rather than of the order (V+1) which one might have expected at first sight) 
irrespective of the number of the layers. It may also be noted that the quanti- 
ties (dv;/dc) will have to be known only for the intermediate layers, but not 
for the core, nor for the outermost layer extending to infinity. 
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THE EMISSION SPECTRUM OF AIF IN THE 
VACUUM ULTRAVIOLET! 


S. M. NaupE anp T. J. Huco 


ABSTRACT 


Two vacuum ultraviolet band systems of aluminum monofluoride, the 
B'Z+-X1Z+ and C'Y+-X'Z+ systems, are described. The rotational analyses 
of a few bands in each system are given and the rotational and vibrational 
constants of the X'Z*+ and B'X* electronic states are determined. Accurate 
T,-values for the All, B'Z+, C'2+, and D'A states are given. 


A. INTRODUCTION 


Recent studies (Rowlinson and Barrow 1953; Dodsworth and Barrow 1955; 
Naudé and Hugo 1953, 1954; Barrow and Rowlinson 1954) of the spectrum of 
the aluminum monofluoride molecule have shown the existence of a number 
of singlet and triplet electronic states. These studies, both in emission and 
absorption, covered the wavelength region from 1250 A to 8800 A, but, in 
spite of the fact that a number of band systems involving transitions to or 
from the ground state were observed, no rotational analysis of the ground 
state has as yet been made. This may be ascribed to the fact that the only 
band system involving the ground state which lies in a region easily accessible 
to high dispersion instruments is the A'II-X'Z+ system at 43950 cm—. The 
vibrational structure of this system is, however, of such a nature that even 
under high dispersion the rotational structure is not resolved because of the 
strong overlapping of successive bands in the sequences. The other systems 
involving the ground state all lie in the vacuum ultraviolet where only low 
dispersion data were available. 

The opportunity was generously offered to one of us to obtain high dis- 
persion spectrograms of the region below 2000 A on the 3-meter vacuum 
spectrograph of the National Research Council of Canada in Ottawa. In its 
fourth order this instrument resolved the fine structure of some of the bands of 
the B'S+-—X'Z+ and C'S+—X'Z+ systems sufficiently well to permit rotational 
analyses. 


B. EXPERIMENTAL 


The AIF spectrum was excited according to the method previously described 
(Naudé and Hugo 1953) and photographed on the 3-meter vacuum spectro- 
graph described by Brix and Herzberg (1954). The reciprocal dispersion was 
0.63 A/mm. Overlapping orders were eliminated by the use of a LiF prism in 
front of the spectrograph slit (Brix and Herzberg 1954). Wavelength standards 
were furnished by second order iron lines. The wavelengths of the iron lines 
were converted to vacuum wave numbers by means of Kayser’s ‘‘Tabeile 
der Schwingungzahlen”’ as corrected by Edlén (1953). 


*Manuscript received October 19, 1956. 
Contribution from the National Physical Research Laboratory, Council for Scientific and 
Industrial Research, Pretoria, Union of South Africa. 
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C. THE B'Zt+-X!12+ SYSTEM 


This system was first observed by Barrow and Rowlinson (1954) in their 
study of the absorption spectrum of AIF in the Schumann region. It forms the 
strongest absorption feature and Barrow and Rowlinson were able to observe 
five sequences, viz. the 2-0, 1-0, 0-0, 0-1, and 0-2 sequences. 

The system is emitted with high intensity from the hollow cathode dis- 
charge. The over-all intensity is comparable to that of the C'S+—X'Z+ system 
which is possibly the strongest emission feature in the Schumann region. In 
the emission spectrum six sequences of single-headed bands degraded to 
shorter wavelengths were observed. The band heads with their estimated 
intensities are listed in Table I. 


TABLE I 
BAND HEADS AND ORIGINS OF THE B!2+—X!2+ system OF AIF. 


Band origin in cm! 


Band Band head Calculated from Calculated from 

v’-v"’ Intensity in cm7! Observed band head Equation (1) 
0-3 1 51925.5 §1933.1 §1932.9 
1-4 2 §2012.7 52020.2 52020.2 
2-5 1 094.6 101.9 101.4 
3-6 0 170.5 177.7 176.6 
0-2 5 697.9 52706.7; 706.6 706.7 
1-3 4 775.8 784.55 784.3 784.5 
2-4 3 848.6 857.0 856.4 
3-5 2 915.2 923.4 922.1 
4-6 1 976.0 984.0 981.4 
0-1 5 53479.8 53489.9 53490.0 53489.9 
1-2 3 547.2 557.1 §58.3 
2-3 2 611.6 621.3 620.7 
3-4 0 669.5 679.0 677.1 
0-0 4 54270.5 54282.4, 54282.6 54282.4 
1-1 2 329.5 341.4 341.4 
2-2 1 383.1 394.7 394.5 
1-0 4 55118.9 55132.0 55134.0 
2-1 2 162.9 177.3 177.6 
3-2 1 201.2 215.2 215.2 
4-3 1 233.7 . 247.4 246.2 
5-4 1 258.6 272.0 270.6 
2-0 2 951.0 969.7 970.2 
3-1 2 979.6 997.8 998.3 
4-2 2 56001.8 56019.6 56020.0 
5-3 1 014.6 031.9 034.9 


With the available resolution only the 0-0, 0-1, 0-2, and 1-3 bands could 
be analyzed, because the other bands are either too weak or too strongly 
overlapped to permit a rotational analysis. The origins of all the bands were, 
however, determined from the observed band heads and either the observed 
or calculated B,-values using the equation 

Vneaa—¥o = —(By’+B,”)*/4(B,’ —B,”). 


The values for the band origins determined in this way are also listed in 
Table I together with the values determined directly for the analyzed bands. 
These values can be represented by the equation 
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vo = 54250.8,+866.60(v’ +3) —7.45(v’ +4)? 
—0.045(v' +4)?— 


801.95(v"’ +3) +4.70(0”’ +3)? 


The vo values calculated from this equation are included in Table I for com- 


parison. 


The wave numbers of the 0-0, 0-1, 0-2, and 1-3 bands are given in Table II. 





1-3 Band , 
Ahead = 1894.8 A 
R(J) P(J) 

52792.46* 
794.13* 
795.63 
797.37 
799.18 
801.07 
802.92 
804.96 
807.02 
809.16 
811.38 
813.75* 
816.04 
818.47* 
820.96 
823.55 
826.19* 
828.94 
831.69 
834.54 52778.64* 
837.48 779.33* 
840.40 780.00* 
843.55 780.86" 
846.67 781.75 
_ 782.67* 
853.16* 783.63* 
856.52* 784.82 
859.98 785.99 
863.45* 787.18* 
867.02 788.53 
870.68 789.88 
874.38 791.34 
878.30* 792.90* 
882.00* 794.41* 
885.98* 796.15 
889.92 797.98* 
894.08* ‘99.70 
898.11* 801.75* 
902.28 803.52 
906.54 805.70* 
910.88* 807.71 
915.24* 809.83* 
919.75* 812.16 
924.21* 814.43 
929.00 816.87 
933.69* 819.36 
821.89* 
824.51 
827.20 
829.97 
832.75 
835.57* 
838.65 
841.70 
844.73* 
848.06 


TABLE II 
WAVE NUMBERS OF THE LINES IN THE B!2+-—X!2+ system OF AIF 
0-1 Band , 0-2 Band 
Ahead = 1869.9 A Ahead = 1897.6 A 
R(J) P(J) R(J) P(J) 
53490.83* 
492.02*  53488.55* 52708.88*  52705.24* 
493.29* 487.54* 710.16* 704.27* 
494.68* 486.57* 711.61* 703.36* 
496.12* 485.68* 713.06* 702.55* 
497.55* 484.86* 714.59* 701.72* 
499.12* 484.10* 716.16* 701.17* 
500.73* 483.37* 717.84* 
502.42* 482.69* 719.60* 
504.13* 721.43 
505.90* 723.33 
507.77* 725.30 
509.71* 727.31 
511.69* 729.41 
513.74* 731.60 
515.87* 733.86 
518.04* 736.18 
520.23* 738.58 
522.56* 741.06 
524.93* 743.62 
527.35* 746.23 
529.82* 748.91* 
532.42* 751.69* 
535.01* 754.53* 
537.70* 757.42* 701.17* 
540.66* 760.41* 701.72* 
543.49* 763.44* 702.55* 
546.44* 482.69* 766.55* 703.36* 
as 483.37* 769.72* 704.27* 
= 484.10* 772.98* 705.24* 
555.41* 484.86* — 706.31* 
558.60* 485.68* 780.00* 707.42* 
— 486.57* 783.17* 708.61* 
565.17* 487.54* 786.80* 709.93* 
567.95* 488.55* 790.50* 711.28* 
571.29* 489.68* 794.41* 712.74* 
574.72* 490.83* 797.98* 714.23* 
578.14* 492.02* 801.75* 715.81* 
581.77 493.29* 805.70* 717.47* 
585.36* 494.68* 809.83* 719.19* 
589.07* 496.12* 813.75* 721.00 
592.79* 497.55* 818.05* 722.89 
596.55* 499.12* 822.27* 724.85 
600.31* 500.73* 826.48* 726.87 
604.23 502.42* 830.94 728.96 
608.24* 504.13* 835.32" 731.14 
— 505.90* 839.83 733.42 
616.55 507.77* 844.43* 735.75 
620.73* 509.71* — 738.15 
624.84* 511.69* 853.89* 740.60 
629.24* 513.74* 858.57* 743.16 
633.59 515.87* 863.45* 745.79 
637.90* 518.04* 868.45* 748.46' 
642.40 520.23* 873.47* 751.26* 
647.11 522.56* 878.30* 754.13* 
651.41* 524.93* 883.77* 756.99* 
656.34 527.35* 888.86 760.06* 
661.09 529.82* 894.08* 763.14* 
665.94* 532.42* 899.56 766.25* 
670.79* 535.01* 904.95 769.51* 
537.70* 910.55* 772.98* 


0-0 Band 
Mhead = 1842.6 A 
R(J) P(J) 
0 54284.06* 
1 285.26*  54281.75* 
2 286.56* 280.68* 
3 287.98 279.65* 
4 288.85 278.71* 
5 290.13 277.77* 
6 291.67 276.92* 
7 293.20 276.03* 
8 294.77 275.31* 
9 296.44 274.46* 
10 298.10 
ll 299.85 
12 301.67 
13 303.50 
14 305.38 
15 307.37 
16 309.38 
17 311.42 
18 313.55 
19 315.73 
20 317.92 
21 320.19 
22 322.57 
23 324.86 
24 327.34 
25 _ 
26 — 
27 334.99 
28 337.62 
29 340.20 
30 343.13 
31 345.97 
32 348.71 
33 351.68 274.46* 
34 354.59* 275.31* 
35 357.70* 276.03* 
36 360.62* 276.92* 
37 363.98* 277.77* 
38 367.00 278.71* 
39 370.12 279.65* 
40 373.50 280.68* 
41 376.91 281.75* 
42 380.19* 282.88* 
43 — 284.06* 
44 387.15* 285.26* 
45 390.65* 286.56* 
46 394.48* 287.98 
47 398.08* 289.39 
48 401.69 290.82 
49 405.50 292.28 
50 — 293.84 
51 413.09 295.46 
52 416.99 297.17 
53 421.16* 298.87 
54 425.05* 300.66 
55 429.13 302.50 
56 433.25* 304.41 
57 437.37 306.30 
58 441.69 308.30 
59 445.97 310.41 
60 450.30 312.52 
61 454.65 314.63 
62 459.10 316.87 
63 463.57 319.23 
64 321.54 
65 323.90 
66 326.39 
67 328.81 


*Lines marked with an asterisk could not be measured accurately. 
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The B, and vo values for these bands were.determined graphically from the 
equation 


R(J-1)+P(J) = 2v0+2(B,’ —B,”’) J?—2(D,'’ —D,") F?(J?+1). 


It was found that, within the present accuracy, the term in (D,’—D,’’) could 
be neglected for all four bands, i.e. D,’ = D,’’ = 9.1.X10-7 (Naudé and Hugo 
1954). In the case of the 0-0, 0-1, and 0-2 bands the accurate By’ value ob- 
tained from the 0-0 band of the F'II—B'=+ system was ‘1sed to determine B,’’. 
In the case of the 1-3 band the B” value was obtained from the 1-3 C'2+—X'Z+ 
band where the B,’ value is known accurately from the C'Z+-A'II system. 
The B, values for the X'Z+ and B'=+ states as obtained from the B'Z+— X'!Z+ 
and C'Z+— X'Z+ systems are given in Table III. The listed By value for X'Z* 


TABLE III 
OBSERVED AND CALCULATED B,-VALUES OF THE X!2+ AND B!Z+ staTEs OF AIF 
X'1Z+ state B'=+ state 
v B, in cm B, in cm" B, in cm B, in cm™! 
(observed) (calculated) (observed) (calculated) 
0 0.5493+0.0001 0.54987 0.5759s8+0.00005 0.57508 
1 0.5442+0.0009 0.54503 0.57127+0.0001 0.57127 
2 0.540;+0.0001 0.54020 
3 0.535,+0.0001 0.53537 
4 0.530;+0.0001 0.5305. 


is the average value from the two band systems. The B, values of the X'Z+ 
and B'=+ states can be represented respectively by the equations 


B, = 0.55223 —0.00483(v+ 3) 
and zB = 0.57834 —0.0047:(v+3) ° 


The values calculated from these equations are also shown in Table III. 


D. THE C'!2+-X'!Z+ SYSTEM 


This system was first observed in absorption by Barrow and Rowlinson 
(1954). In absorption the system is very weak but quite strong in emission 


TABLE IV 
BAND HEADS AND ORIGINS OF THE C!2+—X!2Z+ sysTEM OF AIF 


Band origin in cm=! 





Band Band head Calculated from Calculated from 
v’-v"’ Intensity in cm=! Observed band head Equation (2) 
1-4 2 55563.8 55569.6 55569.7 55570.5 
2-5 0 726.3 732.1 733.2 
0-2 3 56173.2 56179.9 56180.0 
1-3 3 327.6 56334.02 334.1 334.3 
2-4 2 481.3 487.7 487.6 
0-1 5 955.4 962.9 963.1 
1-2 4 57100.5 57107.8 57108.1 
2-3 3 244.7 251.9 252.0 
0-0 5 746.9 57755.64 755.4 755.6 
1-1 2 882.8 891.1 891.1 
10 3 58674.0 58683.6 683.7 
2-1 2 799.3 808.7 808.8 
3-2 0 923.5 932.6 932.9 
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from the hollow cathode discharge. The observed emission spectrum consists 
of the 1-0, 0-0, 0-1, 0-2, and 0-3 sequences. The bands are single-headed and 
degraded to shorter wavelengths. The band heads with their estimated intensi- 
ties are listed in Table IV. The 0-3 band is masked by the stronger bands of 
the 1-0 sequence of the B'2+— X'Z+ system. 

The 0-0, 1-3, and 1-4 bands were analyzed and the wave numbers of the 
lines are given in Table V. The rotational constants and origins of these bands 
were determined from the R(J—1)+P(J) data by using the accurate B,’ 
values of the C'Z+-A'II system (Naudé and Hugo 1953). The term in 


TABLE V 
WAVE NUMBERS OF THE LINES IN THE C!2+-X'!2+ sySTEM OF AIF 




































































0-O Band ,. 1-3 Band . 1-4 Band , 
J Mead = 1731.7 A Ahead = 1795.3 A Mead = 1799.7 A 
R(J) P(J) R(J) P(J) R(J) P(J) 
0 57756.82* 55570.74* 
1 758.07*  57754.52* 572.00*  55568.56* 
2 759.40* 753.48* 573.37* 567.58* 
3 760.82* 752.49* 56339.33* 574.86* 566.65* 
4 762.28* 751.69* 340.96 576.49* 565.95* 
5 763.88* 750.86* 342.54 578.18* 565.39* 
6 765.47* 344,27 580.00* 564.81* 
7 767.19* 346.12 581.86* 
8 769.01* 348.14 583.92* 
9 770.91* 350.10 586.01* 
10 772.83* 352.24 §88.22* 
1l 774.84* 354.51 590.45* 
12 776.94* 356.76 593.04* 
13 778.92* 359.17 595.62* 
14 781.18* 361.64 598.20* 
15 783.46* 364.26 600.99* 564.81* 
16 785.82* 366.97* 603.86* 565.39* 
17 788.28* 369.77 606.77* 565.95* 
18 790.82* 372.70 56329.59* 609.88* 566.65* 
19 793.42 375.63 330.17* 613.04* 567.58* 
20 796.07 378.70 330.97 616.31* 568.56* 
21 798.79 381.83 331.82 619.75* 569.55 
22 801.68 385.10* 332.82* 623.23 570.74* 
23 804.56 388.50* 333.82 626.80 572.00* 
24 807.61 391.78* 334.93 630.45* 573.37* 
25 810.68 750.86* 395.38* 336.18 634.19* 574.86* 
26 813.76 751.69* 399.04* 337.48 638.14* 576.49* 
27 816.97 752.49* 402.83 338.91 642.03* 578.18* 
28 20.24 753.48* 406.69 340.48 646.20* 580.00* 
29 823.64 754.52* 410.62 342.04 650.45 581.86* 
30 827.07 755.60 414.66 343.75 654.78 583.92* 
31 830.59 756.82* 418.90 345.65 659.28 586.01* 
32 834.16 758.07* 423.10 347.55 663.81 588.22* 
33 837.81 759.40* 427.46 349.51 668.38 590.45* 
34 841.57 760.82* 431.89 351.66 673.15 593.04* 
35 845.39 762.28* 436.41* 353.89 677.96 595.62* 
36 849.33 763.88* 441.08* 356.20 682.96 598.20* 
37 853.26 765.47* 445.68* 358.62 688.02 600.99* 
38 857.30 767.19* 450.46 361.16 693.13 603.86* 
39 861.49 769.01* 455.39 363.76 698.49 606.77* 
40 865.67 770.91* 460.32 366.45 703.85 609.88* 
41 869.94 772.83* 465.44 369.18 709.43 613.04* 
2 874.32 774.84* 470.66 372.07 714.94 616.31* 
43 878.73 776.94* 475.90 375.07 720.66 619.75* 
44 _ 779.26* 378.22 
45 887.79 781.48* 381.34 
46 892.47 783.78* 
47 897.24 786.26* 
48 902.09 788.77* 
49 907.05 791.39* 
50 911.99 794.08 
51 916.97 797.77 
52 922.16 799.65 
53 927.34 802.57 
54 932.65 805.58 
55 937.97 808.65 
56 943.47 811.80 
57 815.08 
58 818.37 








*Lines marked with an asterisk could not be measured accurately. 
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(D,’—D,"’) could be neglected again for all the bands. This gives D,’ = D,” 
= 9..X10-7 (Naudé and Hugo 1953). The B;’’ and B,” values are given in 
Table III and the » values in Table IV. As for the B'2+—X'Z+ system the vo 
values of all the observed bands determined from the band head data and B, 
values are compared in Table IV with the values calculated from the equation 
representing the band origins, viz. 

(2) vo = 57687.54+938.39(v' +4) —5.18(v’ +4)? —801.95(0"’ +4) +4.70 (vo +4)’. 


E. DISCUSSION 


The rotational and vibrational constants of the X'Z+ and B'Zt electronic 
states obtained in this study are summarized in Table VI. 


TABLE VI 


VIBRATIONAL AND ROTATIONAL CONSTANTS OF THE 
X'2+ anp B'Z+ states or AIF 


Constant X1z+ Birt 

B, (cm) 0.55223+0.0001 0.5783,+0.0002 
a, (cm~) 0.004s3+0.0001 0.0047:+0.0001 
D (cm~) 9.41077 9..X 1077 
I, (g-cm.?) 50.663 X 10~4° 48.38; X 10~*° 
r. (cm.) 1.6547 10-8 1.617 9X 10-8 
Yoo (cm) — 54282.4, 
T. (cm~) 0 54250.8, 
we (cm) 801.95 866.60 
weXe (cm™) 4.75 7.45 
WeVe (Cm~) —_ 0.045 





It is interesting to compare the wave number interval between vo9C'Z*+ and 
vooB'Z* as obtained from the vacuum ultraviolet bands with that obtained 
from the visible and infrared systems F'II-A'II, C'Z—A'll, and F'II-B'2, for 
this should give an indication of the absolute accuracy of measurement 
attained in the vacuum ultraviolet. These two determinations give respectively 
the values 3473.1, cm and 3473.24 cm~. In the determination of the last 
value the Edlén corrections to the Kayser table were taken into account. 
From these values it seems probable that an absolute accuracy of 0.1 cm™ is 
attained in the vacuum ultraviolet. The relative accuracy is appreciably higher. 

From the data presented here and those previously published (Naudé and 
Hugo 1953, 1954) accurate 7, values can be assigned to the A'Il, B'Z, C'S, 
and D'A electronic states. These values are respectively : 43947.99, 54282.4s, 
57687.54, 61227.0o. 
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THE SPECTRUM OF SILICON HYDRIDE! 


A. E. DouGLas 


ABSTRACT 


The 0-0, 1-0, 2-1, and 2-2 bands of the ?A-*II band system of SiH have been 
photographed under high dispersion and analyzed. The vibrational and rotational 
constants have been determined and the dissociation energy has been found 
to be 3.19+0.25 ev. A list of SiH lines which may occur as interstellar absorption 
lines is given. 


INTRODUCTION 


A band system due to a *A-"II transition of silicon hydride has been known 
for many years and has been photographed under high resolution by Jackson 
(1930) and by Rochester (1936). The analysis of the 0-0 band has been given 
by Jackson and somewhat modified by Mulliken (1931). Rochester has given a 
detailed analysis of the 1-1 band of SiH and the 0-0 and 1-1 band of SiD. 
Recently during an investigation of the Siz molecule (Douglas 1955), it was 
found that SiH bands were emitted strongly from the same discharge tube. 
Spectra from this source which were photographed at a dispersion of 0.35 
A/mm. showed the 0-0, 1-1, 2-2, 1-0, and 2-1 bands of the 2A—*II band 
system. No other band systems which could be attributed to SiH were found. 
Because the SiH molecule is one of considerable astrophysical and chemical 
interest, it was considered desirable to measure and analyze these new bands 
of SiH. This paper gives the measured wave numbers of the lines and the 
molecular constants derived from them. 


RESULTS 


The 0-0, 1-0, 2-1, and 2-2 bands have been measured and analyzed. This 
analysis shows that the analysis of the 0-0 band given by Jackson is essentially 
correct but, because of lack of resolution on his plates, a number of lines of 
low J value were incorrectly assigned. The measurements and the analysis 
of the 1-1 band given by Rochester were found to be quite satisfactory and 
this band was not remeasured. The wave numbers of the lines are given in 
Table I. 

The analysis of the bands was quite straightforward and all of the necessary 
equations have been given by Rochester. The constants which have been 
determined are given in Table II. For the *II state the small term vy in the 
spin splitting has been neglected. For the ?A state, the constants A and y 
were found to be 3.83 and 0.9 cm respectively; values which are very close 
to those given by Rochester. The experimentally observed A doubling cannot 
be represented accurately by the approximate equations of Mulliken and 
Christy (1931) and no improvement can be made on the A splitting constants 
p and g given by Roghester. 
~ Manuscript received September 27, 1956. 
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TABLE II 
ROTATIONAL AND VIBRATIONAL CONSTANTS OF SiH (cmM™) 


































0-0 Band 1-0 Band 2--1 Band 2-2 Band 
cys 7.391 7.390 Tilted 6.961 
p” 3.951074 3.941074 4.1x<10- 3.8X10-4 
B’ 7.287 6.858 6.347 6.347 
D’ 5.23 10-¢ 6.011074 7.31074 7.26X10-4 
vo 24193 .05 25853 .70 25345. 20 23445. 40 
a” 143.35 143.35 143.81 144.28 
war 042.47, wex,”’ = 35.67, we.’ = 1858.13, wx.’ = 98.73, By,” = 7.498—0.214(v+}), 


470 —0.347 (v+4) —0.041(v+})?, Te = 24393.15, rv.” = 1.521, r,’ = 1.523. 


a 
T 
Sil 


DISCUSSION 


For SiH and related molecules, the correlation of the ?II and 2A states 
with the states of the separated atoms has been given by Mulliken (1932). 
The 2II state of the molecule, which is the ground state, dissociates into the 
ground states of the atoms. The ?A state of SiH dissociates into the ground 
state of hydrogen and the 'D state of silicon. Thus, if the energy necessary to 
dissociate the 2A state of SiH can be determined, then this value can be com- 
bined with the excitation of the 'D state of Si and the yoo value of the 7A—II 
band system to give the dissociation energy of SiH. 

It is evident from the constants given in Table II that the 2A state is a 
shallow one. An extrapolation of the two observed vibrational quanta to the 
dissociation limit gives a value of 7840 cm™ for the depth of the ?A potential 
curve. The error in this number may be quite large, perhaps as great as 
+2000 cm-!. The dissociation energy of the *II ground state of SiH, corre- 
sponding to this value for the A state, is 25740+ 2000 cm™. Thus the relatively 
crude estimate of the depth of the A state leads to quite a good estimate of 
the dissociation energy of SiH. The dissociation energy of SiH is very nearly 
the same as that of the similar molecule CH. 

The SiH molecule probably occurs in interstellar space and it is possible that 
interstellar absorption lines of SiH may be found. The wavelengths of the 
possible absorption lines are given in Table III. Of the lines listed in this table, | 
the Ro:(3) line at 4119.482 A should be much the strongest since it is the only 


TABLE III 
POSSIBLE INTERSTELLAR ABSORPTION LINES OF SiH 












Transition v (cm!) » (A) 
R2a(}), 0-0 24268 .02 4119. 482 
Qal(Z), 0-6 247 .04 4123 .052 






247 .23 4123 .022 


Ri(3), 0-0 252.99 4122.040 
253.18 4122.008 


R2(3), 0-0 292.02 4115.417 
292.22 4115.384 


Ra(3), 1-0 25928 . 20 3855 .711 
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line of the 0-0 band arising from the lowest rotational level and the 0-0 band 
is much stronger than the 1-0 band. This line occurs on the published spectra 
of both Jackson and Rochester. Jackson has given an incorrect assignment 
to the line and Rochester has not assigned it. It is interesting to note that for 
the corresponding electronic transition of CH it is not the first Re; line but the 
first Re line which arises from the lowest rotational level of the *II state. This 
difference is quite in accord with theory (see Herzberg 1950, page 234) since, 
for SiH, the coupling of the *II state is case (a) while for CH it is case (6). | 

A short search for other band systems of SiH failed to show any. The 0-0 
band of the *A—*II system is very strong and it appears that, with the source 
used here, any other bands which do occur are less than one-hundredth as 
intense as the 0-0 band. 


I wish to thank Dr. G. Herzberg for his interest in this work and Miss L. 
Howe for carrying out much of the measurement and computation. 
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THE FREEZING POINTS OF HIGH PURITY METALS AS 
PRECISION TEMPERATURE STANDARDS 


I. PRECISION MEASUREMENTS WITH STANDARD 
RESISTANCE THERMOMETERS! 


E. H. McLAREN 


ABSTRACT 


The techniques and difficulties encountered in measuring temperatures to the 
highest precision with platinum resistance thermometers are discussed. It is 
shown that the relative drift of the resistance coils in the Mueller resistance bridge 
used for these measurements is less than a part per million per year. The intrinsic 
resistance of a platinum thermometer is comparatively unstable, and results 
showing some effects of cold work and heat treatment on several thermometers 
are given. 

As each precision temperature determination involves the resistance of the 
thermometer at the triple point of water, extensive measurements have been 
carried out to obtain information on: (a) the reproducibility of temperature in 
particular cells, (b) the variation in temperature among cells, and (c) the long 
term stability of cell temperatures. 

The limiting uncertainties in temperature measurements due to variation in 
the — the thermometers, and the triple point cells are each of the order of 
IoC: 


INTRODUCTION 


From —183° C. to 630° C. the International Temperature Scale is defined 
in terms of standard resistance thermometers that have been calibrated at 
the oxygen, ice, steam, and sulphur points; excellent reviews on temperature 
measurement in this range have been given by Mueller (1941), Stimson (1955), 
and Hall (1955). At the upper end of the range the limit of accuracy in tempera- 
ture determinations is set by the degree of reproducibility of the sulphur 
boiling point; unfortunately the high degree of pressure dependence (10-*° C. 
per micron Hg) and the considerable period, of the order of a day (Stimson 
1955), that is required to establish temperature equilibrium adversely affect 
this reproducibility, which is considered to be of the order of 5X10-*° C. 
(Hall 1955). 

The freezing points of highly purified samples of various metals are free 
from these particular disadvantages; in this laboratory experiments on tin, 
cadmium, and zinc have shown pressure dependences of less than 10-*° C. per 
cm. Hg and, with appropriate techniques, freezes with a duration of one to 
two hours that are readily reproducible to the order of 2X10-4° C. These 
results suggest that the zinc point (419.5° C.) could usefully replace the sulphur 
point (444.600° C.) on the International Temperature Scale and they also 
provide a means of testing, to this order of accuracy, the validity of the 
interpolation formulae used for platinum resistance thermometers. Details of 
these experiments will be presented in Part II of this paper. 

Temperature measurements to the order of 10~*° C. are not only dependent 
on the stability and homogeneity of the measured environment, however, but 

1Manuscript received September 24, 1956. 
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demand very great care in the selection and use of the measuring equipment . 
Part I of this paper describes the techniques that have been employed and the 
limitations that have been determined at this laboratory in the use of resistance 
bridges, platinum resistance thermometers, and water triple point eells. 


RESISTANCE MEASUREMENTS 

Equipment 

The resistance measurements reported here were made on a Leeds and 
Northrup G-2 Mueller Bridge (Mueller 1916; Mueller 1941; Stimson 1955) 
thermostatted to plus or minus one hundredth of a degree at 35° C. and having 
an extra 10-° ohm decade (LeLacheur 1952); a worm drive is used on the slide 
wire connecting the two ratio coils to give fine control on the bridge ratio 
adjustment. Extension leads to the thermometer junction boxes are No. 12 B 
and S gauge copper conductors terminating at copper binding posts and laid in 
half inch B X conduit. The detector is a Leeds and Northrup type HS galva- 
nometer mounted on a Julius suspension in a draft-free compartment and 
having a nominal sensitivity of 0.017 uv./mm. at the 5.5 meter optical path 
distance employed. The equipment is located in an air-conditioned room which 
is semidarkened for taking measurements. 


Stability 

1. Bridge zero 

For high precision measurements it is necessary to determine the correction 
to true zero resistance for the bridge. In Leeds and Northrup Mueller bridges 
the Wheatstone type network is trimmed so that a small positive excess of 
resistance remains in the X arm when the thermometer is shorted out; this 
off-zero resistance is measured and subsequently subtracted from every 
resistance determination. The measurement is made by galvanometer de- 
flections using a branch current of 5 ma. and the uncertainty in the measured 
correction is less than one microhm. Values of the bridge zero correction over 
a two month period are shown in Fig. 1, the fluctuations being associated with 





MICROHMS 








Fic. 1. A series of bridge zero corrections during a 10 week period. 


slight changes in the resistance of connecting conductors and components that 
are outside the thermostatted section of the bridge. An error of 10 microhms 
in bridge zero will introduce an error of about 0.0002° C. in the temperature 
measured at the zinc point with a 25 ohm standard resistance thermometer. 


2. Bridge Ratio 

The bridge ratio, i.e. the ratio of the resistances of the two fixed 1000 ohm 
arms of the bridge, is closely monitored immediately after any precision 
resistance determination and the readings are corrected to unity ratio with an 
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uncertainty of less than one microhm. The ratio coils are mounted outside the 
thermostatted enclosure in the bridge and the greatest depariures from unity 
occur when room temperatures vary abnormally. Occasionally the ratio may 
drift as much as four parts in ten million in a day, which, without correction, 
would correspond to an error of nearly 0.0002° C. in a calculated zinc tempera- 
ture. 


8. Bridge Resistors 

The stability of the resistors in the bridge is of utmost importance in re- 
sistance determinations, and considerable time and effort must be spent in 
checking coil resistances when the greatest possible accuracy is being sought. 

Fig. 2 shows the results over a period of two months of monitoring a Leeds 
and Northrup 10 ohm standard resistor, thermostatted to +0.01 deg., against 


£ 0° oo 
uo 272.2 202, of 2° © 5? Gp Pol og a5 


4 Std. Dev. 1.542 
MAY 





Fic. 2. The results of a series of difference measurements on a standard 10 Q resistor and 
the first resistor in the 10 Q decade of the bridge over a period of 10 weeks. 


the 10 ohm coil in the first decade of the bridge. Each measurement was made 
with a branch current of 5 ma. in a single sequence of NRRN commutations. 
The standard deviation about the average resistance for the period is +1.5 
microhms and is accounted for by the uncertainties attendant on balancing the 
bridge in the four commutations and determining the bridge zero and bridge 
ratio corrections. Resistance is estimated to the nearest microhm in each of 
these operations. 

Systematic internal calibrations of the coils in the bridge have been carried 
out during the past two years, the 10-ohm coils about every month and the 
lower resistance coils at longer intervals. Results obtained from the comparison 
of individual 10-ohm coils against the sum of the 10 one-ohm coils in the bridge 
are given in Fig. 3. These results are given in terms of the bridge unit of re- 
sistance which is arbitrarily defined as one-tenth of the sum of the resistances 
of the 10 one-ohm coils in the bridge. 

Fig. 3 shows that the relative drift of the resistors is less than a part per 
million per year except for the abrupt discontinuity that appeared in June 1955. 
This discontinuity occurred during a general power failure of several hours’ 
duration which resulted in the section of the bridge which is normally thermo- 
statted at 35° C. falling to room temperature. It is possible that the resulting 
change in relative resistance of the 10-ohm coils, varying between 0.8 and 
2.4 parts per million, is due to work hardening of the manganin as the coil 
forms contracted. Provision has since been made for battery operated standby 
power to heat the enclosure. The conclusion drawn from these results is that 
the coils in the Leeds and Northrup bridge, when well aged and continuously 
thermostatted, have a stability comparable with that of the highest quality 
standard resistors. 


Serre ree 
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Fic. 3. Variation in resistance of individual coils in the 10 2 decade against the sum of the 
10 coils in the 1 & decade over a period of three years. 


The importance of careful bridge calibration is evident from the following 
considerations: Temperature is determined from the ratio of the resistance of 
the thermometer at the temperature under investigation to its resistance at 
the ice point. Each of these resistances on the platinum thermometer will be 
balanced by the sum of perhaps 25 coils in the bridge and as the calibration of 
each coil in the bridge has an uncertainty of about +1 microhm the probable 
error in the determination of the resistance of the thermometer is of the order 
of +5 microhms. These combined uncertainties in calibrated resistance values 
result in a probable error of about +0.0001° C. in a temperature measurement 
at the zinc point with a 25.5 ohm standard thermometer. 


STABILITY OF Rrp 


The equation relating the resistance of a platinum wire to the International 
Temperature Scale in the range 0° C. to 630° C. has the form 


R, = Ro(1+At+Be) or R; = (Rrp—AR)[1+At+Be] 


where AR = Ryp—Ro, A and B are temperature coefficients of the particular 
sample of platinum, R; is the resistance of the wire at temperature ¢, and 
Rrp and R, are the resistances at the triple point of water and the ice point 
respectively. For work of the Highest accuracy the resistance Rp on a standard 
resistance thermometer is computed from the value of the resistance at the 
triple point of water, the ice point by definition being 0.01 deg. below the 
triple point. As Ryp enters into every temperature calculation the stability of 
the intrinsic resistance of a thermometer and the reproducibility of the tem- 
perature attained in water triple point cells directly affect all temperature 
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determinations. Definitive information on the constancy of the coefficients A 
and B is not yet available but an investigation to determine how Rrp varies 
with handling procedure and an intercomparison of the temperature realized 
in several triple point cells has been made. 


(a) Cold Work and Heat Treatment on Standard Thermometers 

As the electrical resistivity of a metal is a structure-sensitive property 
(Broom 1954) the specifications of the so-called ‘‘strain-free’” standard re- 
sistance thermometer includes the terms ‘‘well-annealed’”’ (C.I.P.M. 1948). 
This requirement ensures the relief of strain introduced in the drawing and 
fabrication procedures, but unfortunately there remains the possibility of 
introducing additional resistance by accidentally cold working the ther- 
mometer in the course of its day-to-day use. 

Experience with 12 standard thermometers clearly indicates that even with 
very careful handling it is not possible to eliminate changes in resistance due 
to work hardening of the platinum. A gentle tap on the sheath near the coil of 
a Meyers thermometer (Stimson 1955) may add as much as 20 microhms 
(equivalent to 0.0002° C.) to its original resistance of about 25.5 ohms at the 
triple point of water. Subsequent annealing at 450° C. will remove part of the 
added resistance, the amount removed being dependent on the type of defects 
introduced into the wire; recovery from vacancy formation occurs more readily 
than recovery from defects produced by higher types of deformation (Broom 
1954). Annealing at temperatures much higher than 450° C. is impractical for 
thermometers on mica crosses in pyrex sheaths. 

The effect of this work hardening was studied on two standard thermometers: 
S144 made by Meyers and having a single layer bifilar coil about one inch long 
and TM-1 constructed at the National Research Council in the coiled helix 
form (Meyers 1932) with a coil length of about two inches. Each thermometer 
was first annealed at 450° C. until its resistance at the triple point of water 
went through a minimum, after which the cold work was applied by lightly 
striking the tips of the thermometers 40 times against the edge of a wooden 
table; they were then reannealed at 450°C. Fig. 4 shows how Rrp varied 
during the operation. Annealing removed 74% of the added resistance from 
thermometer S144 and only 13% from thermometer TM-1. 

In an investigation to determine the effect of cold working on the tempera- 
ture coefficients of platinum wire, three annealed thermometers were used to 
determine the Rz,/Rrp ratios on four zinc freezes. One of the thermometers, 




















a % = Q 
0-0----0 a 
a $00 "CA ‘ 
~ 20°C ANNEAL! 2 3000 OPC AEA. $e QeeOeeQ ow oe HO 
a ' - 
: “7s f 450°C ANNEAL 
z sor 44 (74% RECOVERY leo 
a ° = B14 (78% RECOVERY) = TM-1 (13% RECOVERY) 
w 200h 2 Nr aa 
% ° OO ie Oca cen ° 2 1000 
te 100 450°C ANNEAL & 
2 ° z ° 
= 2 4 6 2 4 6 8 10 12 os 2 4 6 2 oe 6 6 10 
DAYS HOURS DAYS HOURS 


Fic. 4. Effect on Rrp of cold working and annealing two resistance thermometers. 
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S144, was then work hardened at room temperature until 14,000 microhms had 
been added to Ryp; about 6000 microhms of this increase in resistance were 
removed by reannealing at 450° C. The Rz,/Rrp ratios were again determined 
with the three thermometers on four zinc freezes. Table I shows that the ratio 







TABLE I 














Before cold working S144 After cold working S144 

See eee erat SS ARrp, A(Rzx/Rrp), 
Rrp (b.u.) Rzn/Rrp Rrp (b.u.) Rz»/Rre p.p.m. p.p.m. 
$144 25.49243 2.567938 25.50081 2.567902 +329 —14 


$156 25.48280 2.568059 25.48280 2.568059 0 0 
S163 25.48099 2.567710 25.48099 2.567710 0 0 































for the strained thermometer was lowered by 14 p.p.m. while the ratios for 
the untreated thermometers remained constant, the estimated uncertainty in 
a Rz,/Rrp ratio measurement being less than 1 p.p.m. It appears from this 
experiment that cold working a thermometer permanently alters its tempera- 
ture coefficients, if subsequent annealing is restricted to a temperature of 
450° C.* In an effort to increase the temperature coefficients of thermometer 
S144 an electric current was used to heat the wire to 650° C. for one hour. This 
lowered the ratio Rz,/Rrp by an additional 16 p.p.m., presumably owing to 
active contamination of the hot platinum by the mica cross. 

Continual use of a thermometer at elevated temperatures also contributes 
to the steady increase of the annealed resistance at a given temperature. Fig. 5 
shows the variation of the resistance at the zinc point, Rz», for a Meyers 
thermometer that was immersed continuously for 8 days at about 420° C. 
















Mean Daily Rzy for Thermometer S I56 
Continuously immersed at ~ 420°C 






DATE No.of Freezes 
2i 5 





| —-waarad 









2 23 






October 1952 


Fic. 5. Increase in resistance at Rzn for a Meyers thermometer during heating. 










*Corrucini (1951) demonstrated much greater lowering of temperature coefficients of an- 
nealed platinum wires by rapid quench hardening; reannealing at temperatures between 400° C. 
and 800° C., he was able to regain the original coefficients within the precision of his measure- 
ments, i.e. to within the order of 100 p.p.m. 
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Rrp increased proportionately, at the rate of 10 to 20 microhms per day, so 
that no measurable changes in the ratio Rz,/Rrp occurred during this period. 
With the temperature coefficients of the wire apparently remaining constant 
the increase in the resistance during heat treatment is presumably caused by 
extension or volatilization of the platinum wire (Corrucini 1951). 

Fig. 6 shows the variation of the reference resistance Rrp for a thermometer 
that has been used for several thousands of measurements at the freezing 
points of zinc, cadmium, and tin. Table II lists some resistance temperature 
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Ryp 25.48 XXX b.u. 


200 
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Fic. 6. Variation in the annealed resistance at Rrp for a thermometer used constantly in 
the range 0° to 450° C. over a period of three years. 


TABLE II 
Relative change, 
Sept. 1953 Feb. 1956 p.p.m. 
Rrp 25.48220 b.u. 25.48298 b.u. +31 
Ran 65.43966 b.u. 65.44181 b.u. +33 
Rzn—Rrp 39.95746 b.u. 39.95883 b.u. +34 
Rz2/Rrp 2.568054 2.568060 + 2.3 


relations that were determined for this thermometer in September, 1953, and 
in February, 1956. The change in the ratio Rz,/Rrp implies an increase in the 
temperature coefficients of the platinum wire, presumably due to the long term 
heating. This increase of 2 p.p.m. corresponds to an increase of about 0.002° C. 
in a calculated zinc temperature. 

Further studies on the long term stability of the temperature coefficients of 
platinum thermometers and of factors which may account for the observed 
changes are underway. 


(b) Water Triple Point Cells 

During November and December, 1954, a series of temperature measure- 
ments in several triple point cells prepared by Dr. R. S. Turgel, formerly of 
this laboratory, was designed to determine the variations in temperature be- 
tween cells and any temperature dependence on the age of the ice mantles. 
Before making measurements all the thermometers used were annealed at 
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about 450° C. until recovery from previous cold work was completed and Rrp 
had begun to rise slowly; thereafter these thermometers were handled gently 
and were kept in the temperature range 0° C. to 23° C. The values of Rrp were 
measured with 1 ma. and 2 ma. currents through the thermometer and extra- 
polated to the zero current resistance. 

Frequently it is found that the value of Rrp obtained on the day the ice 
mantle is prepared is lower by as much as 10 microhms than the average of the 
week. For this reason it is common practice here to allow a full day for the 
interface to attain temperature equilibrium before a cell is used in precision 
temperature determinations. 

Fig. 7 shows the results* of monitoring Rrp with a single thermometer on a 


Rrp 25.483XX b.0 


AV. Rrp 25.483399 b.2 
STO. DEV. $2.7 wa 


SPREAD 9 »Q 





Fic. 7. Monitoring Rrp with a single thermometer on a single water triple point cell over a 
period of five days. 


single cell for five days, the ice mantle being prepared on the first day. The 
standard deviation about the mean Rrp was 2.7 microhms with a maximum 
spread of 9 microhms. This is the smallest variation over this length of time 
that has been obtained in this laboratory, but it must be emphasized that the 
thermometer was used under conditions that normally are impractical, i.e. 
with a minimum of handling and in the absence of high temperature environ- 
ments; the standard deviation about the mean Ryp of a thermometer that is 
being frequently handled and used at a variety of temperatures is of the order 
of two or three times this value. 


Rrp 25.483XX b.Q. 


TC. 25@ T.C. 270) 4 
OS 2 
V. Ryp25.48340,0.2. 25.48340, b.Q. 
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Fic. 8. Monitoring Rrp with a single thermometer on two water triple point cells over a 
period of five days. 


*In Figs. 7, 8, and 9 and Tables III and IV the values of Rrp are given in arbitrary bridge 
ohms; i.e. the bridge reading is corrected to true zero and unity ratio but not to bridge units of 
resistance as defined above. This latter correction is unnecessary in view of the small differences 
in temperature between triple point cells. 
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Fig. 8 shows the results of monitoring two cells with the same thermometer 
for a similar period. Both the standard deviation and the spread from the 
mean Rrp for each cell are measurably larger than those shown in Fig. 7. The 
graph of the daily averages of Rrp clearly indicates a difference in temperature 
realized in the two cells. 

An intercomparison of several triple point cells was carried out along the 
following lines. Each week ice mantles were prepared on four cells and were 
allowed one day to attain stability. Two series of measurements were made 
concurrently during the next four days, one of these being single daily de- 
terminations of Ryp on each cell with the same thermometer, S156, and the 
second consisting of three determinations of Ryp with each of four combinations 
of a thermometer and a cell, different combinations being measured each day: 
the use of four thermometers tends to average out errors arising from increases 
in Ryp due to accidental cold working of the thermometers during handling. 
Tables III to VI give the results of the intercomparisons, T.C. 27 being used 




















TABLE III 


SINGLE THERMOMETER INTERCOMPARISON (1954) 
Rrp on thermometer S156 in bridge ohms 













Deviation of 
daily 
average 
Daily av. from weekly 
PC 2e TCH25 TC se T.C. 33 for group average 

















‘ 25.483406 25.483415 25.483404 25.483415 25.483410 —2 uQ 
24 409 413 416 415 413 =P 
25 407 414 412 416 412 0 
26 407 425 414 413 415 +3 


Weekly average 25.483407 25.483417 25.483412 25.483415 25.483412 
Rel. to T.C. 27 0 pQ +10 po +5 pQ +8 uo 
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N 25.483408  25.483409 25.483406 25.483414 25.483409 

Dec. 1 405 413 408 411 409 —1 
2 410 412 410 412 411 ss 
3 411 411 412 414 412 








Weekly average 25.483408 25.483411 25.483409 25.483413 25.483410 
Rel. to T.C. 27 0 nQ +3 pQ +1 wQ +5 wQ 
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25.483413 25.483404 25.483415 25.483418  25.483412 
415 394 417 421 412 —1 
416 404 419 420 415 +2 
423 388 416 423 412 


Weekly average 25.483417 25.483398 25.483417 25.483420 25.483413 
Rel. to T.C. 27 0 pQ —19 pQ 0 pQ +3 nQ 
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Cell 
Nov. 23 27 
24 25 
25 32 
26 33 


Weekly average 


Nov. 30 27 
Dec. 1 22 
2 21 
3 24 


Weekly average 





Cell 

Dec. 7 27 
8 25 

9 29 

10 31 


Weekly average 
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TABLE IV 


FOUR-THERMOMETER INTERCOMPARISON (1954) 
Rrp in bridge ohms 
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$156 Cell 
25.483406 33 
413 27 
412 25 
413 32 
25.483411 
25.483408 22 
413 21 
410 24 
414 27 
25.483411 
$163* Cell 
25.479536 31 
550 27 
540 25 
511 29 
25.479534 


$155 
25.478053 


25.478054 


25.478054 
052 


056 
052 


25.478054 


$155 


25.478046 
053 
062 
057 


25.478054 





Cell 
32 
33 


27 
25 


Cell 


31 
25 





S165 Cell 
25.481604 25 
609 32 
606 33 
611 27 
25.481608 
25.481607 24 
611 27 
607 22 
610 21 
25.481609 
S165 Cell 
25.481603 25 
594 29 
609 31 
616 27 
25.481606 


S214 


25.485514 
511 
511 
510 


25.485512 
25.485512 
505 
514 
509 


25.485510 


$214 


25.485515 
509 
483 
511 


25.485504 


*Thermometer S163 substituted for S156, which was temporarily unserviceable at the begin- 


ning of this week. 


TABLE V 
DEVIATION OF Rrp IN pQ= FOR INDIVIDUAL CELLS FROM 


THE WEEKLY GROUP AVERAGE (1954) 








S156 S155. S165 
7.0.37 -830 —4 —2 uo 
3% +3 +2 +3 
320+] +5 —4 
33. +2 oJ +1 
TT... +80 <3  =-8.8 
23 +42 0 +1 
oy ~¥ =f 
2% +43 +2 +2 
S163* S155 S165 
T.C.27 +2 -1 +30 
31 —23 -8 —12 
29 +6 +3 ~ 3 
25 +16 +8 +10 


Average 
Sum four deviation 
S214 thermometers thermometer 
—2 uw —13 pQ —3 wf 
+2 +9 +2 
-1 +1 0 
—1 +1 0 
—5 uw —12 pQ —3 uo 
+4 + 7 +2 
—1 — 6 —2 
+2 +9 +2 
$214 
+ 72 +11 uo + 32 
—21 —64 -16 
+ 5 +11 + 3 
+11 +45 +11 


Cell Rrp 
rel. to 
4.G. ae 


0 ua 
—19 


0 
++ 3 


*Thermometer S163 substituted for S156, which was temporarily unserviceable at the begin- 
ning of this week. 
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TABLE VI 


DEVIATION OF Rrp IN p& OF THE DAILY GROUP AVERAGE 
FROM THE WEEKLY GROUP AVERAGE (1954) 


Deviation of 
daily average from 
Sum four weekly average 
thermometers (sum/4) 


—8 pQ —2 pQ 
-—2 0 


0 0 
+8 +2 


S155 


a + 22 —8 uh — 3nQ +11 »pQ +2 wQ 0 pQ 
8 +16 —1 —12 + 5 +8 +2 
9 +6 +8 + 3 —21 —4 —1 
10 —23 +3 +10 +7 —3 —l 


*Thermometer $163 substituted for S156, which was temporarily unserviceable at the begin- 
ning of this week. 


as a reference cell. An additional intercomparison between Rrp values of 
National Research Council triple cell 25 and a cell loaned to us by the J. and J. 
Instrument Co., Silver Spring, Md., was carried out in March, 1955, with 
thermometer S156 and the results are presented in Fig. 9. 


N.R.C. (0) 
25.48345,, b.2 25.48345,,b.2 


Ryp 25.483xx b. 2 


2 Rip 


Fic. 9. Intercomparison of Rrp values for a J. and J. cell and an N.R.C. cell. 


Table VII(a) summarizes the intercomparison of National Research Council 
triple point cells by relating all cells to T.C. 27 in terms of the measured 
average differences in Rrp in microhms. The weighted mean differences were 
calculated because the four-thermometer measurements include twice as many 
determinations of Rp on each cell as the single-thermometer ones. The excel- 
lent agreement between the results of the two series of measurements justifies 
the conclusion that there are differences in Rrp between cells of less than 
10 microhms that are real and measurable. 

For a 25.5 ohm standard resistance thermometer at the triple point of water, 
a change in resistance of 10 microhms is equivalent to about 0.0001° C. so that 
the average temperatures of eight of the National Research Council cells were 
well inside 0.0001° C. The remaining cell, yielding a temperature about 
0.0002° C. lower, was discarded as being substandard. 
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TABLE VII 
SumMMARY (1954) 
(a) AR¢p of cells relative to T.C. 27 


Single Four Weighted 
thermometer thermometer mean 
comparison comparison difference 


5 
5 
5 


+444 


+444 
AkrOOOTR— 


++1 0 ++4++ 
wCnrDsoonw-~ 
++ | 
wWwSonnar 
+41 


(b) ARrp with age of ice mantle 


Deviation of daily average from weekly average 
Week 2nd day 3rd day 4th day 5th day 


Single thermometer 1 —2 uQ2 +1 uw 0 pQ +3 uQ 
2 —1 —1 +1 +2 
3 —1 —} +2 =) 


Average deviation —1 0 +1 +1 


Four thermometers —2 0 0 +2 
-1 -1 +1 0 
0 +2 -1 -1 


Average deviation -1 0 0 


In order to obtain information on the long term stability of the cells a second 
intercomparison was made during April and May, 1956. Table VIII shows the 
results of the two intercomparisons. No measurable change in the relative 
differences was observed for six cells, two cells indicated changes only slightly 
greater than the uncertainty of measurement, which is considered to be of the 


TABLE VIII 


Rrp OF CELLS RELATIVE TO T.C. 27 
The limiting uncertainty in a difference measurement 
is of the order of 3 wQ 


Cell Nov. Dec./54 Apr. May/56 Change in Rrp 
0 pQ 0 pQ 


thet 
5 


Ark OOoukr- 


— 

Wo © Pe DO 
1 

Lh mla&lool | 


— 


++ 141444 
~f 


++! 
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order of three microhms, and the remaining cell showed a temperature drop 
of 0.00015° C. This decrease in triple point temperature was caused by air 
entering the cell, possibly through a very fine crack in the pyrex envelope. 

In these experiments there was no indication of a correlation between the 
differences in equilibrium temperatures among the eight National Research 
Council cells and the hydrostatic head of water on the interface. As an addi- 
tional check on this matter temperatures against depths of immersion of a 
thermometer were measured in a cell with an ice mantle 14 in. long. Although 
some evidence was obtained suggesting a trend of decreasing temperatures at 
greater depths, of the same order of magnitude as the standard pressure effect 
on the freezing temperature of water, the reproducibility was poor; it is possible 
that convection currents in the thermometer well result in variations from the 
expected temperature gradient. 

Table VII(b) shows an analysis of the 1954 data with a view to detecting 
any changes of cell temperature with the age of the ice mantles from the second 
to the fifth day. The variations found are well inside the uncertainty of the 
measurements. 


CONCLUSION 


These experiments suggest that the practical limit in precision temperature 
determinations is of the order of 10-*°C, (10-*° C. = 10 uQ for 25.5 ohm 
platinum thermometers) using present-day equipment. The uncertainty due 
to the detection and calibration of resistances of the platinum thermometer is 
somewhat less than 10 microhms and variations in temperature among triple 
point cells appears to be of the order of 10-4° C. Normal handling of well- 
annealed platinum resistance thermometers may introduce errors of this order 
of magnitude due to accidental straining of the wire during the course of a 
measurement. 

These similar limitations placed on temperature measurement by the bridge, 
triple point cells, and the platinum thermometers, resulting from a well- 
balanced experimental technique, combine to form a formidable barrier to 
any substantial increase in precision. 
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PROPERTIES OF GALLIUM INDIUM ANTIMONIDE! 


J. S. BLAKEMORE? 


ABSTRACT 


Electrical conductivity and Hall effect are measured for p-type specimens of 
polycrystalline GaSb.InSb with impurity concentration about 10!7 cm. 
Analysis of these data suggests that un/up = 11, and that the intrinsic gap 
varies linearly with temperature from 0.265 ev. at 0° K. Measurement of the 
photoconductive limit at various temperatures shows that the gap widens on 
heating, though the electrical data seem difficult to reconcile with the large 
gradient of +1.1X10-* ev./° C. indicated by the optical data. 


1. INTRODUCTION 


The study of semiconducting intermetallic compounds is a matter which 
has received appreciable attention only in the last few years, yet already a 
considerable amount of systematic data has been brought to light (e.g. Welker 
1954; Frederikse and Blunt 1955). Much of this work has been concerned with 
compounds formed between elements of the third and fifth groups of the 
periodic table, such as InSb. 

One aspect of semiconductor investigation concerns the properties of alloys 
in various proportions between two related. semiconductors. Thus the proper- 
ties of germanium-silicon alloys (Levitas et a/. 1954) are found to vary con- 
tinuously according to the proportion of silicon in the alloy, whilst Nussbaum 
(1954) has found a similar gradation of properties in Te-Se alloys. When we 
come to alloys between two semiconducting compounds, Busch and Winkler 
(1953) have found that the intrinsic gap in alloys of the type Mg2Ge,Snii_y) 
varies for different values of y between the gap values for Mg:Ge and Mg:Sn. 
Welker (1955) finds that the gap width in ternary alloys between GaP and 
GaAs also varies regularly between the values for the binary phases which 
form the limits of composition.* 

In this paper, results are presented on one composition of the ternary phase 
between InSb and GaSb, that represented by the equimolecular composition 
GaSb.InSb. Goryunova and Fedorova (1955) have found that a wide range of 
solid solution exists between these binaries. These authors found that X-ray 
powder pictures showed bands rather than sharp lines, indicative of a range of 
interatomic spacings. This suggests that in GaSb.InSb we have a disordered 
alloy rather than an ordered compound with the chalcopyrite structure. This 
hypothesis is supported by the nature of the solidus curve for this ternary 
system (Koster and Thoma 1955). 


1Manuscript received September 4, 1956. 
Contribution from the Physics Department, University of British Columbia, Vancouver, 


2Now at Minneapolis-Honeywell Research Center, Hopkins, Minn., U.S.A. 

*These ternary alloys permit a continuous range of composition and are quite distinct from 
the ternaries discussed by Goodman and Douglas (1954), of the type ABX2 where A comes 
from the first group of the periodic table, B from the third, and X from the sixth. The latter 
compounds will achieve the electronic balance for semiconducting behavior only when A and B 
are present in equal atomic proportions. 


91 





92 CANADIAN JOURNAL OF PHYSICS. VOL. 35, 1957 


By analogy with the results of Welker (1955) quoted above we might expect 
GaSb.InSb to exhibit an intrinsic gap somewhere between the limits of 0.23 ev. 
(InSb) and 0.77 ev. (GaSb) at absolute zero with a negative temperature co- 
efficient of some —4X10-‘ ev./° C. (e.g. Roberts and Quarrington 1955). 

It is of interest to know how mobilities and effective masses in GaSb.InSb 
will compare with those in the parent binaries. The valence band properties in 
GaSb and InSb are not very different, the mobility at room temperature 
being about 700 cm.?/volt. second for pure samples of either, whilst the 
effective mass ratio (m,/mpo) differs only by a factor of two between the com- 
pounds. The conduction bands in the two binaries, however, show marked 
differences. That in GaSb is “‘normal’’, with an effective mass ratio of some 0.2, 
whilst the corresponding electron mobility is moderate. But electron mobilities 
in InSb are extremely large, corresponding with a very small electron mass ratio 
of some 0.015 (Frederikse and Blunt 1955). The consequent small density of 
states renders degeneracy likely for n-type samples (Burstein 1954) with 
consequent complexity in interpretation of measurements. It is clearly of 
interest to know how the behavior of the conduction band in GaSb.InSb com- 
pares with the differing attributes of this band in GaSb and InSb. 


2. EXPERIMENTAL PROCEDURES AND RESULTS 


(1) Material Used 

The specimens used all came from a single ingot kindly made available by 
Standard Telecommunication Laboratories Ltd., of Enfield, England. The 
electrical properties showed the material to be p-type, with an effective 
acceptor concentration slightly over 10'7 cm.—*. This crystal was produced by 
the horizontal zone melting method, and was polycrystalline with an average 
grain size of about 0.5 mm. There were a number of cracks near the outer 
surface of the ingot, and it was mechanically rather fragile, presumably owing 
to the polycrystalline character. Cleaved or polished surfaces were fairly 
bright with a slight bluish tint, and the coloration was enhanced on etching 


with agents such as CP4. 


(2) Electrical Properties 

Specimens were cut as rectangular slabs for measurement of electrical con- 
ductivity and Hall coefficient by standard d-c. techniques. Contacts were 
applied by soldering with lead-tin solder or indium, and at high temperatures 
pressure contacts were used. Observations of conductivity were made over the 
temperature range 77° K. to 750° K., whilst the Hall effect measurements 
were made up to 600° K. The Hall coefficient was constant in the extrinsic 
range, indicating that the extrinsic activation energy is very small. The elec- 
trical conductivity varied slowly below room temperature, reflecting a small 
temperature dependence of hole mobility. Results for one specimen in the 
temperature region of maximum interest are shown in Fig. 1. 


(3) Optical Absorption 

Several slices of material were reduced in thickness by grinding and the 
surfaces then polished. Owing to the fragility of thin sections, only one sur- 
vived, which had been reduced only as far as 0.5 mm. The optical transmission 





BLAKEMORE: GALLIUM INDIUM ANTIMONIDE 


a 
° 


HALL COEFFICIENT Ry 
nn 
° 


o 


CONDUCTIVITY 


3-0 38 
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Fic. 1. Electrical conductivity ¢ (ohm™ cm.) and Hall coefficient Rg (cm.~*/coulomb) 
as a function of reciprocal temperature for specimen 108, in which (Na— Na) = 1.7X10!7cm.-3. 


was measured at 295° K. using conventional techniques with a Perkin Elmer 
model 83 monochromator and vacuum thermocouple. Assuming a refractive 
index of approximately four, the reflection coefficient should be 0.36, and this 
figure is used to convert values of transmission factor into the corresponding 
ones of absorption coefficient, shown in Fig. 2 as a function of wavelength. It 
will be noted that the absorption edge is very diffuse, and that the minimum 
absorption coefficient is still quite high. 


250 


ABSORPTION COEFFICIENT 
8 
o 


a 
°o 


4 8 12 16 
WAVELENGTH (microns) 


Fic. 2. Optical absorption coefficient a (cm.~!) as a function of wavelength for a slice 
0.5 mm. thick. 
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(4) Photoconductivity 

The experimental arrangement used for photoconductive observations has 
been described previously (Blakemore 1956). Specimens were prepared with 
four contacts. In view of the very low resistance of these cells, a high ratio 
step-up transformer preceded the tuned amplifier system. 

The photoconductive signal was the same for a light beam interrupted at 
90 c.p.s. and at 870 c.p.s., indicating an electron lifetime of less than 10~* 
seconds. The spectral dependence of photoconduction, expressed as the frac- 
tional increase of conductivity per incident photon, is shown in Fig. 3 for one 


20 
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Fic. 3. Intrinsic photoconductive response plotted on an equal-photon basis, as a function 
of photon energy for three temperatures (specimen 125, in which (Na— Na) = 2. 3X10! cm. =e}, 


specimen measured at three temperatures. This is fully representative of the 
data collected on other cells. The units for the ordinate of this figure are the 
same for the three curves: unity on this is equivalent to an electron lifetime of 
10-* seconds for unit quantum efficiency. 

This figure reveals at once a striking and unexpected fact: that the intrinsic 
response extends further into the infrared at low temperatures, i.e. the intrinsic 
energy gap increases with temperature. 

The position of the ‘“‘edge”’ is quite sharp at 0.36 ev. for 90° K. with all the 
cells studied. At the higher temperatures the edge is more diffuse. By con- 
sidering the points of intersection of lines drawn from the high energy and low 
energy portions of these curves, optical gaps of about 0.51 ev. at 218° K. and 
0.59 ev. at 295° K. are deduced. 

3. DISCUSSION 


The interpretation of electrical conductivity and Hall coefficient to obtain 
values of the electron and hole concentrations m and p, and the corresponding 
mobilities 4, and uw, = un»/b, has been described by Pearson and Bardeen 
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(1949) and by Beckenridge et al. (1954). It need only be briefly indicated here. 
The equations at our disposal are 








_ 32) pont 
(1) Ru = 8e (p-+nb) ’ 
(2) o = eu,(p+nb) , 






and at all temperatures 
(3) (p—n) = (Na—Na), 


where N, and Ng are the densities of acceptor and donor impurity centers, 
assumed completely ionized. In a stable material (V,— 4) will be unchanged 
by heating, but in GaSb.InSb we find that samples heated above 500° K. show 
signs of a decrease in this quantity which is only partly reversed on cooling. 
In consequence, results computed above 500° K. are not taken into account in 
considering the temperature variation of the intrinsic carrier concentration. 

From the limiting value of Ry at low temperatures it emerges that for 
specimen 108 (as illustrated in Fig. 1) the effective acceptor concentration 
(Na—Na) = 1.7X10'7 cm.-*. Comparing this with the negative maximum of 
Ry, which is 













3x (6—1)? 
(4) Raax = — 3, X 45(N— No)’ 
it is found that in this material the mobility ratio b = u,»/u, = 11. These two 
quantities enable us to find the values of 2 and p corresponding with the value 
of the Hall coefficient at any temperature. 

Simultaneous knowledge of the electrical conductivity o then permits 
evaluation of the hole mobility u,, as shown in Fig. 4. The steeply sloping 
portion for high temperatures, where the form of variation is approximately 
u, « T-*7, results from lattice scattering, whilst the more nearly constant 
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Fic. 4. Free hole mobility u, (cm.?/volt.second) as a function of temperature for specimen 
108. 
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mobility at lower temperatures is presumably the result of impurity and dis- 


location scattering. 
Let us now return to consideration of carrier concentrations. For specimen 


108, values of m and p are deduced from the Hall coefficient for temperatures 
above about 280° K. (At lower temperatures Ry lies too close to the exhaustion 
value to permit accurate evaluation of n.) Above this temperature one can thus 
compute the intrinsic concentration »; = (np), which is also given by 
(5) n, = 2(2emokT/h*)?” M3” exp(—E/2kT) 

4.83 X10!5 (MT)*” exp(—E/2kT) cm.- 


where the composite mass factor M = (m,m,/m,*)}, and E is the intrinsic gap 


width. 
The corresponding data for m, on this specimen are shown in Fig. 5, where 


n,/T?” is given as a function of 1000/7. Over the temperature range up to 


2°5 3-0 35 
“1 
1000/T (°K ) 


Fic. 5. Variation of n;/T*? against 1000/T for specimen 108. For n; in units of cm.~?. 


500° K. the behavior is represented quite closely by a straight line on this 
semilogarithmic plot, which suggests that the variation of intrinsic gap with 
temperature is essentially linear. The data of Fig. 5 can be expressed as 


(6) n;, = 2.5X10" T?” exp(—0.265/2kT) cm.-3 
and if we assume that 
(7) E = E,+8T, 


then the identification follows that Ey = 0.265 ev., whilst 
(8) exp(@/2k) = 19.3 M3”, 


so that a knowledge of the effective mass ratio is contingent on that of the 
temperature dependence of the intrinsic gap. The optical transmission and 
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photoconductive studies reported in Section 2 were carried out with this in 
mind. Whilst the optical absorption data of Fig. 2 do not provide an absorption 
edge sufficiently definite to be of much assistance, the photoconductive values 
for the intrinsic edge are reasonably well defined. The critical photon energies 
for the three temperatures do in fact fit quite well with a variation 


(9) E = 0.265+1.1X10-T ev. 


and the occurrence of 0.265 ev. as the extrapolation to 0° K. is indeed remark- 
able. However, the positive temperature coefficient of the intrinsic gap seems 
inordinately large by comparison with the values found for other semicon- 
ductors. Moreover, if this value for 8 is substituted in Equation 8, the result is 
obtained that M = 10, which is larger by at least an order of magnitude than 
would be expected from the values for other intermetallic compounds. More- 
over, such a large effective mass ratio is inconsistent with the fairly high 
mobilities noted both for electrons and for holes in this material. 

In these circumstances no definite conclusions can be drawn concerning the 
place of this ternary alloy in the scheme of semiconducting compounds until 
more information is available, preferably on both n- and p-type crystals of 
widely varying impurity content. 
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SELECTIVE REFLECTION FROM MERCURY VAPOR AT 
HIGH PRESSURES'! 


J. A. GaLt? anp H. L. WELsuH 


ABSTRACT 


Selective reflection from mercury vapor in the region of the 2537 A resonance 
line was investigated at pressures up to 340 atm. using reflection cells of special 
design. The results were interpreted on the basis of the classical theory of re- 
flection from an absorbing medium. By fitting calculated curves to the experi- 
mental reflection contours, values of the oscillator strength, f, and the damping 
constant, y, were determined. The f values so obtained are density-dependent 
and at high pressures are approximately equal to one half the value for the 
free atom. As predicted by theory, the damping constant varies directly as the 
density of the atoms in the vapor. This result contradicts the earlier work of 
Welsh, Kastner, and Lauriston (1950) in which it was concluded that y varies 
as the square root of the density. A subsidiary reflection maximum was observed 
at 2540 A; it is attributed to Hg: molecules which occur in relatively large 
concentrations at high densities. Some preliminary observations on the selective 
reflection at the 1850 A resonance line were made up to 4.4 atm. 


INTRODUCTION 

Selective reflection from a metallic vapor in the neighborhood of its 
resonance lines, discovered by Wood (1909) in mercury vapor at 2537 A, has 
been studied more recently by Welsh, Kastner, and Lauriston (1950) in mer- 
cury and cadmium vapors and by Lauriston and Welsh (1951) in the vapors 
of the alkali metals. In these later investigations light of continuous spectral 
distribution was reflected from the inner surface of a window in contact with 
the vapor and the reflectivity of the interface was determined as a function 
of frequency. Vapor pressures in the range 0.1 to 15 atm. were used for mer- 
cury and cadmium, and 0.02 to 1 atm. for the alkali metals. Selective reflec- 
tion from mercury vapor has been studied at low vapor pressures by Cojan 
(1954), employing ingenious methods to resolve the individual hyperfine 
components of the 2537 A line. 

The classical theory of reflection from an absorbing medium was used 
by Welsh et al. to fit the experimental reflection contours with calculated 
curves obtained by varying the values of the oscillator strength, f, and the 
damping constant, y, in the anomalous dispersion equations. In the medium 
pressure range the experimental contours could be reasonably well reproduced 
by using a constant value of f and values of y which varied as the square root 
of the number of atoms, Na, per cubic centimeter. The f value so determined 
agreed, in general, with the value determined by other methods. At pressures 
above 5 atm., however, it did not appear possible to obtain calculated con- 
tours in agreement with experiment by any variation of f and y. The density 
dependence of the damping constant, y « +/N,, deduced in these selective 
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reflection experiments differs fundamentally from the dependence, y « N,, 
which holds quite generally for emission and absorption lines. This apparent 
difference between the broadening processes involved in selective reflection 
and in absorption was difficult to understand. 

The purpose of the present work was to investigate the apparently anoma- 
lous behavior of the selective reflection at high densities. The mercury line 
at 2537 A was used, since it lies in a convenient spectral region, and since 
mercury vapor pressures of several hundred atmospheres can be produced 
without undue experimental difficulties. 

Selective reflection in the neighborhood of the strong resonance line 
of mercury at 1850 A has not been previously investigated; some preliminary 
observations on this line are reported in Appendix I. 















EXPERIMENTAL PROCEDURE 


To generate vapor pressures of mercury up to, say, 400 atm. requires 
temperatures up to 1100° C. For pressures lower than 50 atm. the mercury 
can be contained in a sealed fused quartz tube with thick walls and heated 
in a simple furnace. At higher pressures it would be necessary to apply an 
external gas pressure to balance the vapor pressure; the reflection cell and 
furnace would thus have to be placed in a pressure vessel. However, since 
at high temperatures fused quartz becomes permeable to gases, the gas in 
the pressure vessel would diffuse through the window and walls of the cell 
and mix with the mercury vapor. A satisfactory design of the reflection cell 
thus presents some difficulties. In the design adopted here the walls of the 
cell were of stainless steel and the window was a long cone of fused quartz. 
A temperature difference of over 1000°C. could be maintained between 
the lower end of the cone, in contact with hot mercury vapor, and the upper 
end, which was sealed to the steel wall at room temperature. 

The essential details of the reflection cell and the pressure vessel are given 
in Fig. 1. The window of the reflection cell was a truncated cone of fused 
quartz of high ultraviolet transparency, 10 cm. long, with polished upper 
and lower surfaces. The cone fitted snugly into a thin steel shell, B, to which 
it was sealed, near the top, by a Teflon packing. A coil of platinum wire em- 
bedded in refractory cement, C, was used to heat the reflection surface at the 
lower end of the quartz cone. A long tubular stem, a continuation of the 
steel shell, connected the space just beneath the reflecting surface with the 
open reservoir of mercury, D. Thermal insulation was provided by Alundum 
grain surrounding the furnace; however, the insulation was quite inefficient 
when the Alundum was permeated by high density gas. A copper cooling ring, 
F, helped to keep the mercury in the reservoir cool by conducting heat to the 
walls of the pressure vessel, which were water-cooled. The parts were as- 
sembled, as in Fig. 1(a), in a steel tube, G, which was placed first in a vacuum 
chamber for outgassing and filling with mercury, and then in the pressure 
vessel shown in Fig. 1(d). 

The pressure vessel was designed to withstand pressures up to 1500 atm. 
A thick fused quartz window, K, was contacted to a holder by a thin film 
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Fic. 1. Apparatus for studying selective reflection from high pressure mercury vapor: 
(a) reflection cell, (b) high pressure vessel. A—fused quartz cone; B—stainless steel shell; 
C—platinum wire furnace; D—liquid mercury; E—Alundum insulating grain; F—cooling 
ring; G—steel tube; H—slots to reduce heat conduction along G; J—electrical connector for 
furnace power; J—constriction in stem to reduce heat losses; K—pressure window; L—Teflon 
seal; M—window holder; N—hardened steel thrust ring; O—water bath; P—gas inlet. 












101 





GALT AND WELSH: REFLECTION FROM MERCURY VAPOR 


of rubber cement. The free area of the window was divided into two parts 
to allow the light to enter and leave by separate paths. The stainless steel 
window holder was sealed to the body of the pressure vessel by a Teflon 
gasket, L, which was compressed by eight thrust bolts carried in a retaining 
cap screwed to the body of the vessel. A high current lead to convey 60-cycle 
power to the furnace entered the bottom of the vessel through a pressure 
seal with Teflon packing. Gas from a hand-operated compressor entered 
the vessel through 1/16 in. steel capillary tubing. Helium was used to furnish 
the balancing pressure notwithstanding its high thermal conductivity, since, 
of the gases commercially available, it had the highest ultraviolet trans- 
parency. A water bath surrounding the pressure vessel dissipated heat con- 
ducted from the furnace in order to keep the mercury reservoir cool and to 
protect the Teflon seals. To facilitate optical adjustments the pressure vessel 
was mounted on a base which permitted fine horizontal motion in two per- 
pendicular directions. 

The action of the apparatus described can now be clarified. The reflection 
cell, previously filled with mercury in the auxiliary vacuum chamber, re- 
mains full until sufficient heat is applied to create a vapor pressure equal to 
the gas pressure. The level of the liquid mercury then positions itself so that 
the pressure of the vapor just beneath the quartz cone is exactly balanced by 
the pressure of the gas, which is measured by a Bourdon gauge. 

With the apparatus described it was not possible to heat the mercury much 
higher than 1100° C. without melting the platinum heater. An attempt was 
made to increase the efficiency of the furnace by introducing mercury into 
the cell from above, thus eliminating conduction losses in the long stem 
below the furnace; this gave only a slight improvement. That most of the 
heat losses occurred by convection and conduction in the gas was indicated 
by the fact that much more power was required to maintain a given tem- 
perature at high gas pressures than at lower pressures. 

At temperatures above 1000°C. recrystallization of the fused quartz 
caused the reflection surface to deteriorate rapidly; the useful lifetime of 
the cell at high temperatures was therefore short and the surface required 
repolishing after every experiment. Another factor contributing to the diffi- 
culties of working at high temperatures was the corrosion of the steel walls 
of the cell by hot mercury. Small quantities of titanium and magnesium 
added to the mercury reduced, but did not eliminate, this difficulty. The 
steel shell had therefore to be renewed frequently. 

A hydrogen ultraviolet continuum of the high intensity necessary for the 
experiments was produced by a discharge, at 3000 v. and 1 amp., between 
large aluminum electrodes in a water-cooled tube of fused quartz. The dis- 
charge, confined to a capillary tube 6 mm. in diameter and 50 cm. long, 
was viewed end-on. To stabilize the intensity of the lamp a small portion of 
its light was directed onto a photocell, the output of which was amplified 
and used to control the lamp current by means of a saturable reactor in series 
with the primary of the lamp transformer. 

By means of a suitable arrangement of mirrors and lenses the light from 
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the hydrogen lamp was allowed to fall on the lower surface of the quartz cone 
(Fig. 1(a)) and the reflected light was directed into the spectrograph. The 
pressure window and the upper surface of the quartz cone were tilted slightly 
so that light reflected from them did not enter the spectrograph. A Hilger 
single prism instrument of the Littrow type was used, giving a reciprocal 
linear dispersion at 2537 A of 43 cm.-!/mm. In place of the usual photo- 
graphic plate a detector unit, consisting of an exit slit and a 1P28 photo- 
multiplier tube, was constructed to move along the focal plane of the spectro- 
graph. The photomultiplier tube was coupled through a preamplifier to an 
electronic pen recorder. 

In studying the selective reflection at 2537 A a spectral region extending 
500 cm.~! on either side of the resonance line was scanned. Silicon emission 
lines, present in the spectrum of the hydrogen lamp, were used for frequency 
calibration. At each vapor pressure several recorder traces of the selective 
reflection were made and averaged to give a mean observed contour. To 
obtain relative reflected intensities the mean contours were modified by using 
background traces obtained by reflection from the cold cell before the mer- 
cury was introduced. 

To establish a scale of percentage reflected intensity, it was assumed that 
the reflectivity in the wings of the pattern approaches asymptotically a 
value Ry at large distances from the resonance frequency. It was further 
assumed that the asymptotic reflectivity is given by R = (n’—n9)?/(n’+1)? 
for perpendicular incidence, where ’ is the refractive index of fused quartz 
and mp» (~1) is the hypothetical refractive index of the mercury vapor in the 
absence of anomalous dispersion at 2537 A. The value of mo as a function 
of the density of the vapor was obtained by omitting the 2537 A term in the 
dispersion equation of Wolfsohn (1933). Since the ultraviolet refractive 
index of fused quartz has been measured only over a limited range of lower 
temperatures (International Critical Tables 1929), the values for high tem- 
peratures had to be obtained by extrapolation; data for the visible, which 
cover a larger temperature range, were used to estimate the proper trend in 
the extrapolation. 

The calculation of N,, the number of free mercury atoms per cubic centi- 
meter, is complicated by the fact that in high pressure mercury vapor there 
is an appreciable concentration of diatomic molecules. The temperature of 
the vapor at the window of the reflection cell was calculated from the balancing 
gas pressure, P, by using the vapor pressure data of Bernhardt (1925), and 
the density of atoms plus molecules, N, was determined from the relation 
N = P/kT. The densities obtained in this way are probably somewhat too 
high because the vapor was slightly superheated. Since NV = N,+N,,, where 
Nm is the number of Hg» molecules per cubic centimeter, the value of N, 
can be obtained only by estimating the molecular density, V,. From the total 
absorption of the band system at 2540 A, ascribed to the Van der Waals 
Hg. molecule, Kuhn (1937) determined the absolute concentration of Hg» 
molecules at a pressure of 8 mm. Hg. In this calculation he assumed that the 
oscillator strength, f, of the transition in the molecule is the same as in the 
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free atom. He then calculated the dissociation energy, Do, of the ground state 
of the molecule from the Gibson—Heitler (1928) formula for the equilibrium 
constant using an estimated value of the moment of inertia, J, and a vibra- 
tional frequency, wo, determined from the structure of the 2341-2311 A 
band system of the molecule. A reversal of this procedure was used in this 
investigation to estimate the molecular concentrations in the high pressure 
vapor. The Gibson—Heitler equation was used in the form 


2 3/2 21/2 71/2 
(Na—Nm) mk T —_ Okt) _-Dear 










Nn 4y/ahI a 
and N,,/Na was calculated at the various temperatures used in the selective 
reflection experiments by substituting Kuhn’s values of wo, 7, and Dy in the 
equation. These calculations indicated that at the highest pressure reached, 
340 atm., about 35% of the atoms are combined in Hg: molecules. 

Measurements of the selective reflection were made in the pressure range 
41 to 340 atm. Accurate observations at lower pressures were precluded 
by the tendency of the mercury to boil violently and splash against the 
window. Table I summarizes the pressures used, the corresponding tem- 
peratures calculated from the vapor pressure curve of mercury, the total 
vapor densities expressed as N/No, where Ny is Loschmidt’s number, and 
the corresponding atomic and molecular densities, N,/No and N»,/No. The 
calculated values of Ro, the asymptotic reflectivity, are also included in the 
table. 

















TABLE I 


‘TEMPERATURE, DENSITIES, AND BACKGROUND REFLECTIVITIES FOR VARIOUS MERCURY VAPOR 
PRESSURES 

























Total Atomic Molecular 
Pressure, Temp., density, density, density, Ro, 
atm. Mm N/No Na/No Nm/ No % 














41 672 11.8 11.1 0.7 3.86 
82 800 20.8 18.9 1.9 3.65 
122 880 29.0 25.7 3.3 3.43 
184 944 41.4 35.3 6.1 3.13 
204 964 45.0 38.0 7.0 3.04 
272 1016 57.7 47.6 10.1 2.75 
340 1062 69.7 55.2 14.5 2.50 












RESULTS AND DISCUSSION 


Graphs of reflectivity as a function of wave number for several mercury 
vapor pressures between 41 and 340 atm. are shown in Figs. 2 and 3. At low 
pressures the curves follow the characteristic selective reflection pattern, which 
exhibits a reduced reflectivity on the low frequency side of the resonance 
line and a greatly increased reflectivity on the high frequency side. The 
main peak, which is high and narrow at low pressures, reaches a maximum 
intensity at about 80 atm. and then decreases somewhat irregularly in in- 
tensity with further increase in pressure; it also broadens and its frequency 
decreases slightly with increasing pressure. In the lower pressure contours a 
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Fic. 2. Selective reflection from mercury vapor at medium pressures. 
Fic. 3. Selective reflection from mercury vapor at high pressures. 


weak subsidiary maximum is apparent on the low frequency side of the main 
peak; this increases in intensity as the pressure rises, gradually filling in the 
minimum of the reflection contour on the low frequency side of the resonance 
line. 

The subsidiary maximum has a wavelength of 2540 A and thus coincides 
with the maximum of the system of discrete bands which has been observed 
in absorption by Wood and Voss (1928), Kuhn and Freudenberg (1932), and 
Kuhn (1937), and in emission and absorption by Mrozowski (1949). This 
band system has been attributed to the transition, 0{-"Z4, in the Van der 
Waais molecule Hge. The subsidiary maximum in the reflection pattern can 
therefore be ascribed with reasonable certainty to Hg2 molecules. The presence 
of the molecular reflection makes the quantitative treatment of the atomic 
reflection more difficult, since the low frequency portion of the atomic reflec- 
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tion is strongly overlapped by the molecular reflection. In addition, the 
calculated number of free atoms, Ng, which is essential for the interpretation, 
can have considerable error because of the uncertainty in estimating the 
concentration of molecules. 

The normal approach to the interpretation of experimental results on 
selective reflection is through the classical anomalous dispersion equations, 
which give the refractive index, m, and the absorption index, k, as a function 
of frequency, v: 





tk, Bae vo—V__ 
(1) ort <* TS (vo—v)’+(y/40)°’ 
Nofe* (7/47) 





2 2nk = Jaye 

(2) 2armvo (vo—v) +(y/40) 

In these equations mp is the hypothetical refractive index in the absence of 
anomalous dispersion, vo the resonance frequency of the virtual oscillator in 
sec.—!, and f the oscillator strength. The reflectivity, R, for normal incidence 
is given by the Fresnel relation 


(3) R 





where n’ is the refractive index of the bounding medium. 

Several attempts were made to fit the experimental contours with reflec- 
tion curves calculated from equations (1), (2), and (3) with different basic 
assumptions. The extensive numerical work involved was carried out with 
the help of the electronic digital computer of the Computation Centre, Uni- 
versity of Toronto. The most obvious method of approach is to use the cal- 
culated values of NV, and a value of f equal to that of the free atom and to 
vary the damping constant, y, until the experimental contour is well repro- 
duced. However, no agreement between the calculated and the experimental 
curves could be obtained by this method; if the peak heights of the two 
curves were made to agree, the calculated contour was always much broader 
than the experimental contour. It became apparent that a satisfactory re- 
production of the experimental curves could be obtained only by using 
Nf in addition to y as a parameter in the curve fitting. Examples of the 
agreement obtained at various vapor pressures are shown in Fig. 4 and the 
values of N,f and y are listed in Table II. It will be noted that the fitting had 


TABLE II 
OPTICAL CONSTANTS DERIVED FROM THE REFLECTION CONTOURS OF Hg 2537 


Frequency 
Pressure, Nof, vs f shift, vo—v0', 
atm. cm.-3 cm.~} 








41 4.0X10" 1.3X10" 0.013 2 
82 8.3 2.5 0.016 6 
122 11.0 3.6 0.016 10 
184 13.0 4.8 0.014 13 
204 14.0 4.0 0.014 15 
272 17.0 5.4 0.013 23 
340 20.0 7.3 33 
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Fic. 4. Comparison of calculated selective reflection curves (dashed) with experimental 
curves. 


to be confined, in general, to the maximum and high frequency wing of the 
pattern because of the molecular reflection in the low frequency wing. No 
significance should be attached to the smaller discrepancies in curve fitting 
since these can be attributed to errors in estimating the background intensity. 

In Fig. 5 the values of y obtained from the curve-fitting procedure are 
plotted against the relative density of atoms, N,/No. The experimental points 
for the four lowest pressures were determined by averaging curves from 
three separate experiments while the points for the three highest pressures 
were obtained from a single experiment. The points define fairly accurately 
a straight line through the origin, indicating that the damping constant is 
proportional to the density. From the slope of the straight line in Fig. 5 a 
damping coefficient y/N, = 4.8X10- sec.—! cm.’ is obtained. There seem to 
be no data available on the broadening of the 2537 A line in emission or 
absorption in pure mercury vapor which can be used to calculate a broadening 
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Fic. 5. Variation of the damping constant, y, with density. In the insert the values of 
y determined by Welsh et al. (1950) are shown on an expanded scale. 


coefficient for comparison. However, Weisskopf (1933) has derived the 
theoretical formula, Avi, = (e?f/2rmyo)Na, for the half width of a spectral 
line broadened by collisions between similar atoms. Since Avy, = y/2z, the 
formula gives y/N, = e?f/mvo for the damping coefficient. If the f value is 
taken as 0.0267, the theoretical damping coefficient is 5.7X10-* sec.—! cm.’, 
in fair agreement with the value obtained from the selective reflection data. 

The Weisskopf formula was derived on the basis of an impact theory 
involving the dipole interaction of similar atoms; the statistical theory of 
Margenau and Watson (1936) yields essentially the same formula. In these 
theories the damping coefficient, y/Ng, is independent of the temperature. 
The classical collision theory of Lorentz, on the other hand, leads to a tem- 
perature dependence of the damping coefficient as given by the formula 
v/Nq = 8p?\/(xkT/u), where uw is the mass of the atom and p the “optical 
collision diameter’’. Thus, in the Lorentz theory y/N,./T, rather than y/N,, 
is constant. Since the selective reflection data have considerable experimental 
error and are confined to the interval 945° K. to 1335° K., it is not possible 
to say for certain which of the theories is in better accord with experiment; 
the evidence is, however, slightly in favor of a temperature-independent 
broadening coefficient. It is of interest to note that the collision diameter 
calculated from the Lorentz formula is 12.2 A. 

In the curve-fitting procedure outlined above, the product N,f was treated 
as a parameter rather than as a known quantity. If the values of N,f thus 
determined are divided by the corresponding calculated values of Ng, the f 
numbers listed in Table II are obtained. The results show that for all the 
pressures used in the experiments the oscillator strength is of the order of 
one half of the value, 0.0267, for the free atom. If the value for 41 atm., which 
appears to be in error, is disregarded, the f value shows a systematic decrease 
with increasing pressure; it appears to decrease from the ‘‘free atom” value 
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fairly rapidly in the pressure range 0-150 atm. and then more slowly. To 
establish the reality of this trend it would be necessary to make a careful 
study of the selective reflection in the pressure range up to 50 atm.; with 
the experimental arrangement used in the present investigation this study 
was unfortunately not possible. 

It cannot be affirmed that a dependence of the oscillator strength on den- 
sity has been established with certainty in these experiments. Nevertheless, 
a careful consideration of the sources of error indicates that the conclusion 
is at least qualitatively correct. An overestimation of N, because of super- 
heating of the vapor at the cell window or because of an underestimation of 
N,» could lead to a lowering of the f value, but not to the extent observed. 
Overlapping of the atomic reflection by the molecular reflection would have 
the effect of raising rather than lowering the f value. The calculation of 
Nm rests on Kuhn’s assumption that the f value in the Hg molecule is the 
same as for the atomic transition. If, on the other hand, the molecular f 
value is lower than that of the free atom, as is indicated by the density depen- 
dence found in this investigation, the f value determined from the selective 
reflection would be still further lowered. 

There appears to be no theoretical evidence for or against a pressure de- 
pendence of the oscillator strength. A theoretical treatment of the variation 
of f at high densities has been given by Méglich and Rompe (1942), but the 
results are not stated in a form which can be applied to the present case. 
There is no reliable experimental evidence that atomic transition probabilities 
are changed by foreign gases at high pressures; this question has been dis- 
cussed in some detail by Kuhn (1937). Also, it will be shown in a later publica- 
tion that very high pressures of foreign gases do not alter the oscillator 
strength determined from selective reflection (Galt and Welsh 1957). However, 
it can be argued that like atoms at high densities might exhibit changes in the 
f values of electronic transitions if coupling effects are present. Because of 
the high opacity of mercury vapor in the region of its resonance line, it 
would be very difficult to establish by absorption measurements that the 
f value decreases. It is, in fact, one of the advantages of the reflection method 
that optical constants can be determined in many cases which would be 
inaccessible to observation by the more conventional transmission methods. 

In fitting a calculated reflection curve to the experimental curve at a given 
pressure, it was found that the resonance frequency v9’ of the calculated 
curve did not coincide with the frequency vo of the Hg 2537 line. The fre- 
quencies vp and vo’ are marked on the graphs of Fig. 4 and the frequency 
shifts vo’ — vo, tabulated in Table II, are plotted as a function of Naf in 
Fig. 6. A shift of the resonance frequency caused by coupling of atoms has 
been predicted by Weisskopf (1933); the coupling shift, which is proportional 
to N.f (Born 1933), is plotted as a dashed line in Fig. 6. At low densities the 
experimental curve approaches the Weisskopf curve asymptotically but at 
higher densities there is a large deviation. The frequency shifts are well 
represented by a relation of the form vo’—v) = aN,+6N,2, where the linear 
term represents the Weisskopf coupling shifts. However, there seems to be 
no theoretical interpretation of the term in N,?. 
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Fic. 6. Shift of the resonance frequency, vo—vo’, as a function of Naf. 


An attempt was made to calculate the effect of the molecular reflection 
with dispersion equations involving two separate resonance terms, one for 
the 2537 A line and one for the 2540 A line. At pressures up to 82 atm. this 
method could be made to account for the shape of the molecular reflection 
with fair accuracy but at higher pressures it was totally inadequate. It appears 
that the molecular band of frequencies can not be accurately represented 
by a single frequency in the dispersion formula. 

The conclusions reached in the high pressure experiments differ radically 
from those of Welsh, Kastner, and Lauriston (1950) and of Lauriston and 
Welsh (1951) at lower pressures; a critical reconsideration of these earlier 
experiments was therefore undertaken. In particular, the validity of the 
relation y « +/N,, which was believed to be established at low pressures, 
now appears questionable. The values of y determined by Welsh, Kastner, 
and Lauriston in the pressure range 1/2 to 8 atm. are plotted on an expanded 
scale in Fig. 5, which also shows the straight line determined in the present 
investigation. The highest pressure points are in reasonable agreement with 
the straight line relation but the low pressure points all lie considerably above 
the line. It is believed that the discrepancy can be attributed to two in- 
dependent sources of error, each of which would tend to increase the values 
of y assigned at the lower pressures without greatly affecting the higher 
pressure values. The first source of error concerns the slit width corrections 
which were applied to the calculated curves before they were compared 
with the experimental curves. The correction which was applied assumed a 
spectral slit width of 1.0 cm.~', a value determined from the geometry of the 
spectrograph. It has since been established (St. John 1952) that aberrations 
in the optics of this spectrograph impose a lower limit of 1.4 cm.~' on the 
spectral slit width which can be realized in practice. The second source of 
error resulted from the assumption that the resonance line could be repre- 
sented by a single frequency, vo. However, the separation of the extreme 
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hyperfine components of the 2537 A line is about 0.7 cm.~! (Schiiler and 
Keyston 1931), which is an appreciable fraction of the width of the reflection 
pattern. In fact, the experimental curve of Cojan (1954) obtained with 
resolving power high enough to show the individual hyperfine components 
is only one third as wide as the curve of Welsh et a/. taken at approximately 
the same pressure (1/2 atm.). To determine whether these two sources of 
error could adequately explain the apparent deviation from the y « N law at 
low pressures, selective reflection curves were calculated using the anomalous 
dispersion equations involving terms for five separate resonance frequencies. 
The numerical value of NV,f was divided among the five individual hyperfine 
components in proportion to the relative intensities of lines as determined 
by Schiiler and Keyston. The curves so calculated were modified for a slit 
width of 1.4 cm.~' and compared with the experimental curve. The best 
agreement of the calculated curve with the experimental curve was obtained 
with a damping constant smaller than that originally assigned and in good 
agreement with the y « N, relation. It is therefore concluded that the measure- 
ments of Welsh, Kastner, and Lauriston are not in disagreement with the 
present results. A re-examination of the selective reflection curves in the 
neighborhood of the 2288 A line of cadmium, obtained by Welsh, Kastner, 
and Lauriston, also indicates the validity of the y « WN, relation; the results 
of these calculations are given in Appendix II. 


APPENDIX I. SELECTIVE REFLECTION FROM MERCURY VAPOR 
AT 1850 A 


A preliminary investigation of selective reflection in the neighborhood of 
the 1850 A resonance line of mercury was made for pressures up to 4.4 atm. 
The mercury was contained in a sealed fused quartz reflection cell and the 
vapor pressure estimated from the temperature of a side tube containing 
liquid mercury. Light from a low voltage hydrogen arc was reflected from the 


200 ————_—— — 0—— 


a 
-—— 


S 


y¥ (Units of 10'? sec) 


% Reflected Intensity 





(50 300 0.4 1.2 


0.8 
av cm" Na/No 


Fic. 7. Selective reflection from mercury vapor at Hg 1850: (a) comparison of calculated 
contour (dashed ) with the experimental curve for 4.4 atm., (b) variation of the damping constant 
with density. 
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inner surface of the cell window into a fluorite spectrograph, and the entire 
optical path was evacuated to eliminate atmospheric absorption. Oil sensitized 
photographic plates, Kodak 103a0 U.V., were calibrated for intensity measure- 
ments by a series of six exposures of equal duration taken through wire 
screens of various mesh sizes. 

The selective reflection contour for a pressure of 4.4 atm. is shown in 
Fig. 7(a). The large extent of the region of selective reflection even at this 
comparatively low pressure is a consequence of the high oscillator strength 
of the 1850 A resonance line. Theoretical curves were calculated to fit the 
experimental curves for pressures of 4.4, 2.2, and 1 atm. In each case the 
value of f used was 1.16, as determined by Wolfsohn (1933) from refractive 
index measurements. The very satisfactory agreement between the calculated 
and experimental curves is illustrated in Fig. 7(a). The values of y deter- 
mined from the curve fitting are plotted against the density, N,/No, in Fig. 
7(b). The relation y « N, is adequately confirmed. 





APPENDIX II. SELECTIVE REFLECTION FROM THE VAPORS OF 
CADMIUM AND THE ALKALI METALS 

Selective reflection in the neighborhood of the cadmium resonance line 
at 2288 A was investigated by Welsh, Kastner, and Lauriston (1950) for 
pressures up to 14.7 atm. They stated that their results confirmed the y « WN, 
relation found for the Hg 2537 line, but noted a marked discrepancy between 
their calculated and experimental contours at higher vapor pressures. With 
the improved computation facilities now available a new analysis of the data 
on the cadmium selective reflection was undertaken. 

The experimental contours, obtained by Welsh et a/. for cadmium vapor 
pressures of 5.4 and 14.7 atm., were fitted by calculated curves. The oscillator 
strength, f, was given the value 1.38, as determined by the earlier authors 
at lower pressures, and only the value of y was varied in the curve-fitting 
procedure. In Fig. 8(a) the calculated contour for 14.7 atm. is compared with 
the experimental curve. The fit is fairly satisfactory and could be improved 
if the experimental curve were rotated slightly; this would merely indicate 
an error in estimating the original background intensity contour. The re- 
calculated values of y are plotted in Fig. 8(b), along with the values assigned 
by the earlier authors. It can be seen that a linear dependence of y on Ng 
represents the data much better than a square root dependence. It is possible 
that the values of y at low densities are spuriously high because of inadequate 
slit width correction but there is no direct evidence to show that this is so. 
Since the higher vapor densities in these measurements were rather uncertain 
it was not considered worth while to try to establish, by more refined calcula- 
tions, a possible variation of the f value with density. 

The selective reflection from the alkali metal vapors was considered by 
Lauriston and Welsh (1951) to support the relation y « +/N,, although here 
also discrepancies were observed at higher pressures. In these experiments 
the vapor densities were deduced from the experimental reflectivity curves 
and, for this reason, it does not seem worth while to attempt a new analysis. 
For both sodium and potassium vapors the reflection maxima decreased in 
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Fic. 8. Selective reflection from cadmium vapor at Cd 2288: (a) comparison of calculated 


contour (dashed) with the experimental curve for 14.7 atm., (6) variation of the damping 
constant with density (+ Welsh, Kastner, and Lauriston (1950), @ recalculated values). 


intensity at the highest vapor pressures; this behavior is in accordance with 
the y « JN, relation but not with the y « /N, relation. 

The formula Aviy,, = (e?f/2amvo)Na, derived by Weisskopf (1933) and by 
Margenau and Watson (1936) for collision broadening by like atoms, in- 
dicates that the ratio yvo/N,f should be a constant for all spectral lines. A 
verification of this conclusion by measurements of half-widths of absorption 
lines is difficult, but the selective reflection data for Hg 2537, Hg 1850, and 
Cd 2288 can be used as a qualitative check. These data are summarized in 
Table III. Since, for these lines, the ratio yvo/N.f varies by a factor of two 
at most, although the f values differ by factors up to 100, it is concluded that 


TABLE III 
VALUES OF yv0/Naf AND p FROM SELECTIVE REFLECTION DATA 


v/ Na, yvo/ l aS» 
f value* cm.’ sec.~! cm.* sec. 


7 0.0267 0.48X10-8 2.110? 
Hg 2537 (oon 0.48 4.0 
Hg 1850 1.16 2'3 3.2 
Cd 2288 1.38* 2.4 2.3 


*The f values marked with an asterisk were determined from selective reflection. 


the formula of Weisskopf is at least qualitatively correct. In the last column 
of Table III collision cross-sections, p, calculated from the Lorentz formula, 
are tabulated. The large cross-sections calculated for Hg 1850 and Cd 2288 
as compared with Hg 2537 illustrate the lack of physical reality of the Lorentz 
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collision theory as applied to broadening by similar atoms. On the other 
hand, the large broadening coefficient for lines with high oscillator strengths 
is a natural consequence of the dipole interaction theory of Weisskopf and 
of the statistical theory of Margenau and Watson. 
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INFLUENCE OF FOREIGN GASES AT HIGH PRESSURES ON 
THE SELECTIVE REFLECTION FROM MERCURY VAPOR! 


J. A. GaLt? AnD H. L. WELSH 


ABSTRACT 


The influence of the foreign gases hydrogen, helium, nitrogen, and argon on the 
selective reflection from mercury vapor at 2537 A was studied at pressures up 
to 1500 atm. The results were interpreted on the basis of the classical theory 
of reflection from an absorbing medium. The damping constant was found to 
vary linearly with foreign gas density as predicted by collision broadening 
theory. Frequency shifts and collision diameters determined from selective 
reflection data agree fairly well with values measured in absorption by other 
workers. 


INTRODUCTION 


Selective reflection of light from high pressure mercury vapor in the neigh- 
borhood of the resonance line, 2537 A, has been studied by Galt and Welsh 
(1957). It was shown that the reflection contours could be interpreted in 
terms of the classical theory of reflection from an absorbing medium. In 
particular, the damping constant, y, was found to be proportional to the 
number of atoms, N,, per cubic centimeter; however, the oscillator strength, 
f, which had to be assigned to the line in the high pressure vapor, is markedly 
smaller than the accepted value for the free atom. The variation of the damp- 
ing constant, y « Ng, in selective reflection is in accordance with most theore- 
tical treatments of collision damping in emission and absorption. However, 
no comparison of the selective reflection data with absorption data is possible 
because the high absorption index of the resonance line makes quantitative 
absorption measurements difficult even at low pressures. On the other hand 
the influence of foreign gases on the absorption line can be studied even at 
high pressures. Thus the absorption of mercury vapor at 2537 A was studied 
by Fiichtbauer, Joos, and Dinkelacker (1923) in the presence of foreign 
gases at pressures up to 50 atm. More recently, Robin and Robin (1951) 
and Robin (i954) have investigated the same absorption line with foreign 
gas pressures as high as 1500 atm. A direct comparison of the effects of foreign 
gases in absorption and in selective reflection is therefore possible. 

The influence of foreign gases on the selective reflection from mercury 
vapor at 2537 A was studied by Schnettler (1930), who measured the in- 
tensity of the resonance line from a mercury arc after it had been reflected 
from a quartz— mercury vapor interface. Since only a slight reduction in 
intensity was observed for foreign gas pressures up to 15 atm., Schnettler 
concluded that, unlike the case of resonance radiation, the mean lifetime of 
the excited atom plays no role in selective reflection; this conclusion was 
later confirmed more directly by Hansen and Webb (1947). The slight de- 
crease in intensity at high pressures observed by Schnettler could be inter- 
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preted as a diminution of the absorption index at the resonance frequency by 
collision broadening of the line. Since nearly monochromatic light was used 
for observing the selective reflection, no measurement of the damping constant 
was possible. 

In the present investigation, light of continuous spectral distribution was 
used to study the effect of foreign gases at pressures up to 1500 atm. on the 
selective reflection at 2537 A. With an adaptation of the apparatus used 
for pure mercury vapor at high pressures, the intensity contour of the selective 
reflection was studied as a function of foreign gas pressure. 


EXPERIMENTAL PROCEDURE 


The reflection cell, the detail of which is shown in Fig. 1, was a modification 
of the cell used for pure mercury vapor, and was operated in the same pressure 


Oh 


y 
y 
y 
Z 


SSO 
aN 


Veda 


Fic. 1. Reflection cell for the study of selective reflection from mixtures of mercury vapor 
and foreign gases. A—fused quartz cone; B—Teflon packing; C—fused quartz shell; B-steel 
tube; E—Alundum insulating grain; F—porcelain insulator for thermocouple wires; G—- 
liquid mercury; H—Teflon seat for C; J—platinum wire furnace; J—supporting tube for A. 
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vessel with the same optical arrangement (Galt and Welsh 1957). The shell, 
C, surrounding the conical fused quartz window was made of fused quartz 
rather than of steel since, in this case, it need not be impermeable to gas. 
Liquid mercury filled the lower end of the shell to within one millimeter of 
the reflecting surface. A channel formed by a narrow flat strip ground on the 
solid quartz cone allowed gas to enter the space between the window and the 
mercury meniscus to mix with the mercury vapor. The mercury and the 
lower end of the conical window were heated to temperatures between 650° C. 
and 750° C. by a coil of platinum wire embedded in refractory cement. To 
measure the temperature of the mercury pool, a Pt: Pt-Rh thermocouple 
was inserted into a small re-entrant tube in the bottom of the fused quartz 
shell. The space surrounding the outer cone was filled with Alundum grain 
for thermal insulation. The solid cone was fastened to a metal sleeve, J, and 
could be easily withdrawn to replenish the mercury without disturbing the 
rest of the apparatus. This sleeve also served as a labyrinth to prevent mercury 
vapor diffusing into the light path between the quartz cone and the pressure 
window. 

The thermocouple wires entered the pressure vessel through two glands 
of the type shown in Fig. 2. The Teflon packing, B, forms a pressure seal 


Fic. 2. High pressure thermocouple lead-in. A—platinum wire; B—Teflon packing; 
C—steel capillary tube hard soldered to A; D—collar soldered to C; E—conventional cone- 
type pressure seal; F—wall of pressure vessel. 


and also provides electrical insulation for the wire, A. The lead should, of 
course, be continuous in passing through the pressure wall to avoid a spurious 
e.m.f. in the thermocouple circuit. However, platinum wire by itself is too 
soft to withstand the pinching-off effect when the packing is compressed. It 
was found that a satisfactory seal could be made by soldering the platinum 
wire into a piece of fine steel capillary, C. To give added support a small 
steel collar, D, was soldered onto the steel capillary. Electrical lead-ins em- 
ploying Teflon packing have been used frequently in this laboratory, and 
appear to be quite reliable even at a pressure of 5000 atm. 

In contrast to the experiments with pure mercury vapor at high pressures, 
the foreign gas experiments were relatively easy to perform. The reflectivity 
curves were more self-consistent than those for pure mercury vapor because 
the complete sequence of curves was obtained without dismantling the 
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apparatus. Also, because of the lower temperatures used, there was no deterior- 
ation of the reflecting surface due to recrystallization of the silica. During an 
experiment the temperature of the liquid mercury was held constant to 
give a vapor pressure in the range 40 to 70 atm., and the total pressure applied 
by the foreign gas, as measured on a Bourdon gauge on the gas compressor, 
was varied between 100 and 1500 atm. The partial pressure of the mercury 
was taken to be the vapor pressure, calculated from the data of Bernhardt 
(1925), at the temperature indicated by the thermocouple, and any change 
in vapor pressure caused by the presence of the gas was neglected. The tem- 
perature shown by the thermocouple was used to calculate also the foreign 
gas densities, although radiation pyrometer measurements indicated that 
the temperature of the reflecting surface was somewhat higher. The densities 
for hydrogen were obtained by extrapolating the isothermal data of Michels 
and Goudeket (1941), and for helium by extrapolating the straight line 
relationship for the temperature variation of the second virial coefficient 
given by Schneider and Duffie (1949). The estimated high temperature virial 
coefficients given in Circular 564 of the National Bureau of Standards (1955) 
were used for nitrogen and argon. 

Intensity traces of the selective reflection were obtained with the recording 
spectrograph and were reduced by a method similar to that used for pure 
mercury vapor (Galt and Welsh 1957). The determination of the absolute 
reflectivity was complicated by the presence of a gas with an appreciable 
index of refraction in contact with the window. The contributions to the 
refractive index, m,—1, from the foreign gas, and ,—1, from the mercury 
vapor, were calculated separately from the partial densities and added. At 
the highest pressure the presence of the gas decreased the background re- 
flectivity, Ro, from 4% to 3% for argon and nitrogen; for hydrogen and 
helium the decrease was much smaller. Most of the experiments were carried 
out with the liquid mercury at a temperature of 685° C., for which the relative 
density, N,/No (No = Loschmidt’s number), of the vapor was 12. The 
relative densities of the foreign gas, N,/No, ranged from 15 to 180. 


RESULTS AND DISCUSSION 


Examples of the changes in the selective reflection from mercury vapor pro- 
duced by the addition of foreign gases at high pressures are shown in Figs. 3 
and 4. The predominant effect is a lowering of the maximum of the reflection 
contour, which tends to obliterate the pattern at the highest pressures. A 
weak subsidiary maximum is observed at 2540 A on the curves for the lower 
pressures; this maximum has been previously interpreted in terms of Hg» 
molecules (Galt and Welsh 1957). 

The main features of the experimental reflection contours can be explained 
satisfactorily by the classical theory of anomalous dispersion. With the 
facilities of the Computation Centre, University of Toronto, theoretical 
curves were calculated to fit the experimental curves using a procedure 
similar to that employed for pure mercury vapor. For the oscillator strength, 
f, the value 0.0267 for the free atom was used, and the number of atoms per 
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Fic. 5. Comparison of calculated (dashed curve) and experimental reflection curves for 
mixtures of mercury vapor and helium. 
cubic centimeter, N,, was calculated as outlined by Galt and Welsh (1957). 
The damping constant, y, was determined from the best agreement of the 
calculated reflection curves with the experimental curves, using a trial and 
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error procedure. The calculated curves for the lowest pressures were modified 
slightly to correct for the instrument. line shape of the spectrograph before 
they were compared with the experimental curves. Fig. 5 shows some of the 
theoretical curves for helium as foreign gas and the experimental curves to 
which they correspond. The fit is poor on the low frequency side because of 
the Hg» reflection which could not be allowed for in the calculations. The 
agreement between the calculated and experimental curves for the other 
foreign gases was similar to that shown for helium. In most cases the cal- 
culations were confined to the lower part of the pressure range since the 
reflection patterns at the highest pressures were comparatively poorly defined 
(Figs. 3 and 4). 

The values of y obtained from the curve-fitting procedure are plotted in 
Fig. 6 as a function of the relative density, N,/No, of the foreign gas. In- 
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Fic. 6. Variation of the damping constant, y, with density for the various foreign gases. 


cluded in the plots are the values of y for N, = 0, calculated from the damping 
coefficient, y/N, = 4.8 X 10-* cm.’ sec.—!, for pure mercury vapor (Galt and 
Welsh 1957). In all cases y increases linearly with N, so that a damping co- 
efficient, y/N,, can be determined for each foreign gas; these are listed in 
Table I. 

The optical collision diameter, p, which can be considered as an indication 
of the average extent of the intermolecular forces, can be derived from the 
slope of the lines in Fig. 6 by means of the relation 





Ss eS 
P4N,V 2xkT’ 


where yu is the reduced mass of the mercury atom and the foreign gas mole- 
cule (Weisskopf 1933). In Table I the values of p obtained from the selective 
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TABLE I 
DAMPING COEFFICIENTS AND COLLISION DIAMETERS FOR Hg 2537 








: Selective reflection Absorption 
Foreign gas Ce oe 
Damping coefficient, Collision diameter, Cee diameter, 


+/ Ng, cm. sec.~! X 10° 


H: 4.47 (685° C.) { ~ 

He 3.09 (685° C.) 5.4t 

N; 2.11 (685° C:) 8.1* 
1.57 (685° C.) 

A { 1:35 (ras C3 { 9.4° 


*Calculated from the data of Fiichtbauer, Joos, and Dinkelacker (1923). 
fCalculated from the data of Robin (1954). 


reflection measurements are compared with values calculated from the 
absorption data of Fiichtbauer, Joos, and Dinkelacker (1923) and of Robin 
(1954). The agreement, although not good, can be considered satisfactory 
evidence that damping in selective reflection is produced by the same collision 
processes that are responsible for the broadening of the absorption line. 

One experiment was carried out at a higher temperature, 748° C., with 
argon as the foreign gas. As is evident from Fig. 6, the damping coefficient 
and hence the collision diameter (Table I) is smaller at the higher tempera- 
ture. This temperature effect is not what one would expect and may perhaps 
be spurious. The matter could not be investigated further because the experi- 
ments were abruptly terminated by the bursting of the pressure window 
and the destruction of essential parts of the apparatus. 

As the pressure of the foreign gas increases, the peak of the reflection 
shifts towards higher frequencies for hydrogen and helium (Fig. 3) and 
towards lower frequencies for nitrogen and argon (Fig. 4). However, when 
theoretical curves are fitted to the experimental contours, it is seen that the 
displacement of the origin of the reflection pattern, »’, from the frequency 
of the unperturbed resonance line, », is always towards lower frequencies. 
For nitrogen and argon the shift, v»—»’, is a linear function of the gas density 
(Fig. 7). For hydrogen and helium the shifts appear to increase more rapidly 
than linearly at the higher densities. The shifts for hydrogen, argon, and 
nitrogen are in the same direction and of the same order of magnitude as the 
frequency shifts of the maximum of the absorption line, as measured by 
Fiichtbauer, Joos, and Dinkelacker (1923). The shift observed with helium 
as foreign gas is in the direction opposite to that observed by Robin (1954) 
for the absorption line; however, the absorption measurements may be 
complicated by the presence of HgHe molecules, which do not seem to be 
present to any extent at the high temperature used in the reflection experi- 
ments. 

The effect of foreign gases at high pressures on the selective reflection from 
mercury can thus be interpreted in a straightforward way in terms of collision 
damping. The main difference between the results with foreign gases and 
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Fic. 7. Shift of the resonance frequency, vo — yo’, as a function of gas density for the various 
foreign gases. 


the results with pure mercury vapor at high pressures is that, for the former, 
the oscillator strength is that of the free atom while, for the latter, an oscillator 


strength decreasing with increasing mercury vapor pressure must be assumed. 
This effect in pure mercury vapor must be a consequence of coupling of the 
mercury atoms, and it appears that the coupling is not present when foreign 
gases at 100 atm. or more are mixed with the vapor. The range of foreign 
gas pressures up to 100 atm., which could not be investigated in the present 
apparatus, should show effects indicative of the decoupling of the mercury 
atoms with increasing foreign gas pressure. 
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NOTES 


A NOTE ON A PROOF BY DYSON 


R. O. A. ROBINSON 


In the second edition of his Cornell notes Dyson (1952) gives a simplified 
proof ci the important theorem regarding the expansion in normal form of 
any simple product O of electron—positron and electromagnetic field operators. 
The object of this note is to point out an omission from the proof and to 
rectify it. 

The theorem is stated as follows: ‘‘We define a factor-pairing of O by picking 
out from O a certain even number of factors, either all or none of any inter- 
mediate number, and associating them together in pairs. To each factor- 
pairing ” corresponds a normal constituent O, obtained as follows: for each 
pair of factors PQ which is paired in , O, contains the numerical factor (PQ) 
(the vacuum expectation value of the product PQ), the order of P and Q being 
maintained as it was in O. The remaining unpaired factors Ri, Ro, ... appear 
in O, rearranged in normal form, viz. 


(1) On = +(PQ)(P’'Q')... N(RiR2...), 


the sign in front being + or — according to the even or odd character of the 
permutation of the electron—positron operators from the order in which they 
are written in O to the order in which they are written in equation 1. With this 
definition of O, we have the following theorem. Every operator product O is 
identically equal to the sum of the O, obtained from all its factor-pairings.” 

Let O be a product of m factors. Since from the definition of the normal 
product the theorem is trivially true for m = 1 and m = 2, it is only necessary 
to prove it true for arbitrary m assuming it true for m—2. This is done in two 
steps. First Dyson shows that if O’ is an operator involving m—2 factors, the 
theorem is true for all 


(2) O = (PQ#¥QP)O’ 


where P and Q are field operators and where the plus sign appears only if P 
and Q are both electron—positron operators. Further, by using the commutation 
rules it is always possible to write 


(3) O = N(O)+2 


where > is a sum of terms of type 2. The theorem is trivially true for N(O). 
Since it is also true for ©, the conclusion is made that it is true for O. 

However, in order to make this complete it must be shown that expansions 
in accordance with the theorem are additive, i.e. Expansion of A + Expansion 
of B = Expansion of (A + B), A and B being any two products. For example, 
let P:iQi, P2Q2, and PQ be products of two fermion operators satisfying the 
relation 
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(4) P,Q: + P2Q2 = PQ 


and let O be any product. Expanding the products, PQO, and cancelling the 
appropriate terms, we must verify that 


(5) (P,L)N(Qi, O—L)—(QiL)N(Pi1, O—-L)+(PiL)(QiM)N(O—L—M) 
+(P2L)N(Q2, O—L) —(QeL)N(P2, O—L)+(P2L)(Q2M)N(O-L—M) 
=(PL)N(Q, O—L)—(QL)N(P, O—L) £(PL(QM)N(O—L—M) 


where L, M are factors of O, and O—L, O—M, O—L—M represent the 
products obtained by omitting from O one or both factors. Both sides of 
equation 5 are to be summed over L and M. The + signs depend on the 
ordering of the fermion operators of O, but the point of importance is that, 
for any given ordering, they are either all + or all —. 

The complexity of equation 5 disappears if attention is restricted to the 
requirements of equation 3. We are concerned only with combinations effected 
by using the commutation relations. Thus factors are always combined in 
pairs as in 4 and the only case to be considered for fermion operators is: 


(6) P, = Q,, P, = Q,, P and Q both c-numbers. 


Subject to these conditions, the right-hand side of equation 5 vanishes, since 
the vacuum expectation value of single operators is always zero. Substituting 
for P2 and Q2 and using the relation 
(7) Do (QL)PiM)N(O-L—M) = — 2) (PiL)(QiM)N(O-L—M) 
L,M L.M 

the left-hand side also vanishes. The minus sign on the right-hand side of 
equation 7 results from the interchange in position of the fermion operators* 
L and M in corresponding normal constituents in accordance with equation 1. 

In the case where products of two boson operators are combined, we choose, 
analogously to equation 4, 


(8) P:Q:—P2Q2 = PQ. 

The conditions 6 remain unchanged while equation 5 becomes 

(9) (PiL)N(Qi, O—L)+(QiL)N(Pi, O—-L)+(PiL)(QM)N(O-L—M) 

—(P2L)N(Qo, O—L) —(Q2L)N(P2, O—L) —(P2L)(Q2M)N(O-—L—M) 

=(PL)N(Q, O—L)+(QL)N(P, O-—L)+(PL)(QM)N(O—L—M). 

In this case LZ and M must both be boson operators, since otherwise the vacuum 

expectation values vanish, and so the + signs do not arise. It is easy to see 

that equation 9 holds subject to the conditions 6. Finally, the proof remains 

essentially unchanged if instead of expanding products such as PQO we con- 


sider more general products of the form OPQO’. This completes the theorem 
as required. 


*If L and M are not both fermion operators the corresponding terms in 7 vanish on account 
of the vacuum expectation values. 
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PREPARATION OF PURE LITHIUM SAMPLES 
FOR ELECTRICAL MEASUREMENTS* 


W. B. PEARSON 


The normal method of preparing samples of the alkali metals for electrical 
measurements is to run the molten metal under vacuum into a soft-glass mold 
which has platinum leads sealed in at the required positions. It has, however, 
not been possible to use this method for preparing lithium samples because 
of the rapid reaction of molten lithium with the glass. Some attention has 
therefore been given to making special non-reactive glasses, and MacDonald 
and Stanworth (1950) have for instance shown that it is possible to prepare 
pure lithium samples in a special borate-clad glass, C 10. 

We now find that pure lithium samples can be prepared in soft-glass molds 
by the following method. The lithium is first immersed in mineral oil (Stanolax) 
and cut to the appropriate size but then, instead of the Stanolax being washed 
off with benzene in the normal manner, it is rapidly mopped off the surface 
with cloth or absorbent paper before the metal is placed in the melting tube. 
The film of oil remaining on the surface of the lithium is sufficient to cover the 
surface of the melting tube and mold and prevent reaction, provided that the 
lithium is driven down into the mold with a small pressure of helium gas as 
soon as it has melted, and also that the specimen tube is immediately removed 
from the hot oil bath and placed in a cold oil bath. This can be done quite 
conveniently and a pure specimen with a clean shiny surface is obtained. If, 
however, the metal is not quenched, it starts to react with the glass in a few 
seconds and the shiny metal-glass interface turns dark brown. 

The residual resistance ratio of a specimen prepared in this manner for Hall 
effect measurements was 1.9X10-*, compared with 1.5 10-* for an extruded 
wire of lithium of rather similar initial purity. 


MacDonaL.p, D. K. C. and STaNwortTH, J. E. 1950. Proc. Phys. Soc. B, 63, 455. 
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HEIGHT-GAIN IN THE FORWARD-SCATTERING OF 
RADIO WAVES BY METEOR TRAILS* 


C. O. HINEs AND M. O’Grapy 


Radio signals in the V.H.F. band can be transmitted over long paths by 
means of forward-scattering from meteor trails (Forsyth and Vogan 1955). 
Even if the parent meteors were incident on the atmosphere uniformly and 
isotropically, the trails which are effective in the scattering process would not 
be distributed uniformly. Instead, various observational effects tend to empha- 
size the contribution from certain regions of the sky and suppress that from 
others. The factors involved are fairly well understood, and contour charts 
of the relevant distribution functions have been presented (Eshleman and 
Manning 1954; Hines and Pugh 1956). In particular, regions somewhat to 
either side of the great circle path tend to dominate the distribution. 

The observed occurrence rate of useful meteors and the mean received 
signal level will, of course, depend on the particular antenna systems employed. 
The dependence can be predicted -theoretically by combining the antenna 
illumination pattern at the meteoric height § with the pertinent distribution 
function, and then integrating the product, for each assumed antenna system 
in turn. 

Typical results of this process have been presented (Hines, Forsyth, Vogan, 
and Pugh 1955), to illustrate the effect of various changes in the azimuthal 
gain pattern. In those calculations, simplified single-lobed vertical gain 
patterns were adopted for both transmitter and receiver, with elevations such 
that the maxima of the two antenna beams intersected at the § level, and § 
was taken to be 100 km. 

It is by no means evident that the optimum elevation angle for meteoric 
scattering is the one previously adopted. Moreover, while 100 km. is a typical 
height for meteors, it is by no means exclusive; fairly wide variations must be 
anticipated, particularly towards lower values. At the suggestion of Dr. P. A. 
Forsyth, a new set of computations was undertaken to investigate the effect 
of variations in elevation angle and in §. The results of the computations are 
summarized here. 

A relatively refined antenna pattern was assumed in this work, where the 
amplitude gain function was of the form 


sisl(iak/s) cin o] sin[{(wa/) sin e] sin[(4b/A) cos e sin a] 


sin ée cos e sin a 


This is the Fraunhofer form appropriate to an effective aperture of vertical 
dimension a and horizontal dimension }, radiating at wavelength A, and 
situated at a height h above a perfectly conducting ground. The normal to the 
aperture lies in the great circle direction towards the opposite end of the trans- 
mission path; ¢ measures elevation above the local horizontal plane, and a 


*The work was performed under project PCC. No. D48-95-11-01. 
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measures azimuth from the great circle direction. The antennas contemplated 
were five-element Yagis mounted horizontally, with vertical beam width 
(between half-power points) of 64° and horizontal beam width of 50°, corre- 
sponding to a = 0.84 and b = 1.05\. For the range of h’s considered, the 
first factor in the gain expression gives rise to the dominant variation with 
elevation angle. It results, of course, from the reflective interference introduced 
by the ground, and more than one of its lobes must be taken into account in 
obtaining a proper estimate of height-gain in meteoric scattering. 

The values § = 80, 90, 100, and 110 km. were adopted in succession. In 
each case, the azimuths and elevations were determined for a large number of 
points at the © level, and the corresponding antenna gains were calculated 
both for the transmitting and for the receiving antenna, assuming a trans- 
mitter—-receiver separation of 1000 km. This was repeated for several values of 
h in turn, to give the effect of changing antenna heights. The resultant set of 
illumination patterns was then combined, each pattern in turn, with the 
contour charts giving the distribution of observable meteors and of mean 
signal sources, and the products in turn were integrated over the © level. 
(The charts employed, those drawn as solid lines by Hines and Pugh, were 
based on an assumed § of 100 km.; but they are not expected to be as strongly 
dependent on § as are the gain functions, so their continued use seems 
justified.) 

The results of this synthesis are displayed by the solid-line curves in Fig. 1. 
The first column of graphs indicates the variation, with antenna height and 
with meteoric height, in the number of meteors detected with received signals 
exceeding some arbitrary counting level. The second depicts the corresponding 
variations in mean signal level, or in total duration above a counting level. 
The scales of ordinates in the two columns are arbitrary, but they are con- 
sistent from one row to the next. The heights, /o, at which the antennas would 
have to be placed to get maxima intersecting at the appropriate § level are 
indicated by arrowheads on the horizontal axis. 

On the whole, the graphs speak for themselves, but a few points may be 
emphasized. The dominant feature of each curve is the hump which appears, 
centered on a height near ip, and attributable to a degree of intermeshing of 
the two main antenna beams at the meteoric level. As / increases beyond ho, 
subsidiary antenna lobes play a more important part, and the fall-off is not 
as rapid as might have been expected; in fact, a secondary maximum arises 
for large $. This behavior suggests that the form of the curves is dominated 
by the variations in antenna illumination alone. The modifying effects of the 
distribution functions are appreciable, however, as is indicated by the dif- 
ferences between corresponding members of thé two columns. 

In addition to the symmetrical situation so far considered, with the two 
antennas at the same height, calculations were made on the assumption of 
different antenna heights. Typical curves are shown for the © = 90 and 
100 km. levels. These obtain when one antenna is maintained at a height of 
1.3 (dotted curves) or 2.2\ (broken curves), the other antenna being allowed 
to vary in height. 

Meteor trails are produced, on the average, at greater heights in the morning 
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Fic. 1. Height-gain in meteoric forward-scattering. The ordinates indicate, on arbitrary 
scales, (a) the detection rate of meteors above a counting level and (6) the mean signal strength 
obtained by meteoric scattering, for various heights ) of the meteor level and for a range of 


heights h of the antennas. 

hours and at lesser heights in the afternoon hours. Accordingly, with fixed 
heights for the antennas, a diurnal variation will be introduced in the apparent 
meteoric rate and mean signal strength as a result of § changes alone, quite 
apart from actual variations in the absolute meteoric rate. In this connection, 
it should be noted that the mean signal level is inversely proportional to the 
ambipolar diffusion coefficient, which may be expected to increase by a factor 
of four as the § level increases by 10 km.; this change has not been incorporated 


in the graphs. 
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